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Abstract 

This  study  produced  an  image  processing  algorithm 
development  tool  on  a  VAX  computer  system  using  the  recent 
advances  in  an  Image  Algebra  developed  by  G.  Ritter  et  al  at 
the  University  of  Florida.  The  image  algebra  provides  the 
basis  for  a  hardware  and  software  independent  environment  for 
the  expression  of  practically  all  image  processing  algorithms. 

The  goals  of  this  project  were  twofold.  The  first  goal 
was  the  implementation  of  the  image  algebra  operators  in  a 
high  level  language  to  achieve  hardware  independence.  The 
second  goal  was  the  design  and  implementation  of  a  flexible 
preprocessor  that  could  translate  image  processing  algorithms, 
written  in  the  image  algebra  language,  into  a  high  level 
computer  language  which  could  be  compiled  and  executed  on  the 
VAX  computer. 

The  implementation  was  achieved  in  the  PASCAL  computer 
language.  All  of  the  basic  image  algebra  operators  and  the 
preprocessor  were  successfully  programmed  on  the  VAX  computer, 
but  a  complete  software  independence  of  the  image  algebra  was 
not  achieved.  This  version  also  produces  very  large  blocks  of 
executable  code  for  relatively  simple  algorithms. 

Examples  of  the  power  and  simplicity  of  the  image  algebra 
langauge  and  preprocessor  environment  are  included. 
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A  PASCAL  IMPLEMENTATION  OF  THE  IMAGE  ALGEBRA 


I .  Introduction 


Background 

Developers  of  image  processing  algorithms  in  military, 
industrial,  and  academic  organizations  have  built  large 
computer  programs  to  facilitate  their  development  work.  Each 
program  is  the  result  of  a  constant  evolution  of  image 
processing  operations  within  each  organization.  Unfortunately, 
in  the  process  of  building  excellent  "in-house"  development 
tools,  each  organization  has  created  an  image  processing 
environment  different  from  every  other.  Some  image  processing 
operations  are  unique  to  one  environment  while  others  may  be 
common  to  a  number  of  environments.  But,  since  there  is  no 
common  mathematical  basis  for  all  of  these  operations,  it  is 
often  impossible  to  analyze  and  compare  one  environment  with 
another . 

This  situation  creates  two  major  problems  for  the 
customers  of  image  processing  software,  especially  for  the 
Department  of  Defense:  (1)  increased  costs  through  funding  of 
research  and  development  in  several  organizations  separately 
producing  similar  or  identical  processes,  and  (2)  difficulty 
in  analyzing  and  comparing  the  performance  of  competing 
algorithms  (1:1-2). 
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Realizing  this  situation  would  create  incompatibility  and 

unnecessary  confusion  between  competing  automatic  target 

recognition  (ATR)  algorithms,  the  Air  Force  Armament 

Laboratory  (AFATL)  at  Eglin  AFB ,  Florida  contracted  the 

University  of  Florida  to  develop  an  Image  Algebra  (IA).  The 

goal  the  Image  Algebra  Project  is 

...the  development  of  a  complete,  unified  algebraic 
structure  that  provides  a  common  mathematical  envi¬ 
ronment  for  image  algorithm  development,  optimization, 
comparison,  coding,  and  performance  evaluation.  (1:2) 

The  primary  goal  of  the  AFATL  was  a  concise  mathematical 

structure  for  representing  image  transformations  that  was 

simple  to  learn  yet  powerful  enough  to  perform  any  image 

processing  operation.  Toward  this  goal,  the  AFATL  outlined  a 

number  of  desirable  properties  that  a  useful  image  algebra 

should  possess.  These  properties  guided  the  development  of 

the  image  algebra  (1:3-4): 

1.  Basic  operations  that  can  model  any  image. 

2.  Elemental  operations  that  can  be  combined  to 
express  any  gray  level  image  transformation 

3.  Elemental  operations  that  are  few,  simple,  and 
easily  learned  by  potential  users. 

4.  Theorems  that  enable  the  simplification  and 
optimization  of  algorithms  through  the  use  of 
identities  involving  the  operators. 

5.  Notation  that  provides  an  understanding  of  image 
manipulations  and  is  capable  of  suggesting  new 
techniques . 

6.  Notation  that  allows  programming  languages  to 
substitute  concise  image  algebraic  expressions  for 
large  blocks  of  code. 

7.  Notation  that  allows  the  use  of  libraries  of 
image  transformations. 

8.  Machine  and  language  independence. 

9.  Compatibility  with  both  sequential  and  massively 
parallel  processors. 
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The  result  of  this  work  by  the  University  of  Florida  is 
an  algebraic  structure  consisting  of  two  operands,  images  and 
templates,  and  eight  binary  operators.  Operations  may  occur 
between  two  images,  an  image  and  a  template,  or  two  templates. 
The  result  of  these  operations  may  be  '■n  image,  a  template,  or 
a  scalar. 

Problem 

The  task  of  this  thesis  project  is  to  implement  the  image 
algebra,  as  conceived  by  the  University  of  Florida,  on  the  VAX 
11/780  located  at  the  Air  Force  Institute  of  Technology.  The 
implementation  is  to  preserve  the  desirable  properties  of  the 
image  algebra  outlined  by  AFATL  and  listed  above.  Consistency 
between  the  Air  Force  Institute  of  Technology  image  algebra 
(AFITIA)  and  the  University  of  Florida  image  algebra  (UFIA)  is 
to  be  maintained. 

Scope 

The  description  of  the  image  algebra  continues  to  evolve 
at  the  University  of  Florida.  Consequently,  this  project  will 
not  reflect  the  more  recent  changes  to  the  structure.  This 
imi  ementation  is  based  upon  the  description  of  the  image 
algebra  in  the  Image  Algebra  Tutorial  and  in  the  Final  Report 
from  Phase  I  of  the  Image  Algebra  Project  (1,2). 


of  the  author.  All  of  the  PASCAL  source  code  listed  in  the 
appendices  was  originally  developed  in  TurboPascal  (3), 
translated  to  VAX  PASCAL  (4),  and  ported  to  the  VAX  11/780  at 
the  Air  Force  Institute  of  Technology  (AFIT). 


General  Approach 

The  general  approach  taken  in  this  thesis  project  is  the 
development  of  a  computer  program  from  the  bottom  up.  This  is 
contrary  to  accepted  program  development  methods  but  it  is 
desirable  in  this  case.  Consistency  between  the  operands  and 
operators  of  the  UFIA  description  and  AFITIA  implementation  is 
a  primary  concern  in  this  project.  Consequently,  an 
implementation  driven  by  the  structure  of  the  low  level 
operands  and  operators  is  preferable. 

The  first  task  involves  implementing  the  image  algebra 
(IA)  operands  in  PASCAL  data  structures  and  translating  the  IA 
operators  to  PASCAL  procedures  and  functions.  Based  upon  this 
design,  the  next  task  is  the  programming  of  an  IA  preprocessor 
that  translates  an  IA  description  of  any  image  processing 
operation  into  PASCAL  source  code  for  subsequent  compilation. 
The  final  task  is  to  install  the  image  algebra  routines  and 
preprocessor  into  an  environment  on  the  VAX  that  allows  a  user 
to  quickly  implement  any  image  processing  operation. 
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Sequence  of  Presentation 

The  following  report  is  divided  into  four  major  sections. 
Section  II  defines  the  basic  operands  and  operators  of  the 
image  algebra.  Section  III  discusses  the  design  of  a  general 
implementation  of  the  image  algebra  to  achieve  a  useful  image 
processing  development  tool.  Section  IV  covers  the  major  part 
of  this  project:  the  design  and  implementation  of  an  image 
algebra  at  AFIT.  Section  V  compares  three  different  image 
processing  algorithms  implemented  in  both  PASCAL  source  code 
and  the  image  algebra  language.  Finally,  Section  VI  contains 
some  general  remarks  on  the  accomplishments  of  this  project 
and  suggests  some  follow  up  projects.  The  appendices  contain 
the  source  code  to  all  of  the  AFIT  image  algebra  programs. 
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After  studying  the  underlying  operations  of  hundreds  of 


existing  image  processing  algorithms,  the  University  of  j 

Florida  investigators  determined  all  of  these  routines  could  j 

t 

be  performed  with  a  small  set  of  operands  and  operators.  In 
fact,  they  proved  this  set  of  operands  and  operators,  in 
conjunction  with  a  programming  language  such  as  FORTRAN,  is 
sufficient  to  perform  any  image-to-image  transformation  (1:48— 

62).  This  implies  the  image  algebra  operands  and  operators 
can  be  used  to  program  all  image  transformations. 

The  next  sections  describe  the  image  algebra  operands, 
operators,  and  notation  as  developed  by  the  University  of 
Florida.  These  sections  provide  the  reader  with  enough 
background  to  understand  the  subsequent  sections  of  this 
report.  The  reader  is  directed  to  the  Image  Algebra  Tutorial 
for  a  more  complete  mathematical  description  of  the  operands 
and  operators  (1:9-39). 

Image  Algebra  Operands 

The  UFIA  description  contains  six  explicit  operands:  the 
set  of  real  numbers,  the  set  of  complex  numbers,  any  finite 
subset  of  k-dimens ional  space  called  X,  the  power  set  on  X, 
the  set  of  all  images  on  X,  and  the  set  of  templates  on  X. 

(1:6)  Although  these  six  operands  completely  define  a 
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heterogeneous  algebra,  there  are  only  two  operands  of  primary 
interest  in  this  project:  images  and  templates. 

Images  are  multiple  dimension  arrays  of  arbitrary 
integral  size  and  may  be  integer,  real,  or  complex  valued. 

The  points  in  images  are  referred  to  as  pixels  and  each  pixel 
in  the  array  has  a  unique  spatial  identifier.  In  common  image 
processing  tasks,  images  represent  two  dimensional  arrays  of 
gray  values  describing  scenes  collected  from  vidicons,  laser 
rangefinders,  synthetic  aperture  radar,  infrared  imagers,  or 
other  sensing  devices. 

Images  in  the  UFIA  are  denoted  with  capital  letters  from 
the  beginning  of  the  alphabet  such  as  A,  B,  or  C.  Arbitrary 
pixels  in  n-dimensional  images  are  depicted  with  n-dimensional 
vectors  such  as  k.  That  is,  an  arbitrary  pixel  in  image  A  is 
denoted  A(k) .  The  gray  level  at  a  particular  image  pixel  is 
depicted  with  a  lower  case  letter  of  the  image  designator: 
a(k)  represents  the  gray  level  at  pixel  k  in  image  A  (1:9). 

To  perform  many  image  processing  tasks,  it  must  be 
possible  to  selectively  choose  and  weight  image  pixels  within 
small  neighborhoods.  The  image  algebra  provides  templates  for 
this  job.  Templates,  sometimes  called  windows  or  masks,  are 
multiple  dimension  arrays  of  integral  size  and  may  be  integer, 
real,  or  complex  valued.  Templates  may  be  of  arbitrary  size 
and  configuration  but  they  are  usually  much  smaller  than 
images.  Image  algebra  templates  consist  of  two  elements: 
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Figure  2.1  :  Examples  of  Image  Algebra  Templates 


(1)  a  configuration  of  pixels,  and  (2)  a  weight  at  each  pixel 
in  the  configuration.  The  template  configuration  describes 
the  relative  orientation  of  all  template  pixels  where  the 
center  pixel  is  commonly  defined  as  the  origin.  The  center 
pixel  can  be  defined  as  any  pixel  in  the  configuration:  it  is 
not  necessarily  the  physical  center  of  the  configuration. 

Only  pixels  in  the  template  configuration  may  have  non-zero 
values;  any  pixel  not  in  the  configuration  has  zero  weight 
(1:23).  Some  examples  of  two  dimensional  templates  are  shown 
in  Figure  2.1.  The  center  of  each  template  is  marked  with 
asterisks.  For  the  purposes  of  the  UFIA  and  this  report,  some 
of  the  template  configurations  have  special  names:  Figure 
2.1(a)  is  a  horizontal  template  of  size  five,  2.1(d)  is  a 
vertical  template  of  size  five,  2.1(b)  is  a  three-by- three 
Moore  template  (2:152),  and  2.1(c)  is  a  Von  Neumann  template 
of  radius  one  (2:152). 


The  UFIA  template  notation  is  similar  to  the  UFIA  image 


notation.  Templates  are  denoted  with  capital  letters  from  the 
end  of  the  alphabet  such  as  R,  S,  or  T.  Arbitrary  pixels  in 
m-dimensional  templates  are  depicted  with  m-dimensional 
vectors  such  as  y.  For  example,  arbitrary  pixel  locations  in 
template  S  are  denoted  by  S(y).  The  weight  of  a  template 
pixel  is  represented  by  a  lower  case  letter  of  the  template 
designator:  s(y)  denotes  the  weight  at  pixel  y  in  template  S 
(1:23). 

There  are  two  types  of  templates  defined  by  the  UFIA: 
variant  and  invariant.  Variant  templates  may  alter  the 
weights  of  their  pixels  depending  upon  the  location  of  the 
template  in  an  image;  that  is,  the  pixel  weights  in  variant 
templates  are  functions  of  the  template  position  in  an  image. 
As  the  name  implies,  invariant  templates  are  not  allowed  to 
alter  their  pixel  weights  during  any  operation  with  an  image. 

Image  Algebra  Operations 

Due  to  various  combinations  of  the  two  operands  and  eight 
operators  of  the  image  algebra,  there  are  fourteen  elemental 
operations  as  described  by  the  University  of  Florida.  The 
result  of  these  operations  can  be  an  image,  a  template,  or  a 
scalar.  There  are  unary  and  binary  operations  on  both  images 
and  templates  and  binary  operations  between  images  and 
templates,  but  these  fourteen  operations  are  not  necessarily 
the  only  ones  supported  by  the  image  algebra.  According  to 


the  University  of  Florida, 

In  general,  all  elementary  functions  supported 
by  all  the  common  higher  level  languages  such  as 
FORTRAN,  PASCAL,  etc.  (i.e.  sine,  cosine,  logs,  expo¬ 
nential,  etc.),  are  accepted  operations  of  the  image 
algebra .  (5:15) 

Figure  2.2  shows  the  names,  notation,  and  mathematical 
description  of  the  fourteen  elemental  operations  in  the  image 
algebra.  The  vector  k  represents  arbitrary  pixels  in  an  n- 
dimensional  image,  and  the  vectors  x,  £,  and  z  represent 
arbitrary  pixels  in  their  associated  m-dimens i onal  templates. 
The  ordered  pair  (k,c(k))  represents  an  arbitrary  image  pixel 
and  its  associated  gray  value  in  the  image  C.  Similarly,  the 
ordered  pair  (z.,t(z))  represents  an  arbitrary  template  pixel 
and  its  associated  weight  in  template  T.  The  zero  vector,  0, 
denotes  the  center  pixel  of  a  template.  In  the  dot  product 
operation,  f  is  a  scalar.  The  union  of  two  images  (A  U  B)  or 
two  templates  (R  U  S)  is  defined  as  the  image  or  template 
resulting  from  the  union  of  all  of  the  image  or  template 
pixels  in  the  two  images  or  templates.  The  intersection  of 
two  images  (A  fl  B)  or  templates  (R  fl  S)  is  defined  as  the 
image  or  template  resulting  from  the  intersection  of  all  of 
the  image  or  template  pixels  in  the  two  images  or  templates. 

Image-Image  Operations.  There  are  five  elemental  binary 
operations  between  images  called  image-image  operations.  Four 
of  these  operations  (addition,  multiplication,  maximum,  and 
exponentiation)  result  in  another  image.  These  operations  are 
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Image-Image  Operations 


operation  notat i on  mathematical  description 

addition  C  =  A  +  B  s  {(k,c(k)>  :  c(k)  =  a(k)  +  b(k)  ) 

where  C  is  defined  for  A  U  B 

multiplication  C  =  A  *  B  =  {(k,c(k))  :  c(k)  =  a(k)  *  b(k)  ) 

where  C  is  defined  for  A  n  B 

maximum  C  =  A  v  B  =  {(k,c(k))  :  c(k)  =  max [ a ( k ) , b ( k ) ]  } 

where  C  is  defined  for  A  U  B 

exponentiation  C  =  A  **B  s  {(k,c<k))  :  c(k)  =  a(k}  **  b(k)  ) 

where  C  is  defined  for  A  0  B 

dot  product  f  -  A  *  B  =  {  f  :  f  =  £  a(k)b(k)  } 

where  C  is  defined  for  A  fl  B 


Image-Template  Operations 

operation  notation  mathematical  description 

circle-plus  C  =  A  *  R  =  {(k,c(k))  :  c(k)  =  £  a(k+x)*r(x)  } 

where  C  is  defined  for  A(k)  n  R(0) 
circle-maximum  C  =  A  <S>  R  s  {(k,c(k))  :  c(k)  =max  [  a(  k+x  )  *r  (  x  )  ]  } 

where  C  is  defined  for  A(k)  n  R(0) 
square-maximum  C  =  A  H  R  =  {(k,c(k))  :  c(k)  =max [ a ( k+x ) +r ( x ) ] } 

where  C  is  defined  for  A(k)  fl  R(0) 


Template-Template  Operations 

operation  notation  mathematical  description 

addition  T  =  R  +  S  s  {(z,t(z))  :  t(z)  =  r(z)  +  s(z);  } 

where  T  is  defined  for  R  U  S 

multiplication  T  =  R  *  S  =  {(z,t(z))  :  t(z)  =  r(z)  *  s(z)  } 

where  T  is  defined  for  R  fl  S 

maximum  T  =  R  v  S  =  {(z,t(z))  :  t(z)  =  max [ r < z )  , s  (  z ) ]  } 

where  T  is  defined  for  R  U  S 

circle-plus  T  =  R  ®  S  =  {(z,t(z))  :  t(z)  =  £  r(x)*s<£)  ) 

where  z^=x  +  x.  and  T  is  defined  for  R(0)  U  S(j^) 

circle-maximum  T  =  R  ®  S  =  {(z,t(z))  :  t(z)  =  max[r(x)*s(jr)  ]  } 

where  z=x+^  and  T  is  defined  for  R ( 0 )  U  S(^) 
square-maximum  T  =  R  0  S  =  { (z, t( z) )  :  t(z)  =  max[r(x)+s(i) ]  ) 

where  z=x+£  and  T  is  defined  for  R ( 0 )  U  S(£) 


Figure  2.2:  Image  Algebra  Elemental  Binary  Operations 
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Image  A  Image  B 


Figure  2.3:  Example  of  Image-Image  Operations 


performed  on  a  pixel-to-pixel  basis:  a  pixel  of  one  image  is 
added,  multiplied,  maximized,  or  exponentially  multiplied  with 
its  spatially  corresponding  pixel  in  another  image(l:12).  For 
example,  Figure  2.3  shows  two  images,  A  and  B,  of  dimension 
5x5.  If  addition  (C  =  A  +  B)  is  performed  on  the  two  images, 
the  resulting  image  is  the  sum  of  the  two  images  at  each  pixel 
location:  C  (  1 1 )  =  A  (  1 1  >  *  Bill),  C(12)  =  A(12)  +  B(12),  .... 

C(55)  =  A(55)  +  B(55).  Similarly,  the  image  resulting  from 
the  maximum  operation  (C  =  A  v  B)  on  these  two  images  would 
be:  C  <  1  1  )  =  max { A  < 1 1 )  , B( 1 1  )  }  ,  C ( 1 2 )  =  max { A ( 1 2 ) , B (1 2 ) )  ,  . .  .  , 

C(55)  =  max { A ( 5 5 ) , B ( 5 5 ) } .  Multiplication  and  exponentiation 
are  analogous  operations. 

The  fifth  elemental  binary  operation  is  the  dot  (inner) 
product.  The  result  of  this  operation  is  a  scalar,  and  it 
provides  the  image  algebra  with  a  method  of  mapping  images  to 
the  real  numbers.  Thi3  operation  is  a  sum  of  the  product  of 


spatially  corresponding  pixels  in  two  images.  Using  the  two 
images  from  Figure  2.3,  the  result  =  A(11)*B(11)  +  A(12)*B(12) 
+  ...  +  A( 55 ) *B( 55 ) . 

These  five  elemental  image  operations  can  be  used  to 
define  other  operations  between  images  such  is  subtraction, 
division,  and  minimum  (1:13). 

Image-Template  Operations.  There  are  three  elemental 
binary  operations  between  an  image  and  a  template  called 
image- template  operations:  circle-plus,  circle-maximum,  and 
square-maximum.  The  result  of  these  operations  is  an  image. 
Each  pixel  of  the  resultant  image  is  a  weighted  function  of 
the  original  image  pixel  and  its  neighbors  delimited  by  the 
template  configuration  centered  on  this  pixel  (1:28).  All 
template  pixels  that  lay  outside  of  the  image  boundaries  at 
any  image  pixel  are  ignored. 

The  circle-plus  operation  performs  a  convolution:  each 
pixel  is  the  sum  of  the  products  of  template  pixels  and  image 
pixels  delimited  by  the  template  configuration  centered  over 
the  image  pixel  (1:28).  For  example,  Figure  2.4  shows  an 
image  A  of  dimension  5x5  and  template  Y  of  dimension  3x3  with 
its  center  defined  as  y(5).  If  the  circle-plus  operation  is 
performed  between  this  image  and  template,  one  calculation  of 
the  operation  is  performed  on  the  image  pixel  at  row  3,  column 
4,  shown  as  pixel  A(34)  in  the  figure.  For  the  calculation  at 


this  pixel,  template  Y  is  overlaid  on  image  A  with  the  center 
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Figure  2.4:  Example  of  Image-Template  Operations 

of  the  template,  y(5),  centered  over  pixel  A(34).  The 
operation  computes  the  sum  of  the  products  of  each  template 
pixel  and  the  image  pixel  it  overlays.  In  this  case,  A(34)  = 
y ( 1 ) *A( 23 )  +  y ( 2 ) *A( 24 )  +  y(3)*A(25)  +  y(4)*A(33)+  y(5)*A(34) 

+  y ( 6 ) *A ( 35 )  +  y ( 7 ) * A ( 4  3 )  +  y<8)*A(44)  +  y(9)*A(45).  This 
calculation  is  repeated  for  each  image  pixel  to  complete  the 
entire  circle-plus  operation. 

With  the  circle-maximum  operator,  each  image  pixel  x  is 
replaced  with  the  maximum  product  of  image  pixel  gray  values 
within  the  neighborhood  of  x  (defined  by  the  template  centered 
over  x)  and  the  template  pixel  weights  that  overlay  each  image 
pixel.  With  the  square-maximum  operator,  each  image  pixel  x 
is  replaced  with  the  maximum  sum  of  image  pixel  gray  values 
within  the  neighborhood  of  x  (defined  by  the  template  centered 
over  x)  and  the  template  pixel  weights  that  overlay  each  image 
pixel .  (1:28) 
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The  circle-maximum  and  square-maximum  operators  can  both 
perform  gray  level  dilation  depending  upon  the  weights  of  the 


template.  By  choosing  a  template  of  unity  weight  for  the 
circle-maximum  operation  and  a  template  of  zero  weight  for  the 
square-maximum  operation,  a  gray  level  dilation  is  performed. 
Figure  2.5(c)  shows  the  result  of  dilating  the  image  in  2.5(a) 
with  the  square-maximum  operator  and  template  shown  in  2.5(b). 
As  an  example  of  the  intermediate  calculations  in  the  square 
maximum  operation,  the  template  is  shown  centered  over  two 
pixels  in  2.5(a).  The  result  of  the  square-maximum  operation 
for  these  two  pixels,  the  maximum  value  after  summing  each 
template  weight  with  the  image  pixel  it  overlays,  is  shown  in 
2.5(c)  . 

These  three  elemental  image- template  operations  can  be 
used  to  define  other  operations  between  images  and  templates 
such  as  circle-minimum  and  square-minimum. 

Template-Template  Operations.  There  are  six  elemental 
binary  operations  between  two  templates  called  template- 
template  operations:  addition,  multiplication,  maximum, 
circle-plus,  circle-maximum,  and  square-maximum.  These 
operations  are  similar  to  their  image-image  and  image- template 
counterparts;  the  major  difference  is  the  region  over  which 
the  operations  are  defined.  Addition  and  maximum  operations 
are  defined  over  the  union  of  the  two  template  configurations, 
but  multiplication  is  defined  only  over  the  intersection  of 
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Figure  2.6:  Template  Configuration  Resulting  from 
T  =  R(+)S,  R ( v ) S ,  or  R[v]S 


the  two  templates.  Circle-plus,  circle-maximum,  and  square- 
maximum  operations  are  defined  over  the  configuration  formed  by 
the  union  of  the  pixels  of  one  template  centered  over  all  of 
the  pixels  of  the  other  template  (1:30).  For  example,  Figure 
2.6  shows  a  Von  Neumann  template,  R,  with  its  center  defined  as 
r(3),  and  a  horizontal  template,  S,  with  its  center  defined  as 
s(3).  The  template  resulting  from  the  circle-plus,  circle- 
maximum,  and  square-maximum  operations  is  shown  as  template  T 
with  center  t(24).  Template  T  is  formed  by  the  union  of  the 
configurations  constructed  by  overlaying  the  center  pixel  of 
template  R  on  each  pixel  of  template  S:  T(x)  =  R(0)  U  S(y).  In 
this  example,  the  pixels  at  t(21),  t(12),  t(22),  t(32),  and 
t ( 2 3 )  are  the  result  of  pixels  r(l)  through  r(5)  when  template 
R  is  centered  at  pixel  s(l). 


The  computation  of  the  resultant  weights  in  template  T 


under  the  operations  of  addition,  multiplication,  and  maximum 


is  simple.  Since  template  addition  and  maximization  are 
performed  over  the  union  of  the  two  templates,  the  resulting 
template  is  the  addition  or  maximum  of  pixels  from  each 
template  with  the  same  configuration  coordinates.  Template 
multiplication  is  performed  over  the  intersection  of  the  two 
templates:  each  resulting  pixel  is  the  product  of  a  pixel  from 
each  template  with  the  same  configuration  coordinates  (1:25). 
Recall  that  template  pixels  outside  the  configuration  have 
weights  of  zero. 

The  computation  of  template  T  resulting  from  the  circle- 
plus,  circle-maximum,  and  square-maximum  operations  between 
templates  R  and  S  is  more  difficult.  Conceptually,  template  T 
is  the  result  of  two  operations  between  the  operand  templates  R 
and  S.  First,  the  template  T  configuration  is  computed  from 
the  union  of  template  R  centered  over  each  cell  of  template  S, 
as  discussed  above.  Second,  the  weight  of  each  cell  in  the 
configuration  of  template  T  is  computed.  For  convenience, 
designate  each  cell  of  templates  R,  S,  and  T  as  R(ri,rj), 
S(si,sj),  and  T(ti,tj)  respectively,  where  (ri,rj),  (si,sj), 
and  (ti,tj)  denote  the  cell’s  offset  from  the  center  of  the 
template.  The  center  is  designated  R(0,0),  S(0,0),  and  T(0,0). 
Thus,  the  weight  of  each  cell  in  T  can  be  computed  as  follows 


(1)  circle-plus:  T(ti,tj)  =  L  [ R ( ri , r j ) *S ( s i  ,  s j  )  ]  where 

=  { ri  +  si , rj  +  s j  )  and  R(0,0)  is  centered  over 

S ( si , s j ) 

(2)  circle-maximum:  T(ti,tj)  =  max [ R ( r i , r j ) *S ( s i , s j ) ] 
where  (ti,tj)  =  ( ri  +  si , rj  +  s j  )  and  R{0,0)  is  centered  over 
S ( si , s j  ) 

(3)  square-maximum:  T(ti,tj)  =  max [ R ( r i , r j ) +S ( s  l  ,  s j  )  ] 
where  (ti,tj)  =  ( ri +si , r j  +  s j  )  and  R ( 0 , 0 )  is  centered  over 
S  (  s  i  ,  s  j  ) 

A  more  rigorous  mathematical  description  of  these  operations 
can  be  found  in  the  Image  Algebra  Tutorial  (1:30). 

These  six  elemental  template  operations  can  be  used  to 
define  other  operations  between  templates  such  as  subtraction, 
division,  minimum,  circle-minimum,  and  square-minimum  (1:25). 

Algorithm  Optimization 

Although  the  convolution  operations  between  two  templates 
are  complicated  and  appear  useless  in  real  image  processing 
algorithms,  they  form  a  very  important  part  of  the  IA.  These 
convolution  operators  are  important  for  the  composition  and 
decomposition  of  templates  which  provide  tools  for  program 
optimization.  (1:39) 

The  decomposition  of  templates  can  be  used  to  break  a 


large  template  into  an  equivalent  set  of  smaller  templates 
which  yield  a  decrease  in  the  computation  of  an  expression. 

For  example,  the  circle-plus  operation  can  be  used  to  calculate 


the  local  average  gray  level  in  an  image  (dimension  256x256) 
with  a  template  (dimension  3x3)  through  C  =  A  <b  R .  If  template 
R  is  decomposed  into  two  smaller  templates,  S  (dimension  3x1) 
and  T  (dimension  1x3),  then  the  average  image  can  be  calculated 
by  C  =  A  *  R  =  A  ©  ( S  ©  T )  =  ( A  4>  S )  4>  T .  This  equation  shows 
a  number  of  methods  for  calculating  the  average,  each  one  with 
a  very  different  computational  load.  Since  each  pixel  in  C  = 

A  4*  R  requires  the  summation  of  nine  products,  the  entire 
operation  needs  589,824  multiplications  and  524,288  additions. 
However,  since  each  pixel  in  C  =  ( A  <b  S )  4*  T  requires  the 

summation  of  three  products  for  each  circle-plus  operation, 
this  entire  operation  requires  393,216  multiplications  and 
262,144  additions.  Thus,  the  second  formulation  needs  30 
percent  fewer  multiplications  and  50  percent  fewer  additions: 
template  decomposition  can  provide  algorithm  optimization. 


III.  General  Implementation  of  the  Image  Algebra 

The  ultimate  goal  of  the  image  algebra  project  is  a 
standardization  of  the  image  processing  development  tools  used 
by  many  different  companies  and  agencies  around  the  country. 
Because  each  organization  has  a  lot  of  money  invested  in  its 
own  image  processing  tools,  the  only  viable  way  of  obtaining 
standardization  is  by  convincing  each  group  that  the  image 
algebra  is  an  efficient  and  powerful  development  tool. 

The  image  algebra  developed  by  the  University  of  Florida, 
and  briefly  presented  in  Section  II  of  this  report,  provides  a 
strong  foundation  for  an  image  processing  environment  that 
could  become  a  standard.  The  mathematical  structure  and  small 
number  of  elemental  operators  allow  a  user  to  model  any  image, 
perform  any  gray  level  transformation,  and  optimize  algorithms 
through  identities  without  forcing  the  user  through  a  long 
learning  process.  Additionally,  the  notation  provides  an 
understanding  of  the  basic  operations  involved  in  any  image 
processing  algorithm  described  with  the  image  algebra. 

In  essence,  the  University  of  Florida  has  achieved  the 
first  six  goals  of  a  desirable  image  algebra  outlined  by  the 
Air  Force  and  discussed  in  Section  I.  They  have  succeeded  in 


developing  an  image  algebra  that  is  both  powerful  and  concise. 
But,  if  the  image  algebra  is  to  become  a  desirable  development 
tool,  its  implementation  must  retain  this  power  and  brevity  as 


well  as  promote  a  complete  image  processing  environment.  The 
implementation  of  the  image  algebra  is  as  important  to  its 
acceptability  as  its  mathematical  structure  and  power.  Thus, 
what  attributes  must  the  image  algebra  implementation  possess 
in  order  to  become  a  useful  and  desirable  image  processing 
development  tool? 

Image  Algebra  Notation 

The  implementation  of  the  image  algebra  notation  should 
be  concise  and  flexible.  Increases  in  the  production  of  image 
processing  algorithms  can  be  realized  through  a  notation  that 
allows  programmers  to  represent  complicated  operations  with 
single  character  operators.  In  addition,  the  notation  must 
allow  a  programmer  the  flexibility  to  manipulate  the  image 
algebra  operators  as  if  they  were  ordinary  algebraic  operators 
in  a  computer  program. 

The  image  algebra  notation  must  provide  access  to  the 
program  statements  supported  by  many  high  level  languages. 
Program  statements  are  those  structures  that  provide  program 
control  during  execution  of  the  routine.  They  include 
assignment  (A=B),  conditional  branching  ( if-then-else ) ,  and 
repetitive  looping  (for-do,  while-do)  statements.  Because 
each  computer  language  supports  and  implements  these  program 
structures  differently,  the  image  algebra  notation  must 
support  program  statements  suitable  for  translation  to  any 


computer  language. 


The  implementation  should  support  arbitrary  names  for 


image,  template,  and  scalar  operands,  and  program  variables. 
Operand  and  variable  names  should  not  be  limited  to  certain 
lengths  or  certain  letters  of  the  alphabet. 

The  notation  must  support  all  of  the  elemental  operations 
of  the  image  algebra,  and  the  computer  representation  of  each 
operator  should  closely  resemble  its  written  representation. 

A  user  should  be  able  to  enter  an  image  processing  algorithm 
into  a  computer  exactly  as  it  is  written  on  paper. 

Computer  and  Language  Independence 

One  of  the  objectives  of  a  standardized  image  algebra 
processing  environment  is  the  ability  to  transfer  image 
processing  algorithms  from  one  computer  system  to  another 
without  (ideally)  modification.  The  purpose  of  this  goal  is 
to  provide  users  of  the  image  algebra  with  the  ability  to 
share  image  processing  routines  without  forcing  everyone  to 
enter  each  algorithm  by  hand.  The  preferred  method  would  be 
tape  transfers  between  computer  sites.  This  goal  requires 
that  the  image  algebra  environment  be  machine  and  language 
independent . 

There  are  three  basic  approaches  for  making  the  image 
algebra  both  machine  and  language  independent.  First,  require 
every  user  of  the  image  algebra  to  employ  identical  computer 
hardware  and  software.  Second,  implement  the  image  algebra  in 
a  high  level  language.  Third,  design  a  high  level  language, 


based  upon  the  image  algebra  notation,  that  can  be  translated 
to  any  computer  system.  Because  the  ultimate  goal  of  the 
image  algebra  project  is  standardization  of  the  entire  image 
processing  community,  the  third  method  is  the  most  logical 
approach . 

The  first  two  approaches  present  problems  to  wholesale 
acceptance  of  the  image  algebra.  Foremost,  it  is  unlikely 
that  any  organization  would  significantly  modify  its  large 
investment  in  computer  resources  to  accommodate  a  new  program. 
High  level  languages  approach  independence,  but  no  high  level 
language  is  completely  machine  independent.  Further,  a  choice 
of  one  language  or  another  will  present  compatibility  problems 
to  one  group  or  another  because  they  do  not  use  that  specific 
language  . 

The  key  to  general  acceptance  of  the  image  algebra  is  an 
environment  that  can  interface  the  image  algebra  description 
of  an  algorithm  to  the  present  computer  resources  of  each 
organization.  This  interface  would  allow  the  sharing  of  image 
processing  routines  by  translating  image  algebra  descriptions 
into  a  language  supported  by  that  organization’s  computer 
system.  Subsequent  compilation  of  the  translated  routine 
would  allow  its  execution  on  the  host  computer  without 
mod i f icat ion . 

The  real  benefit  of  this  design  is  derived  through  the 
flexibility  of  tailoring  the  image  algebra  operations  to  each 
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organization’s  computer  system  while  retaining  machine  and 


language  independence  through  the  interface.  This  allows  each 
organization  to  optimize  the  IA  operations  while  maintaining 
complete  IA  language  compatibility. 


S 


Image  Algebra  Preprocessor 


The  interface  is  provided  by  a  preprocessor  that 
translates  image  algebra  expressions  into  equivalent  blocks  of 
code  in  FORTRAN,  PASCAL,  ADA,  or  some  other  language.  The 
translated  code  is  subsequently  compiled  and  linked  with 
libraries  of  image  algebra  operations  into  executable  modules. 

The  requirements  for  implementing  the  image  algebra  with 
this  type  of  interface  are  minimal.  Each  organization  needs 
both  a  high  level  language  description  of  the  image  algebra 
operations  and  a  preprocessor.  Each  organization  tailors  the 
description  of  the  image  algebra  operations  to  their  specific 
computer  system  and  then  places  them  into  the  computer’s 
libraries.  Once  the  preprocessor  and  operators  are  installed, 
the  computer  can  be  programmed  to  automatically  translate, 
compile,  and  link  each  routine  into  an  executable  module. 

This  design  requires  the  image  algebra  preprocessor  to 
be  flexible.  The  preprocessor  must  be  capable  of  translating 
image  algebra  expressions  and  other  program  statements  into 
equivalent  code  in  several  different  languages.  This  does  not 
imply  that  every  preprocessor  must  be  capable  of  translating 
image  algebra  notation  into  more  than  one  language.  Each 
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computer  system  will  employ  a  preprocessor  that  translates  the 
image  algebra  expressions  and  program  statements  into  the  high 
level  language(s)  supported  by  that  computer. 

It  is  anticipated  that  image  processing  algorithms  will 
possess  statements  containing  multiple  operations  per  line.  A 
scalar  algebraic  equation  like  c  =  a  *  b  -  (c  +  d)/e  has  an 
equivalent  representation  in  the  image  algebra  as  C  =  A  *  B  - 
(C  +  D)/E.  Rather  than  forcing  the  image  processing  engineer 
to  write  this  equation  as  a  series  of  statements  with  one 
operation  per  line,  it  would  be  desirable  to  enter  the  entire 
equation  on  one  line.  Consequently,  the  preprocessor  should 
be  able  to  translate  multiple  image  algebra  expressions  per 
line. 

There  are  a  number  of  ways  to  implement  the  image  algebra 
preprocessor.  One  method  is  a  preprocessor  that  translates 
each  IA  expression  into  the  equivalent  source  code  of  a  high 
level  language  by  directly  substituting  blocks  of  code  for 
each  IA  expression.  This  method  allows  each  program  module  to 
be  self-contained  (free  from  external  procedures)  but  each 
program  module  may  contain  redundant  blocks  of  code.  This 
method  creates  modules  that  are  likely  to  be  larger  than 
necessary,  but  they  may  execute  more  rapidly  because  they 
avoid  the  inherent  overhead  of  procedure  calls. 


Another  method,  the  one  implemented  in  this  project,  is  a 
preprocessor  that  translates  each  IA  expression  into  a  series 


of  equivalent  high  level  language  procedure  calls.  Although 
each  program  module  must  be  linked  with  the  library  of  image 
algebra  operations,  it  is  a  self-contained  executable  module 
that  prevents  redundant  blocks  of  code.  This  method  creates 
modules  that  are  relatively  compact  but  they  execute  more 
slowly  because  they  rely  almost  entirely  upon  procedure  calls. 
The  major  advantages  of  this  technique  are  the  simple  design 
of  the  preprocessor  and  the  standard  interface  between  the 

The  operation  of  the  preprocessor  implemented  at  AFIT  in 
this  project  is  demonstrated  by  Figures  3.1  and  3.2.  Figure 
3.1  shows  an  IA  source  file  containing  a  number  of  operations. 
The  preprocessor  translates  the  IA  expressions  into  a  series 
of  PASCAL  procedure  calls  as  shown  in  Figure  3.2.  At  this 
point,  it  is  not  necessary  for  the  reader  to  understand  all  of 
the  notation  contained  in  these  two  figures.  That  will  be 
explained  in  the  next  section. 
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Image_and_Template_Operat ions 

(*  This  is  the  program  comment  section.  Each  line  is 
(*  marked  with  the  PASCAL  comment  delimiters  at  the  end 
(*  The  program  comment  section  in  the  IA  file  is  not 
(*  passed  to  the  PASCAL  source  file. 

(*  The  following  lines  contain  the  other  IA  program 
(*  sections.  The  program  code  section  includes  two 
(*  template  definitions. 

(*  constant  declaration  section  *) 
const  MaxLoop  =  3 
(*  type  declaration  section  *) 
type  itype  =  integer 
ttype  =  integer 

(*  variable  declaration  section  *) 
var  R,S,T  :  template 
A, B,C  :  image 
j  :  integer 

( *  program  code  section  * ) 
begin 

(*  define  templates  S  and  T  *) 
invariant  template  S 
begin 

S(x.y)  =4 
S ( x , y+  1  )  =  5 
end 

invariant  template  T 
begin 

T ( x , y-  1  )  =  1 
T  (  x  ,  y  )  =2 
T ( x+ 1 , y )  =  3 
end 

R  :=  S  (+)  T 

Get  Image  ( A ,  ’  ’  ) 

for  j  :=  1  to  MaxLoop  do 

(*  loop  the  statements  delimited  by  "begin"  and  "end" 
beg  i  n 
B  :=  A  >  =  j 

C  : =  (C  +  B)  [v]  R 
end 

Put  Image  ( C  ,  ’  ’  ) 
end . 

Figure  3.1:  Sample  Image  Algebra  File 


[inherits  ' iiaoper.env'  , ’iio.env’  ) ] 

program  image_and_template_operat i ons  ( i nput , output ) ; 
const  maxloop  =  3; 
var  r,s,t  :  template; 
a , b , c  :  image ; 
j  :  integer; 
begin 

reset  ( input ) ; 
rewrite  (output); 

s.fig(.l.).r:=  0;  s.fig(.l.).c:=  0;  s.fig(.l.).w:=4; 

s.fig(.2.).r:=  0;  s.fig(.2.).c:=+l;  s.fig(.2.).w:=5; 

s . num  : =  2 ; 


t.fig(.l.).r:=  0; 

t.fig( .1.  )  . c  :  =-  1 ; 

t . f ig ( . 1 . ) 

.  w  :  =  1  ; 

t . f ig (  . 2 .  )  . r  :  =  0; 

t.fig(.2.).c:=  0; 

t. fig( .2. ) 

.  w  :  =  2  ; 

t.fig(  .3. ) . r : =+ 1 ; 

t.fig(.3.).c:=  0; 

t . f ig (  .3.  ) 

.  w :  =  3  ; 

t . num  : =  3 ; 

( *****  r: =s( + ) t  *****) 

TempCi rc lePlus  (s,t,tRl); 
r  : =  tR 1  ; 
getimage  ( a ,  ’  ’  )  ; 
for  j  : =  1  to  maxloop  do 
begin 

b: =a>= j  ttttt) 
maxval  :=  MaxVallmage  (a); 
Charlmage  ( a , j , maxval , iRl ) ; 
b  :=  iRl; 

{*****  c : = ( c+b ) [ v ] r  **»**) 
ImageAdd  (c,b,iRl); 
IraageTempSquareMax  ( iRl , r , iR2 ) ; 
c  :=  iR2; 
end ; 

putimage  ( c  ,  ’  ’  )  ; 
end . 


Figure  3.2:  PASCAL  Source  Code  After  Preprocessing  the 
Sample  Program  of  Figure  3.1 


AFIT  Implementation  of  the  Image  Algebra 


The  image  algebra,  as  conceived  by  the  University  of 
Florida,  can  be  the  basis  for  a  powerful  image  processing 
development  environment.  It  has  the  capabilities  to  perform 
any  gray  level  image  transformation,  and  its  notational  and 
operational  simplicity  are  surprising.  This  image  algebra  has 
the  potential  to  become  a  widely  accepted  image  processing 
too  1  . 

As  discussed  in  Section  III,  the  most  critical  component 
of  the  image  algebra  will  be  its  computer  implementation, 
especially  the  user  interface.  If  the  image  algebra  is  to 
gain  wide  spread  acceptance,  it  must  have  a  user  interface 
that  is  both  simple  and  capable  of  utilizing  the  full  power  of 
the  image  algebra. 

This  section  discusses  the  design  and  implementation  of 
the  image  algebra  at  AFIT.  The  goal  is  an  image  processing 
development  tool  with  a  simple  user  interface  that  allows  full 
access  to  the  power  of  the  image  algebra.  This  implementation 
is  intended  to  achieve  three  of  the  desirable  properties  of  a 
useful  image  algebra  outlined  in  Section  I:  (1)  substitution 
of  concise  image  algebra  expressions  for  large  blocks  of 
program  code,  (2)  notation  that  allows  the  use  of  libraries  of 
image  transformations,  and  (3)  machine  and  computer  language 
i ndependence  . 
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Development  Environment  and  File  Structure 


Before  discussing  the  explicit  design  of  the  AFIT  image 
algebra,  it  is  helpful  if  the  reader  understands  the  overall 
AFITIA  structure  and  how  it  is  incorporated  into  the  VAX 
11/780  processing  environment  using  VAX  PASCAL. 

The  major  components  of  the  AFITIA  development  tool  are 
derived  from  five  PASCAL  source  code  files:  (1)  two  image 
algebra  operations  files  named  ilAOPER.PAS  and  rIAOPER.PAS, 

(2)  two  input/output  operations  files  named  ilO.PAS  and 
rIO.PAS,  and  (3)  the  preprocessor  file  named  PREPROC . PAS . 

The  first  two  files  contain  the  source  code  for  all  of  the 
operand  type  declarations  and  operations  implemented  in  the 
AFIT  image  algebra  with  integer  and  real  operands.  The  second 
two  files  contain  the  source  code  for  all  of  the  input/output 
operations  with  integer  and  real  operands.  Two  files  each  for 
the  IA  operations  and  input/output  routines  are  required 
because  of  the  strong  data  typing  used  in  PASCAL:  some  IA 

operations  can  not  be  programmed  in  PASCAL  to  manipulate  both 
real  and  integer  operands.  Therefore,  a  PASCAL  implementation 
of  the  IA  requires  a  set  of  IA  operations  each  for  real  and 
integer  operands.  The  "i"  and  "r"  preceding  the  IA  operation 
and  input/output  file  names  above  denote  the  type  of  operands, 
integer  or  real,  used  in  that  file's  procedures.  A  copy  of 
the  source  code  for  the  AFITIA  operations  with  real  operands 
is  in  Appendix  A,  and  a  copy  of  the  source  code  for  the  AFIT 


input./output  operations  with  real  operands  is  in  Appendix  C. 
The  fifth  PASCAL  file  contains  the  AFITIA  preprocessor  used 
for  translating  image  algebra  expressions  into  a  series  of 
PASCAL  procedure  and  function  calls.  A  copy  of  the  source 
code  for  the  AFITIA  preprocessor  is  contained  in  Appendix  B. 

Before  any  translation  of  image  algebra  expressions  can 
be  accomplished  on  the  VAX  11/780,  the  development  environment 
muot  possess  eight  files  derived  from  the  above  source  code 
files.  The  first  two  are  the  object  code  descriptions  of  the 
image  algebra  operations,  ilAOPER.OBJ  and  rIAOPER.OBJ,  created 
by  compiling  ilAOPER.PAS  and  rIAOPER.PAS,  respectively.  The 
next  two  are  the  object  code  descriptions  of  the  input/output 
operations,  ilO.OBJ  and  rIO.OBJ,  created  by  compiling  ilO.PAS 
and  rIO.PAS,  respectively. 

The  other  four  files,  ilAOPER.ENV,  rIAOPER.ENV,  ilO.ENV, 
and  rIO.ENV,  contain  the  procedure  declarations  from  the  image 
algebra  operation  and  input/output  files.  These  files  are 
required  to  ensure  proper  compilation  of  a  translated  IA  file 
because  PASCAL  requires  that  all  procedures  and  functions  be 
declared  to  the  compiler  before  they  are  used.  This  method 
permits  VAX  PASCAL  to  separately  compile  and  link  IA  source 
files,  IA  operation  files,  and  input/output  files.  Like  the 
object  code  files,  these  files  are  created  during  compilation 
of  the  PASCAL  source  code  files  by  using  the  "environment" 
compiler  attribute  in  VAX  PASCAL. 


These  eight  files  plus  an  executable  version  of  the 
preprocessor  are  all  that  is  necessary  to  translate  an  image 
algebra  description  of  any  image  processing  algorithm  into  an 
executable  module.  The  AFITIA  development  tool  uses  a  command 
file  called  TIA.COM  to  automate  this  translation  process.  A 
copy  of  this  command  file  is  in  Appendix  D.  Figure  4.1  shows 
a  block  diagram  of  the  entire  AFITIA  development  tool,  the 
interrelationship  between  all  of  the  files  used  to  build  an 
executable  image  processing  module,  and  the  flow  of  execution 
in  the  automatic  translator  command  file.  The  execution  of 
this  command  file  is  outlined  below. 

Starting  at  the  top,  the  file  containing  the  image 
algebra  language  description,  the  IA  file,  is  copied  to  an 
intermediate  working  file  named  TRANSLAT.IA.  The  name  of  the 
file  may  be  any  valid  filename,  and  it  is  assumed  to  have  a 
” . IA"  extension . 

Next,  the  preprocessor  is  invoked  to  translate  the  IA 
expressions  in  this  working  file  into  PASCAL  source  code. 
First,  the  entire  TRANSLAT.IA  file  is  converted  to  lower  case 
letters,  all  comments  are  removed,  and  the  result  is  placed  in 
TRANSLAT . TMP .  Second,  the  preprocessor  examines  TRANSLAT.TMP 
and  determines  which  type  of  the  image  and  template  operands, 
integer  or  real,  are  used  in  this  routine.  The  default  data 
type  is  real.  Based  upon  this  determination,  the  preprocessor 


writes  a  statement  to  another  working  file,  TRANSLAT . PAS ,  that 


Figure  4.1:  AFITIA  Development  Environment  and  File  Structure 
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allows  the  inclusion  of  the  proper  environment  files,  ilO.ENV 
or  rIO.ENV,  and  ilAOPER.ENV  or  rIAOPER.ENV,  for  compilation. 
Finally,  the  preprocessor  continues  examining  TRANSLAT . TMP , 
translates  each  line  of  the  file  containing  IA  expressions, 
and  outputs  the  translated  code  to  TRANSLAT . PAS . 

Once  the  entire  algorithm  is  translated,  the  PASCAL 
compiler  is  invoked  to  operate  on  the  TRANSLAT. PAS  file,  and 
the  compiled  (object)  code  is  returned  in  TRANSLAT . OBJ .  The 
object  file  is  subsequently  linked  with  the  proper  image 
algebra  operation  and  input/output  object  files  resulting  in 
an  executable  file  named  TRANSLAT.EXE.  This  executable  file 
is  then  renamed  to  the  file  name  of  the  original  IA  file  with 
a  ".EXE"  extension,  and  all  of  the  TRANSLAT.***  files  are 
deleted  from  the  directory. 


Second,  the  AFITIA  operands  are  flexible  enough  to  accommodate 


many  different  image  and  neighborhood  configurations.  Third, 


the  AFITIA  can  perform  any  image  transformation  because  its 


basic  operations  are  identical  to  those  of  the  UFIA.  Recall 


that  the  University  of  Florida  proved  this  set  of  operations 


was  capable  of  performing  any  image- to- image  transformation 


when  used  in  conjunction  with  controlled  program  execution. 


The  following  sections  discuss  each  of  these  aspects  of  the 


AFITIA  language  in  more  detail. 


Program  Control.  The  UFIA  elemental  operations  described 


in  Section  II  can  serve  as  the  basis  for  an  IA  language,  but 


these  operations  alone  cannot  form  a  complete  IA  language 


capable  of  describing  all  transformations.  Many  image 


processing  algorithms  require  additional  program  statements, 


such  as  conditional  branching  and  repetitive  looping,  to 


control  the  execution  of  an  algorithm.  Consequently,  a 


complete  image  algebra  language  must  support  both  the 


operations  defined  in  the  UFIA  and  the  programming  statements 


found  in  the  more  popular  high  level  computer  languages 


Since  the  image  algebra  will  be  used  with  many  different 


computer  languages,  ti.e  image  algebra  language  must  define  and 


support  program  control  statements  that  can  be  translated  into 


any  computer  language.  Designing  a  language  to  be  independent 


of  the  target  language  of  the  preprocessor  will  be  necessary 


because  different  computer  languages  support  different  control 


aft* 


statements.  For  example,  PASCAL  supports  " i f- then-else"  and 
"case"  conditional  statements  and  "for",  "while"  and  "repeat" 
repetitive  statements;  FORTRAN  has  " i f - then-e 1 se  "  ,  "do",  and 
"do-while"  statements;  the  C  language  supports  "  i f- then-else "  , 
"while",  "for",  "do-while",  and  "switch”  statements;  and  PL/I 
has  " if- then-else" ,  "select",  "do",  "while",  and  "until" 
control  statements.  (3:57-61;  7:9-3,9-9,9-12;  8:67-85,151-174; 
9:88-101,103-127)  Although  there  is  a  common  set  of  control 
statements,  there  are  differences  that  make  it  difficult  to 
implement  a  universal  notation. 

Therefore,  this  project  did  not  attempt  to  design  generic 
control  statements  that  could  be  translated  to  equivalent 
statements  in  any  high  level  computer  language.  Instead,  the 
AFITIA  language  retains  the  PASCAL  constructs  in  order  to  form 
an  image  algebra  language  that  supports  repetitive  looping  and 
conditional  branching.  Consequently,  this  implementation  of 
the  image  algebra  language  is  ,iot  completely  machine  or 
language  independent,  but  it  does  progress  in  that  direction. 

Operand  Data  Structures  and  Notation.  The  operands  of 
the  AFITIA  must  be  capable  of  modeling  and  transforming  any 
practical  image.  Consequently,  the  data  structures  of  the 
operands  must  be  flexible  enough  to  handle  many  different 
image  formats  and  neighborhood  configurations. 

In  the  UFIA,  the  number  of  images  and  templates,  and 
their  respective  dimensions,  are  limitless.  Unfortunately,  an 
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implementation  of  the  image  algebra  on  a  computer  requires 
some  compromises.  Implementing  the  image  algebra  on  the  AFIT 
computer  requires  that  the  dimensions  of  the  UFIA  operands  be 
restricted  because  of  limited  computer  resources.  The  size  of 
the  operands  supported  by  the  AFITIA  should  be  sufficient  for 
most  image  processing  tasks. 

In  this  version  of  the  AFITIA,  images  are  limited  to 
two-dimensional  arrays,  and  the  number  of  images  available  to 
the  user  is  software  limited  to  100.  Two-dimensional  arrays 
are  used  to  minimize  the  memory  and  program  requirements. 
Images  may  be  named  with  any  combination  of  characters  (except 
v),  numerals,  and  underscores  up  to  eighty  characters  long.  A 
software  switch,  iType,  included  in  the  IA  source  file  can  be 
used  to  implement  the  images  with  real  or  integer  gray  values. 
The  default  switch  setting  is  real.  Complex  valued  images  are 
not  supported  in  this  version  of  the  AFITIA.  Each  dimension 
of  the  image  may  be  any  size  up  to  a  maximum  value  set  by  the 
user  with  two  other  software  switches  found  in  the  IAOPER.PAS 
files,  MaxImageRow  and  Max ImageCo 1 umn . 

The  data  structure  of  the  image  operand  is  a  record  with 
three  fields:  row  dimension,  column  dimension,  and  an  array  of 
integer  gray  levels  indexed  by  row  and  column  position.  This 
data  structure  allows  the  AFITIA  to  manipulate  multiple, 
arbitrarily  sized  images.  Additionally,  since  the  dimensions 
of  the  image  can  be  initialized  by  a  procedure  that  reads  in 
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the  input  image,  the  user  is  freed  from  tracking  ail  of  the 
image  dimensions  during  execution. 

The  image  notation  used  in  this  project  and  by  the 
University  of  Florida  is  similar.  Images  are  denoted  with 
capital  letters  from  the  beginning  of  the  alphabet  such  as  A, 
B,  and  C.  Arbitrary  pixels  in  an  image  are  represented  with 
the  vector  notation  k,  and  the  gray  value  at  an  arbitrary 
pixel  in  image  A  is  denoted  by  A(k).  Specific  pixels  in  the 
two-dimensional  AFITIA  images  are  denoted  with  row  and  column 
qualifiers  to  the  image  letter  designator  such  as  A ( 2 , 3 )  or 
B  (  7  ,  4  )  . 

In  this  version  of  the  AFITIA,  one-dimensional  arrays  are 
used  for  storing  template  configurations  and  weights.  One¬ 
dimensional  arrays  were  chosen  for  the  template  operands  in 
order  to  simultaneously  increase  configuration  flexibility  and 
minimize  computer  storage  requirements.  Li^e  image  operands, 
templates  may  be  named  with  any  combination  of  letters  (except 
v),  numerals,  and  underscores  up  to  eighty  characters  in 
length.  A  software  switch,  tType ,  included  in  the  IA  source 
file  can  be  used  to  implement  the  templates  with  real  or 
integer  weights.  The  default  switch  setting  is  real.  Complex 
valued  templates  are  not  supported  in  this  version  of  the 
AFITIA.  The  size  of  the  templates  is  presently  limited  by 
another  software  switch,  MaxTempCell,  in  the  IAOPER.PAS  file 


to  100  pixels  in  each  configuration.  All  pixels  not  within 
the  configuration  are  defined  to  have  a  value  of  zero. 

The  data  structure  of  the  template  operand  is  a  record 
with  two  fields.  The  first  field  contains  the  number  of 
pixels  in  the  configuration.  The  second  field  is  an  array 
with  each  array  cell  representing  a  pixel  in  the  template 
configuration.  Each  array  cell  contains  a  record  with  three 
fields  of  information  about  that  pixel:  the  (1)  row  and  (2) 

column  offset  from  that  template’s  center  pixel,  and  (3)  the 
weight  of  that  pixel.  This  data  structure  permits  arbitrary 
row  and  column  offsets  for  each  pixel  and,  therefore,  allows 
arbitrary  template  configurations  to  be  implemented  while 
minimizing  the  computer  storage  requirements. 

The  template  notation  used  in  this  project  and  by  the 
University  of  Florida  is  similar.  Templates  are  denoted  with 
capital  letters  from  the  end  of  the  alphabet  such  as  R,  S,  and 
T.  Arbitrary  pixels  in  a  template  are  represented  with  the 
vector  notation  z  where  z  stands  for  the  m-dimens ional  offset 
from  the  template’s  center  pixel.  (Recall  that  the  center 
pixel  is  not  necessarily  the  physical  center  of  the  template’s 
configuration.)  The  weight  at  an  arbitrary  pixel  in  template 
R  is  denoted  by  R(z).  Specific  pixels  in  AFITIA  template 
configurations  are  represented  by  row  and  column  offsets  from 
the  template’s  center  pixel.  For  example,  R(-2,3)  is  a  pixel 
two  rows  above  and  three  columns  to  the  right  of  the  center 
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pixel,  and  the  center  pixel  of  template  S  may  be  represented 
by  S ( 0 , 0 )  or  S ( 0 ) . 

Operator  Implementation  and  Notation.  The  five  elemental 
image-image  operations,  three  elemental  image- template 
operations,  and  six  elemental  template- template  operations  of 
the  UFIA  are  preserved  in  the  AFITIA.  Each  operation  is 
implemented  as  a  PASCAL  procedure  or  function.  Each  basic 
operation  is  performed  by  calling  the  appropriate  subroutine, 
passing  the  necessary  operands,  and  obtaining  the  result 
returned  from  the  subroutine.  The  source  code  for  all  of  the 
AFITIA  operations  is  in  Appendix  A. 

This  implementation  of  the  operators  directly  supports 
the  concise  notation  feature  described  above.  The  lar^e 
blocks  of  code  needed  to  program  each  operation  are  replaced 
by  a  short  subroutine  call. 

The  interface  between  the  program  calling  the  procedure 
and  the  procedure  itself  is  designed  to  provide  maximum 
program  flexibility  while  ensuring  data  and  variable 
integrity.  Using  the  parameter  passing  schemes  provided  by 
PASCAL,  each  procedure  obtains  its  operands  through  value 
parameters  and  the  results  are  returned  to  the  calling  program 
through  variable  parameters.  Local  variables  are  used  in  each 
procedure  to  prevent  the  procedures  from  modifying  any 
variables  not  within  their  scope.  This  scheme  allows  each 
procedure  to  perform  its  operation  without  modifying  the  data 
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in  either  operand,  and  it  allows  the  result  to  be  returned  to 
an  arbitrary  variable  name  in  the  calling  program. 

As  noted  above,  the  AFITIA  is  comprised  of  a  collection 
of  PASCAL  procedures  and  functions  that  allow  the  substitution 
of  short  subroutine  calls  for  the  computer  code  of  each 
elemental  operation.  This  attribute  partially  supports  the 
sixth  image  algebra  guideline  outlined  in  Section  I.  However, 
requiring  a  user  to  write  image  processing  algorithms  as  a 
collection  of  procedure  calls  is  not  compatible  with  the  goal 
of  writing  algorithms  using  the  image  algebra  notation. 

Writing  computer  programs  in  the  image  algebra  notation 
requires  a  mapping  of  image  algebra  operator  symbols  to  ASCII 
characters.  Because  many  of  the  operator  symbols  in  the  image 
algebra  cannot  be  represented  by  a  single  ASCII  character,  the 
AFITIA  uses  short  strings  of  characters  that  closely  resemble 
the  image  algebra  notation.  Figure  4.2  displays  all  of  the 
operators,  their  UFIA  notation,  and  their  equivalent  AFITIA 
notation.  In  the  figure,  A,  B,  and  C  are  defined  as  images, 

R,  S,  and  T  are  defined  as  templates,  and  f  is  defined  as  a 
scalar.  Notice  the  desirable  similarity  between  the  operator 
symbols  of  the  UFIA  and  AFITIA.  This  scheme  avoids  memorizing 
a  mapping  between  image  algebra  operator  symbols  and  ASCII 
characters.  Additionally,  this  should  ease  the  learning 
necessary  to  write  computer  programs  using  the  image  algebra. 


operation 

Image-Image  Operations 
UFIA  notation 

AFITIA  notation 

addition 

C  =  A  +  B 

C  =  A  +  B 

subtraction 

C  =  A  -  B 

C  =  A  -  B 

mul t ipl ica t ion 

C  =  A  *  B 

C  =  A  *  B 

division 

C  =  A  /  B 

C  =  A  /  B 

maximum 

C  =  A  v  B 

C  =  A  v  B 

minimum 

C  =  A  a  B 

C  =  A  ~  B 

exponential 

C  =  A  **B 

C  =  A  *  *B 

dot  product 

f  =  A  *  B 

f  =  A  .  B 

Image-Template  Operations 

operat ion 

UFIA  notation 

AFITIA  notation 

circle-plus 

C  =  A  <b  R 

C  =  A  (  +  )  R 

c i rc 1 e-max i mum 

C  =  A  v  R 

C  =  A  (v)  R 

c i rc le-min imura 

C  =  A  <&  R 

C  =  A  (  *  )  R 

s qua re -max i mum 

C  =  A  53  R 

C  =  A  [v]  R 

square-minimum 

C  =  A  H  R 

C  =  A  [  ~  ]  R 

Template-Template  Operations 

operation 

UFIA  notation 

AFITIA  notation 

addition 

T  =  R  +  S 

T  =  R  +  S 

subtraction 

T  =  R  -  S 

T  =  R  -  S 

multiplication 

T  =  R  *  S 

T  =  R  *  S 

division 

T  =  R  /  S 

T  =  R  /  S 

maximum 

T  =  R  v  S 

T  =  R  v  S 

minimum 

T  =  R  a  S 

T  =  R  ~  S 

circle-plus 

T  =  R  e  S 

T  =  R  (+)  S 

circle-maximum 

T  =  R  &  S 

T  =  R  (v)  S 

circle-minimum 

T  =  R  ®  S 

T  =  R  (  "  )  S 

square -maxi mum 

T  :  R  B  S 

T  =  R  [v]  S 

square-minimum 

T  =  R  a  s 

T  =  R  [“]  S 

Other  Operations 

operation 

UFIA  notation 

AFITIA  notation 

ass ignment 

A  =  B 

A  :  =  B 

R  =  S 

R  :=  S 

characteristic 

A  =  B 

A  =  B 

images 

A  >  B 

A  >  B 

( image- image ) 

A  >  B 

A  >=  B 

A  <  B 

A  <  B 

A  <  B 

A  <=  B 

A  *  B 

A  <  >  B 

characteristic 

A  =  f 

A  =  f 

images 

A  >  f 

A  >  f 

( image-scalar ) 

A  >  f 

A  >=  f 

A  <  f 

A  <  f 

A  <  f 

A  <=  f 

A  *  f 

A  <  >  f 

absolute  value 

:  a: 

:  a: 

:  r  : 

:r: 

Figure  4.2: 


AFIT  Image  Algebra  Operator  Syntax 


AFIT  Image  Algebra  Preprocessor 


If  the  image  algebra  notation  is  to  be  used  for  writing 
computer  programs,  then  an  interface  between  the  image  algebra 
language  and  the  host  computer  language  must  be  provided. 

This  interface  must  be  capable  of  translating  image  algebra 
operands  and  operator  symbols  into  the  more  popular  high  level 
computer  languages. 

Consequently,  the  final  component  of  the  AFIT  image 
algebra  environment  is  a  preprocessor  that  translates  image 
algebra  expressions  into  to  a  sequence  of  PASCAL  procedure 
calls  required  to  execute  the  algorithm.  This  preprocessor 
allows  a  user  to  write  an  image  processing  algorithm  entirely 
with  image  algebra  operands,  operators,  and  common  programming 
constructs,  and  with  little  knowledge  of  computer  programming 
syntax.  The  source  code  for  the  AFITIA  preprocessor  is  in 
Appendix  B. 

The  AFITIA  preprocessor  is  designed  to  translate  any 
unambiguous  image  algebra  expression,  using  the  notation  shown 
in  Figure  4.2,  into  an  equivalent  sequence  of  PASCAL  procedure 
calls.  This  preprocessor  also  supports  any  of  the  PASCAL 
programming  constructs  such  as  assignment,  branching,  and 
repetitive  looping.  The  development  of  an  image  processing 
algorithm  using  the  AFIT  image  algebra  notation  with  the  AFIT 


preprocessor  requires  the  user  to  follow  a  few  guidelines  when 


writing  the  program. 
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AFIT  Image  Algebra  Syntax 

The  AFITIA  preprocessor  expects  the  image  algebra  file  to 
be  comprised  of  up  to  six  basic  sections:  (1)  program  name, 

(2)  program  comment,  (3)  constant  declaration,  (4)  type 
declaration,  (5)  variable  declaration,  and  (6)  program  code. 
Each  section  is  delimited  by  a  keyword  or  symbol,  but  all 
sections  need  not  be  present  to  form  a  working  algorithm. 

Only  the  sections  necessary  for  proper  compilation  are 
required.  Additionally,  the  order  of  the  sections  in  the 
image  algebra  file  is  very  flexible.  The  only  requirements 
are  that  the  program  name,  if  it  is  included,  must  be  first 
section  and  the  program  code  must  be  last  section.  Any 
comments  embedded  in  the  IA  file  are  ignored.  Figure  4.3  shows 
an  IA  program  with  all  six  sections  included. 

Program  Name.  An  image  algebra  program  may  be  named  in 
the  IA  source  file,  but  it  is  not  required.  The  preprocessor 
interprets  the  first  non-blank  line  of  the  file  as  the  program 
name  unless  the  first  line  contains  a  comment  delimiter  or  a 
word  reserved  to  head  any  of  the  major  file  sections.  In  the 
latter  case,  the  program  is  assumed  to  be  unnamed. 

Program  Comment.  The  comment  section  is  an  optional 
section  for  all  of  the  comments  the  user  desires  to  imbed  in 
the  IA  file  about  the  program’s  operands,  execution,  input, 
and  output.  Each  comment  line  is  delimited  by  "(*"  on  the 
left  end  and  "*)"  on  the  right  end.  The  length  and  content  of 
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Image_and_Teraplate_Operat ions 

(*  This  is  the  program  comment  section.  Each  line  is 
(*  marked  with  the  PASCAL  comment  delimiters  at  the  ends. 
(*  The  program  comment  section  in  the  IA  file  is  not 
(*  passed  to  the  PASCAL  source  file. 

{*  The  following  lines  contain  the  other  IA  program 
(*  sections.  The  program  code  section  includes  two 
(*  template  definitions. 

(*  constant  definition  section  *) 
const  MaxLoop  =  3 

(*  type  definition  section  *) 
type  itype  =  integer 
ttype  =  integer 

(*  variable  declaration  section  *) 
var  R,S,T  :  template 
A,B,C  :  image 
j  :  integer 


( *  program 
begin 
( *  def ine 
invariant 
begin 
S(x,y) 

S { x , y+ 1 ) 
end 

invariant 
begin 
T ( x , y-  1  ) 
T  (  x  ,  y  ) 

T ( x+ 1 , y ) 

end 


code  section  *  ) 

templates  S  and  T 
template  S 

=  4 
=  5 


template  T 

=  1 
=  2 
=  3 


*> 


R  :  =  S  (  +  )  T 
Getlmage  ( A , ’ ’ ) 
for  j  : =  1  to  MaxLoop  do 
(*  loop  the  statements  delimited 
beg  i  n 
B  :=  A  >  =  j 

C  :  =  (C  +  B)  ( v]  R 
end 

Putlraage  ( C  ,  ’  ’  ) 
end . 


by 


"beg i n " 


and 


"end"  *) 


Figure  4.3:  Sample  Image  Algebra  File 


the  program  comment  section  is  arbitrary  because  the  AFIT 
preprocessor  ignores  ail  of  the  comments  in  the  IA  file.  This 
implies  that  comments  may  be  placed  anywhere  in  the  IA  file. 
The  comment  section  is  terminated  when  the  preprocessor  finds 
any  of  the  reserved  words  used  to  head  another  section  of  the 
I A  f ile . 

Constant  Declaration.  The  constant  declaration  section 
is  used  to  declare  and  equate  identifiers  with  constant 
values.  This  optional  section  is  reserved  for  initializing 
user-defined  constants  to  be  used  in  the  program.  The  entire 
section  is  headed  by  the  reserved  word  "const",  followed  by  a 
list  of  constants  assigned  to  identifiers.  Each  constant  is 
assigned  to  an  identifier  or  value  with  an  equal  sign.  This 
is  shown  in  the  constant  declaration  section  of  Figure  4.3. 
Note  the  reserved  word  "const"  is  required  only  in  the  first 
line  of  the  section.  The  constant  declaration  section  is 
concluded  when  the  preprocessor  encounters  a  reserved  word 
used  to  head  one  of  the  other  IA  file  sections. 

Type  Declaration.  The  type  declaration  section  is  an 
optional  part  of  the  image  algebra  file  which  is  used  to 
declare  user-defined  data  types.  The  syntax  of  the  type 
declaration  section  is  identical  to  that  of  the  constant 
declaration  section  except  the  reserved  word  "type"  is  used  in 
place  of  "const".  See  Figure  4.3  for  an  example.  The 


software  switches,  iType  and  tType ,  are  placed  here  to  define 


the  data  type  of  the  image  gray  levels  and  template  weights, 
respectively-  The  default  type  is  real.  The  type  declaration 
section  is  terminated  when  the  preprocessor  finds  a  reserved 
word  used  to  head  one  of  the  other  IA  file  sections. 

Variable  Declaration.  The  variable  declaration  section 
is  a  required  part  of  the  image  algebra  file.  All  variable 
operands  used  in  the  program,  such  as  images,  templates,  and 
scalars,  must  be  declared  (identified)  to  the  compiler  in  this 
section.  The  variable  declaration  section  is  headed  by  the 
reserved  word  "var"  on  the  first  line  followed  by  one  or  more 
lines  listing  one  or  more  variables  separated  by  commas.  Each 
line  of  variables  is  terminated  with  a  colon  and  a  data  type. 
The  data  types  may  be  either  user  defined  types  from  the  type 
declaration  section,  types  defined  by  the  image  algebra  such 
as  image  or  template,  or  they  may  be  standard  PASCAL  data 
types  like  real,  integer,  or  character.  Figure  4.3  shows  some 
examples  of  variable  declarations  for  both  IA  defined  and 
standard  PASCAL  data  types.  The  variable  declaration  section 
is  concluded  when  the  preprocessor  finds  a  reserved  word  used 
to  head  one  of  the  other  IA  file  sections. 

Program  Code.  The  program  code  section  is  the  last  part 
of  the  image  algebra  file.  This  is  where  the  image  algebra 
operands  are  manipulated  by  the  image  algebra  operators  and 
the  flow  of  program  execution  is  controlled.  The  program  code 
section  is  headed  by  the  reserved  word  "begin"  and  terminated 
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with  the  reserved  word  "end."  In  between  these  two  words  can 


be  as  many  image  algebra  and  program  control  statements  as 
necessary  to  accomplish  the  image  processing  task. 

The  syntax  of  the  image  algebra  expressions  is  relatively 
straightforward.  Any  expression  may  be  entered  into  the  IA 
file  as  it  is  written  on  paper  using  the  AFITIA  symbols  shown 
in  Figure  4.2.  Few  restrictions  are  placed  on  program  lines 
containing  image  algebra  expressions.  Each  line  may  contain 
multiple  image  algebra  operations.  The  preprocessor  is  set  up 
to  interpret  each  line  with  a  lef t-to-right  hierarchy  within 
parentheses;  that  is,  the  order  of  operation  is  determined 
from  left  to  right  within  any  level  of  parentheses.  One 
restriction  is  that  all  parentheses  be  written  using  the 
and  "}"  symbols.  Examples  of  image  algebra  expressions  with 
multiple  operations  per  line  can  be  found  in  the  sample  image 
algebra  source  file  shown  in  Figure  4.3. 

The  initialization  of  template  operands  is  a  simple  and 
flexible  process.  A  template  definition  is  headed  by  the 
reserved  words  "invariant  template"  followed  on  the  same  line 
by  the  name  of  the  template.  The  next  line  contains  the 
reserved  word  "begin"  to  denote  the  start  of  the  definition. 
The  configuration  and  weights  are  initialized  on  successive 
lines  as  pairs  of  row  and  column  offsets  with  a  weight 


assigned  to  each  offset  pair.  The  syntax  of  each  line  is 
" ( r  +  x , c  +  y )  =z"  where  x  and  y  are  integer  values  and  z  is  an 
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integer  or  real  value  depending  upon  the  data  type  of  the 
template  weights.  Template  definitions  are  terminated  with 
the  reserved  word  "end".  Figure  4.3  shows  the  definition  of 
two  templates,  S  and  T. 

Some  additional  syntax  guidelines  are  carried  over  from 
the  PASCAL  language.  These  rules  affect  program  execution 
because  of  the  way  the  PASCAL  compiler  interprets  program 
control  statements.  Basically,  PASCAL  requires  that  each 
group  of  logically  related  expressions  be  delimited  with  a 
"begin"  and  "end"  statement.  Logically  related  expressions 
are  those  lines  of  program  code  under  the  control  of  the  same 
program  statement  such  as  a  loop  or  conditional  branch.  An 
example  of  this  syntax  requirement  can  be  found  in  the  "for" 
loop  at  the  end  of  Figure  4.3. 

Program  control  statements  may  not  contain  image  algebra 
operators  but  they  may  contain  explicit  IA  procedure  calls. 
The  preprocessor  is  designed  not  to  interpret  any  program 
lines  containing  reserved  words  or  standard  identifiers  from 
the  PASCAL  language.  Therefore,  any  image  algebra  operators 
on  the  same  line  as  a  reserved  PASCAL  word  or  PASCAL  standard 
identifier  are  passed  directly  to  the  TRANSLAT.PAS  file  and 
may  result  in  errors  when  they  are  compiled.  However,  any 
AFITIA  procedure  calls  in  a  program  control  statement  can  be 
compiled  successfully. 
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A  list  of  reserved  words  and  standard  identifiers  from 
both  the  image  algebra  and  PASCAL  languages  is  shown  in  Figure 
4.4.  The  reserved  words  are  primitive  components  of  both  the 
image  algebra  and  PASCAL  which  cannot  be  used  as  identifiers 
in  an  IA  program.  On  the  other  hand,  the  PASCAL  standard 
identifiers  listed  in  Figure  4.4  may  be  redefined  in  an  IA 
program,  but  the  predefined  facility  of  that  identifier  is 
then  lost.  (3:37-38)  Image  algebra  standard  identifiers 
should  be  treated  as  reserved  words . 

Image  and  template  input/output  operations  are  handled 
through  PASCAL  procedure  calls.  The  input/output  routines 
supported  by  the  AFITIA  are  shown  in  Figure  4.5.  Both  read 
operations  prompt  the  user  for  the  name  of  the  file  from  which 
to  read.  Both  write  operations  write  the  image  to  a  file 
named  by  the  user.  If  the  name  of  the  file  in  the  procedure 
call  is  left  blank,  the  user  is  prompted  for  the  file  name 

during  execution.  Appendix  C  contains  the  source  code  for  the 
AFITIA  input/output  procedures. 

The  AFITIA  contains  some  extensions  to  the  UFIA  set  of 
operators  to  ease  the  programming  of  algorithms  in  the  AFIT 
image  algebra  language.  In  addition  to  the  fourteen  elemental 
operators  outlined  in  Section  II,  the  AFITIA  includes  three 
functions  and  one  procedure  for  image  operands,  and  four 
procedures  for  template  operands.  These  operations  are  shown 


VAX  PASCAL  Reserved  Words 


and 

array 

begin 

case 

const 

di  v 

do 

aownto 

else 

end 

externa 1 

file 

for 

forward 

f  unct i on 

goto 

if 

i  n 

label 

mod 

nil 

not 

of 

or 

packed 

procedure 

program 

record 

repeat 

set 

then 

to 

type 

until 

var 

varying 

wh i  le 

with 

VAX  PASCAL  Standard 

Identi f iers 

arctan 

boolean 

byte 

char 

chr 

close 

cos 

eof 

eo  In 

exp 

false 

index 

input 

integer 

length 

In 

odd 

ord 

output 

pred 

read 

readln 

readvar 

reset 

rewri te 

round 

sin 

sqr 

sqr  t 

succ 

text 

true 

trunc 

wr  i  te 

writeln 

writeva 

Image  Algebra  Reserved  Words  and  Standard  Identifiers 


ConfigTempHConst 
Conf i gTempVNConst 
Getlmage 
ImageEqual 
i  R 1 
i  R4 

Max ImageColumn 

maxval 

M i nVa 1 Image 

template 

tR2 

tR5 


Conf igTempVConst 

Constlmage 

GetTemplate 

ImageType 

iR2 

iR5 

MaxImageRow 

MaxVal Image 

Putlmage 

TemplateType 

tR3 

tType 


Conf igTempMooreConst 

do  tval 

image 

invariant 

iR3 

iType 

MaxTempCe 1 1 

minval 

PutTemp] ate 

tRl 

tR4 

v 


Figure  4.4:  Reserved  Words  and  Standard  Identifiers 
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Input  /  Output  Operati ons 


operation 


AFITIA  notation 


read  image  A  from  the 
external  file  ’name* 


Getlmage  (A, 'name' ) 


write  image  A  to  the 
external  file  ’name* 


Putlmage  (A, 'name') 


read  template  R  from  the 
external  file  'name' 


GetTemplate  ( R , ’ name ’ ) 


write  template  R  to  the 
external  file  ’name’ 


PutTemplate  (R, ’name’ ) 


Image  Operations 


operation 


AFITIA  notation 


is  image  A  =  image  B  ? 


minimum  gray  value  of  image  A 

maximum  gray  value  of  image  A 

initialize  image  A  to  dimension 
(row  x  col)  with  constant 
pixel  gray  value  of  X 


ImageEqual  (A,B) 
MinVallmage  (A) 

MaxVal Image  (A) 
Constlmage  ( A , row , col , x ) 


Template  Operations 


operation 

configure  R  as  a  horizontal 
template  of  dimension  (1  x  col) 
with  constant  pixel  weights  of  x 

configure  R  as  a  vertical 
template  of  dimension  (row  x  1) 
with  constant  pixel  weights  of  x 

configure  R  as  a  Moore  template 
of  dimension  (row  x  col)  with 
constant  pixel  weights  of  x 

configure  R  as  a  Von  Neumann 
template  of  radius  (rad)  with 
constant  pixel  weights  of  x 


AFITIA  notation 


Conf igTempHConst  (R,col,x) 


Conf igTempVConst  (R,row,x) 


Conf igTempMooreConst  (R, 
row , col , x ) 


Conf i gTerapVNConst ( R , rad , x ) 


Building  an  Executable  Module 

Once  the  image  processing  routine  is  properly  described 
in  the  image  algebra  language,  it  can  be  transformed  into  an 
executable  module  with  ease  by  invoking  the  VAX  command  file 
TIA.COM  followed  by  the  name  of  IA  file  to  be  transformed. 

This  command  file  automatically  executes  the  preprocessor  to 
translate  the  IA  file  into  a  PASCAL  source  code  file,  executes 
the  PASCAL  compiler  to  transform  the  source  code  into  an 
object  file,  links  the  object  file  with  the  image  algebra  and 
I/O  routines  contained  in  IA0PER.0BJ  and  10. OBJ  to  produce  an 
executable  module,  and,  finally,  renames  the  executable  module 
with  the  name  of  the  original  IA  file.  TIA.COM  is  listed  in 
the  appendix. 

Figure  4.6  shows  the  result  of  preprocessing  the  sample 
program  from  Figure  4.3.  After  compiling  and  linking  this 
source  code,  it  becomes  a  self-contained  executable  module. 
Execution  of  the  image  processing  task  is  accomplished  by 
running  the  executable  module  with  the  VAX  command  "RUN 


f i lename" . 


[inherit* ' iiaoper.env’ , ' iio.env' ) ] 

program  image_and_tempiate_operat ions  ( input , output ) ; 
const  maxloop  =  3; 
var  r,s,t  :  template; 
a,b,c  :  image; 
j  :  integer; 
begin 

reset  ( input ) ; 
rewrite  (output); 

s.fig(.l.).r:=  0;  s.fig(.l.).c:=  0;  s.fig(.l.).w:=4; 

s.fig(.2.).r:=  0;  s.fig(.2.).c:=+l;  s.fig(.2.).w:=5; 

s .  num  : =  2 ; 

t. fig(.l.).r:=  0;  t.fig(.l.).c:=-l;  t.fig(.l.).w:=l; 

t.fig(.2.).r:=  0;  t.fig(.2.).c:=  0;  t.fig(.2.).w:=2; 

t.fig(.3.).r:=+l;  t.fig(.3.).c:=  0;  t.fig(.3.).w:=3; 

t . num  : =  3  ; 

( ttttt  r:=s(+)t  ttttt) 

TempC i rc leP lus  (s,t,tRl); 
r  : =  tR 1 ; 
getimage  ( a , ’ ’ ) ; 
for  j  : =  1  to  maxloop  do 
begin 

(****«  b: =a>= j 

maxval  :=  MaxVallmage  (a); 

Char  Image  ( a , j , maxval , iRl )  ; 
b  :=  iRl; 

(ttttt  c : = ( c+b ) [ v ] r  ****%) 

ImageAdd  (c,b,iRl); 

ImageTempSquareMax  (iRl,r,iR2); 
c  :=  iR2 ; 
end ; 

put  .mage  (c,"); 
end . 


Figure  4.6:  PASCAL  Source  Code  After  Preprocessing  the 
Sample  Program  of  Figure  4.3 


To  demonstrate  the  power  of  the  image  algebra  language 
and  the  simplicity  of  the  AFITIA  development  tool,  this 
section  develops  three  image  processing  algorithms.  The  first 
algorithm,  a  mean  filter,  demonstrates  the  image  algebra’s 
power  for  describing  linear  image  processing  tasks.  The  other 
two  algorithms,  a  median  filter  and  a  local  mode  filter, 
exhibit  the  image  algebra’s  power  for  representing  nonlinear 
image  processing  tasks. 

A  comparison  of  the  two  implementations,  one  in  the  image 
algebra  and  one  in  PASCAL,  of  each  of  the  following  algorithms 
demonstrates  the  tremendous  potential  productivity  increases 
by  a  programmer  using  the  image  algebra.  Due  to  the  concise 
notation  and  powerful  operators,  an  image  processing  engineer 
can  produce  reliable  image  transformation  tools  with  minimal 
development  time.  Further,  these  comparisons  highlight  the 
ease  of  writing  image  processing  algorithms  with  the  image 
algebra  language.  They  also  indicate  some  areas  where  the 
AFITIA  needs  improvement. 

Mean  Filter 

The  mean  filter  is  a  linear  operation  used  for  noise 
suppression  in  images.  The  filter  modifies  each  image  pixel 
to  reflect  the  mean  gray  value  within  a  small  neighborhood  of 
each  image  pixel  where  the  neighborhood  is  defined  by  the 


program  Mean_Filter  ( input , output ) ; 

(*  Let  image  A  be  a  two  dimensional,  integer  valued  array  *) 
(*  of  dimensions  Image_Rows  by  Image_Columns .  *) 

(*  Let  filter  F  be  a  3x3  array  of  constant  weight  =  1.  *) 

(*  The  result  of  the  filter  operation  is  placed  in  image  R  *) 
type  Image  =  array  (  . 1 .  . 256 , 1 . . 256  .  )  of  integer; 

Template  =  array  (  .  1  .  .  3  ,  1  .  .  3  .  )  of  integer; 
var  A,R  :  Image; 

F  :  Template; 

ImageRows , Image_Columns , i , j , k , 1 , sum  :  integer; 

^Include  ’IO_PAS.PAS' 
begin  (*  program  Mean_Filter  *) 
reset  ( input )  ; 
rewrite  (output); 

(*  read  image  A  from  an  external  file  *) 

Get  Image  (A,  ’ input . img ’  ,Image_Rows, Image_Columns )  ; 

(*  initialize  the  filter  weights  to  1  *) 
fori  :=  1  to  3  do 
for  j  :=  1  to  3  do  F(.i,j.)  :  =  1; 

(*  for  each  pixel  in  image  A  *) 
for  i  : =  1  to  Image_Rows  do 
for  j  : =  1  to  Image_Columns  do  begin 
( *  sum  the  gray  values  in  the  neighborhood  * ) 
sum  : =  0 ; 

for  k  :=  -1  to  1  do 
for  1  :=  -1  to  1  do 

if  (i+k>=l)  and  ( i +k < = Image_Rows )  and 
(j+l>=l)  and  ( j  +  1 < = Image_Columns ) 
then  sum  :=  sum  +  A(  . i+k , j  +  1 . )  *F( .k  +  2, 1+2.  )  ; 

(*  divide  the  sum  by  the  number  of  neighborhood  pixels  *) 
(*  and  place  this  value  into  the  resulting  image  *) 

R (  .  i  ,  j  .  )  :  =  round  ( sum/ 9 )  ; 

end ;  ( *  for  j  * ) 

(*  write  the  resulting  image  to  OUTPUT. IMG  *) 

Put  Image  (R, ’output . irag '  , I mage  Rows , Image_Co lumns )  ; 
end.  (*  program  Mean_Filter  *) 

Figure  5.1:  PASCAL  Implementation  of  a  Mean  Filter 


configuration  of  the  filter.  This  filter  has  the  desirable 
effect,  of  reducing  image  noise,  but  it  introduces  blurring  at. 
step  or  ramp  edges  (6:330-331). 

Figure  5.1  shows  a  PASCAL  implementation  of  the  mean 
filter,  including  the  necessary  input/output  operations  to 


make  it  useful  as  an  executable  module.  The  algorithm  for 
this  filter  is  simple.  Following  the  declaration  sections, 
the  routine  reads  an  image  from  a  user-specified  external 
file.  The  external  file  is  assumed  to  be  headed  by  the  row 
and  column  dimensions  of  the  image  followed  by  the  image  data 
arranged  in  a  two  dimensional  (row  x  column)  array.  After 
obtaining  the  input  image,  the  program  steps  through  each 
pixel  of  the  image.  At  each  pixel  of  the  image,  designated  by 
A ( i , j ) ,  all  of  the  pixel  gray  values  lying  under  the  filter 
configuration  centered  over  A(i,j)  are  summed  and  divided  by 
the  number  of  pixels  in  the  filter  configuration.  This 
quotient  replaces  the  gray  value  at  pixel  A(i,j). 

Figure  5.2  shows  the  mean  filter  algorithm  described  in 
the  AFITIA  language.  The  operation  of  this  algorithm  is 
subtle  compared  to  the  direct  PASCAL  implementation.  The 
declaration  and  initialization  sections  are  similar.  After 
reading  the  image  from  an  external  file,  this  algorithm 
replaces  each  pixel  gray  value  with  the  sum  of  pixel  gray 
values  in  the  neighborhood  of  each  image  pixel.  Once  this 
operation  is  completed  by  the  circle-plus  operator,  a  simple 
scalar  division  is  needed  to  complete  the  averaging  process: 
divide  the  image  by  the  number  of  pixels  in  the  neighborhood. 

A  comparison  of  the  Figures  5.1  and  5.2  demonstrates  the 
advantages  of  the  image  algebra.  As  shown,  the  image  algebra 
routine  contains  5  lines  of  executable  (non-comment)  code 


Mean_F i 1 te  r 

(*  Let  image  A  be  a  two  dimensional,  integer  valued  array  *) 
(*  of  dimensions  Image_Rows  by  Iraage_Columns .  *) 

(*  Let  filter  F  be  a  3x3  array  of  constant  weight  =1.  *) 

(*  The  result  of  the  filter  operation  is  placed  in  image  A  *) 
type  itype  =  integer 
ttype  =  integer 
var  A  :  image 

F  :  template 

beg  i  n 

(*  read  image  A  from  an  external  file  *) 

Getlmage  ( A , ’ input . img ’ ) 

(*  set  the  filter  weights  to  1  *) 

Conf igTempMooreConst  (F, 3,3,1) 

(*  sum  the  gray  values  in  the  neighborhood  *) 

A  :  =  A  (  +  )  F 

(*  divide  the  sum  by  the  number  of  neighborhood  pixels  *) 

(*  and  place  this  value  into  the  resulting  image  *) 

A  :  =  A  /  9 

(*  write  the  resulting  image  to  OUTPUT. IMG  *) 

Putlmage  { A ,’ output . img ’  ) 
end  . 

Figure  5.2:  Image  Algebra  Implementation  of  a  Mean  Filter 

compared  to  16  lines  of  code  in  the  PASCAL  routine.  That  is 
68%  less  coding  required  to  implement  the  same  routine  in  the 
image  algebra.  Further,  the  object  code  of  the  IA  program 
requires  only  1 . 5K  bytes  of  storage  compared  to  the  2 . 5K  bytes 
of  storage  required  by  the  object  code  of  the  direct  PASCAL 
implementation.  Additionally,  if  the  user  desires  to  alter 
the  size  of  the  neighborhood  of  each  pixel,  the  PASCAL  program 
requires  recoding  of  both  the  filter  definition  and  the  limits 
of  program  iteration,  but  the  I A  program  requires  only  a 
change  of  the  dimensions  used  in  the  template  definition. 

Subsequent  inspection  of  the  executable  code  size  and  the 
execution  time  of  these  two  programs  shows  the  disadvantages 


of  this  version  of  the  image  algebra.  The  executable  code 
size  of  the  IA  program  needs  18. 5K  bytes  of  storage  compared 
to  only  3K  bytes  of  storage  for  the  PASCAL  program.  The 
execution  time  on  a  256x256  image  is  103.41  seconds  for  the  IA 
program  and  92.52  seconds  for  the  PASCAL  program. 

Median  Filter 

The  median  filter  is  a  nonlinear  operation  used  for  noise 
suppression  in  images.  The  filter  modifies  each  image  pixel 
to  reflect  the  median  gray  value  within  a  small  neighborhood 
of  each  image  pixel  where  the  neighborhood  is  defined  by  the 
configuration  of  the  filter.  This  filter  reduces  the  noise  in 
an  image.,  and  it  usually  does  not  affect  step  or  ramp  edges 
(6:330-331  )  . 

Figure  5.3  shows  a  PASCAL  implementation  of  the  median 
filter  with  the  input/output  routines  necessary  to  make  the 
program  into  an  executable  module.  The  algorithm  for  this 
filter  is  straightforward.  Following  the  declaration 
sections,  the  routine  reads  an  image  from  a  user-specified 
external  described  in  the  mean  filter  above.  After  obtaining 
the  input  image,  the  program  steps  through  each  pixel  of  the 
image.  At  each  pixel  of  the  image,  designated  by  A(i,j),  all 
of  the  pixel  gray  values  lying  under  the  filter  configuration 
centered  over  A ( i , j )  are  put  into  a  histogram.  The  program 
searches  for  the  median  entry  in  the  histogram  and  places  the 
gray  value  at  this  location  into  the  image  at  pixel  A(i,j). 


5-5 


program  Medi an_Fi 1  ter  ( input , output )  ; 

(*  Let  image  A  be  a  two  dimensional,  integer  valued  *) 
(*  array  of  dimensions  Image_Rows  by  Image_Colurans .  *) 
(*  Let  filter  F  be  a  3x3  array  of  constant  weight  =  1.  *) 
(*  Let  H  contain  the  histogram  of  the  neighborhood  *) 
{*  assuming  a  maximum  of  32  gray  levels  (0-31).  *) 
(*  The  result  of  the  filter  operation  is  image  R.  *) 


type  Image  =  array  ( . 1 . . 256 , 1 . . 256 . }  of  integer; 

Template  =  array  (  .  1  .  .  3  ,  1  .  .  3  .  )  of  integer; 
var  A,R  :  Image; 

F  :  Template; 

H  :  array  (.0..31.)  of  integer:  (*  histogram  *) 

Iraage_Rows , Image_Columns , i , j , k , 1  :  integer; 

number , med_val  :  integer; 

XInclude  ’IO_PAS.PAS' 
begin  (*  program  Median_Fi Iter  *) 
reset  ( input ) ; 
rewrite  (output); 

(*  read  image  A  from  an  external  file  *) 

Get  Image  (A,  ’input.img’  , Image_Rows , I  mage  Co lurans )  ; 

(*  set  the  filter  weights  to  1  * ) 
for  i  :=  1  to  3  do 
for  j  :=  1  to  3  do  F(.i,j.)  :  =  1; 

(*  for  each  pixel  in  image  A  *) 
for  i  : -  1  to  Image_Rows  do 
for  j  : =  1  to  ImageColumns  do  begin 
(*  build  a  histogram  of  the  neighborhood  *) 
for  k  :=  0  to  31  do  H( .k. )  :=  0; 
for  k  :=  -1  to  1  do 
for  1  :=  -1  to  1  do 

if  (i+k>=l)  and  ( i +k < = Image_Rows )  and 
(j+l>=l)  and  ( j + 1 < = ImageColumns ) 
then  H (  . A (  . i +k  ,  j  +  1  .  ) .  )  :=  H (  . A (  . i +k , j  + 1  .  )  .  )  +  1; 

(*  find  the  number  of  non-zero  histogram  values  *) 

number  : =  0 ; 

for  k  :=  1  to  31  do 

if  H(  .k.  )  <>  0  then  number  : =  number  +  1; 

(*  the  median  gray  value  is  the  median  entry  in  H  * ) 
med_val  :  =  round  (number/2); 

(*  find  and  place  the  median  gray  value  nto  image  R  *) 
number  : =  0 ; 

for  k  : =  1  to  31  do  begin 
if  H (  . k  .  )  <  >  0  then  number  : -  number  +  1 ; 
if  number  =  med_val  then  R(.i,j.)  :=  H(.k.); 
end ;  ( *  for  k  * ) 

end ;  ( *  for  j  * ) 

(*  write  the  resulting  image  to  an  external  file  *) 
Putlmage  (R, 'output. i mg’ , Image  Rows, Image _Columns) ; 
end.  (*  program  Med i an_F i 1 te r  *) 

Figure  5.3:  PASCAL  Implementation  of  a  Median  Filter 


Median_F i 1  ter 

(*  Let  image  A  be  a  two  dimensional,  integer  valued  array  *) 
(*  of  dimensions  Image_Rows  by  Image_Columns .  *) 

(*  Let  filter  F  be  a  3x3  array  of  constant  weight  =1.  *) 

(*  The  result  of  the  filter  operation  is  in  image  R.  *) 

var  il,A,B,C,R  :  image 
F  :  template 
x,dot  :  integer 
begin 

(*  read  image  A  from  an  external  file  *) 

Getlmage  ( A , ’ input . img ’ ) 

(*  initialize  F=1  *) 

Conf igTempMooreConst  (F,3,3,l) 

(*  initialize  11=1,  R=0  the  same  size  as  A  *) 

Constlmage  ( i 1 , A . row , A . co 1 , 1 ) 

Constlmage  ( R , A . row , A . co 1 , 0 ) 

(*  calculate  the  median  frequency  threshold  *) 

C  :=  { i  1  (  +  )  F}  /  2 

for  x  :=  Mi nVal Image  {A}  to  MaxVallmage  {A}  do  begin 
dot  : =  { A=x ) . i 1 

if  dot  >  0  then  begin 

(*  calculate  the  frequency  of  image  pixels  *) 

(*  under  template  F  with  gray  value  >=  x  *) 

B  :=  { A  >  =x }  <  +  )  F 

(*  update  the  pixel  gray  value  if  frequency  is  *) 

(*  greater  than  the  median  frequency  threshold  *) 

R  : =  R  v  {  { {B-C} >0} *x) 
end  (*  if  dot  then  *) 
end  (*  for  x  do*) 

(*  write  the  resulting  image  to  an  external  image  *) 
Putlmage  ( R output . img ’ ) 
end . 

Figure  5.4:  Image  Algebra  Implementation  of  a  Median  Filter 


Figure  5.4  shows  the  median  filter  algorithm  written  in 
the  image  algebra  language.  This  operation  of  this  algorithm 
is  not  as  obvious  as  the  direct  PASCAL  implementation.  After 
the  usual  declarations  and  obtaining  the  input  image,  A,  this 
algorithm  initializes  three  more  images  of  identical  size  to 
A.  Image  il,  which  is  used  by  the  dot  product  operation  to 
count  pixels,  is  set  to  one  at  every  pixel.  Image  C  is  set  to 
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the  median  frequency  threshold  at  each  pixel.  This  means  each 
pixel  in  C  assumes  a  value  corresponding  to  the  median  number 
of  pixels  within  the  neighborhood  of  that  image  pixel.  Image 
R,  the  result  image,  is  initialized  to  zero  at  every  pixel. 

The  remainder  of  the  algorithm  is  iterated  by  gray  level 
from  the  minimum  gray  value  to  the  maximum  gray  value  in  image 
A.  During  each  iteration,  image  B  is  set  to  the  frequency  of 
gray  levels  within  the  neighborhood  of  each  pixel  greater  than 
the  present  gray  level  of  the  iteration.  Then,  each  pixel 
that  has  a  frequency  greater  than  the  median  frequency 
threshold  is  updated  with  the  present  gray  level  of  the 
iteration.  Thus,  each  pixel  is  updated  with  the  a  new  gray 
value  until  its  frequency  is  greater  than  or  equal  to  the 
median  frequency  threshold. 

A  comparison  of  the  Figures  5.3  and  5.4  demonstrates  the 
advantages  of  the  image  algebra  description.  As  shown,  the 
image  algebra  version  contains  13  lines  of  executable  (non¬ 
comment)  code  compared  to  24  lines  in  the  PASCAL  version.  The 
result  is  better  than  45%  less  programming  required  by  the 
image  algebra  implementation.  Additionally,  the  object  code 
of  the  IA  program  requires  only  2 . 5K  bytes  of  storage  compared 
to  the  3K  bytes  of  storage  required  by  the  object  code  of  the 
direct  PASCAL  implementation.  As  in  the  IA  implementation  of 
the  mean  filter,  the  alteration  of  the  pixel  neighborhood  in 


the  IA  version  of  the  median  filter  requires  only  a  change  in 
the  dimensions  of  the  template  definition. 

Further  inspection  of  the  executable  code  size  and  the 
execution  time  of  these  two  programs  shows  the  disadvantages 
of  this  version  of  the  image  algebra.  The  executable  code 
size  of  the  1A  program  requires  19K  bytes  of  storage  compared 
to  only  3K  bytes  of  storage  required  by  the  PASCAL  program. 

The  execution  time  on  a  256x256  image  is  231.36  seconds  for 
the  IA  program  and  167.16  seconds  for  the  PASCAL  program. 

Local  Mode  Filter 

The  local  mode  filter  is  a  nonlinear  operation  useful  for 
noise  suppression  in  images  also.  The  filter  modifies  each 
image  pixel  to  reflect  the  most  frequent  (mode)  gray  value 
within  a  small  neighborhood  of  each  image  pixel  where  the 
neighborhood  is  defined  by  the  configuration  of  the  filter. 
This  filter  can  remove  noise  from  an  image  without  introducing 
errors  in  step  or  ramp  edges.  Furthermore,  it  is  usually  more 
responsive  to  local  image  context  than  the  median  filter.  For 
example,  assume  the  image  depicted  on  the  next  page  in  Figure 
5.5  is  to  be  filtered  and  the  image  pixel  under  modification 
is  in  the  center.  Without  additional  knowledge  of  the  image 
gray  values,  the  logical  choice  for  the  new  gray  value  would 
be  either  one  or  seven.  The  local  mode  filter  will  chose  one 
of  those  values.  However,  the  median  filter  will  chose  three, 
essentially  ignoring  the  context  of  the  image  neighborhood. 


Figure  5.5:  Image  for  Local  Mode  Filter  Example 


Figure  5.6  shows  a  PASCAL  implementation  of  the  local 
mode  filter  with  the  input/output  routines  necessary  to  make 
the  program  into  an  executable  module.  The  algorithm  for  this 
filter  is  nearly  identical  to  the  median  filter.  The  only 
difference  is  that  the  image  pixel  gray  value  is  updated  with 
the  mode,  rather  than  the  median,  of  the  histogram. 

Figure  5.7  shows  the  same  algorithm  written  in  the  image 
algebra.  This  implementation  is  similar  to  the  image  algebra 
description  of  the  median  filter.  After  the  declarations  and 
obtaining  the  input  image,  A,  this  algorithm  initializes  three 
images  to  zero:  B,  which  is  used  to  hold  the  frequency  of 
pixels  within  each  neighborhood  at  a  gray  level;  C,  which  is 
used  to  hold  the  maximum  frequency  of  pixels  within  each 
neighborhood  at  a  gray  level;  and  R,  which  is  used  as  the 
resultant  image.  Image  C  is  set  to  the  median  frequency 
threshold  at  each  pixel. 

The  remainder  of  the  algorithm  is  iterated  by  gray  level 
from  the  minimum  gray  value  to  the  maximum  gray  value  in  image 
A.  During  each  iteration,  image  B  is  set  to  the  frequency  of 
gray  levels  within  the  neighborhood  of  each  pixel 
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program  Local_Mode_Fi 1  ter  ( input , output ) ; 

(*  Let  image  A  be  a  two  dimensional,  integer  valued  *) 

(*  array  of  dimensions  Image_Rows  by  Image_Columns .  *) 

(*  Let  filter  F  be  a  3x3  array  of  constant  weight  =  1.  *) 

(*  Let  H  contain  the  histogram  of  the  neighborhood  *) 

(*  assuming  a  maximum  of  32  gray  levels  (0-31).  *) 

(*  The  result  of  the  filter  operation  is  image  R.  *) 

type  Image  =  array  (  . 1 . . 256 , 1  . .  256  .  )  of  integer; 

Template  =  array  (  .  1  .  .  3  ,  1  .  .  3  .  )  of  integer; 
var  A,R  :  Image; 

F  :  Template ; 


H  :  array  (.0..31.)  of  integer:  (*  histogram  *) 

Image_Rows , Iraage_Co 1 urons , i , j , k , 1  :  integer; 

max_val  :  integer; 

XInclude  ’ I0_PAS.PAS’ 

begin  (*  program  Local_Mode_Fi Iter  *) 
reset  ( input ) ; 
rewrite  (output); 

(*  read  image  A  from  an  external  file  *) 

Get  Image  (A,  ’input.img’  , I mage  Rows , Image_Colurans )  ; 

(*  set  the  filter  weights  to  1  *) 
for  i  :=  1  to  3  do 
for  j  :=  1  to  3  do  F(  .  i  ,  j  .  )  :=  1; 

( *  for  each  pixel  in  image  A  *  ) 
for  x  : =  1  to  ImageRows  do 
for  j  : =  1  to  Image_Columns  do  begin 
(*  build  a  histogram  of  the  neighborhood  *) 
for  k  :  =  0  to  31  do  H(.k.)  :=  0; 

for  k  : =  - 1  to  1  do 
for  1  :=  -1  to  1  do 

if  (i+k>=l)  and  ( i +k <  =  Image  Rows  )  and 
(j+l>=l)  and  ( j + 1 < = Image_Columns ) 
then  H (  . A (  . i +  k , j  + 1  .  )  .  )  :=  H (  . A (  . i +k  ,  j  +  1 .  )  .  )  +  1; 

(*  find  the  mode  of  the  histogram  *) 
maxval  : =  H (  . A (  .  i  ,  j  .  )  .  )  ; 
for  k  :=  1  to  31  do 

if  H(.k.)  >  maxval  then  max_val  :=  H(.k.); 

(*  find  and  place  the  mode’s  gray  value  into  image  R  * ) 
for  k  :=  1  to  31  do 

if  H(.k.)  =  maxval  then  R(.i,j.)  :=  H(.k.); 

end ;  ( *  for  j  * ) 

(*  write  the  resulting  image  an  external  file  *) 

Put  Image  (R,  'output. img’  , Image _Rows, I mage_Co 1 urans )  ; 
end.  (*  program  Local_Mode_Fi 1  ter  *) 

Figure  5.6:  PASCAL  Implementation  of  a  Local  Mode  Filter 


Local_Mode_F i 1  ter 

(*  Let  image  A  be  a  two  dimensional,  integer  valued  *) 

(*  array  of  dimensions  Image_Rows  by  Image_Columns .  *) 

(*  Let  filter  F  be  a  3x3  array  of  constant  weight  =  1.  *) 

(*  The  result  of  the  filter  operation  is  image  R.  *) 

var  A,B,C,R  :  image 
F  :  template 
x,dot  :  integer 
beg  i  n 

(*  read  image  A  from  an  external  file  *) 

Getlmage  ( A ,  ’ input . img '  ) 

(*  initialize  the  filter  F=1  *) 

Conf igTempMooreConst  (F,3,3,l) 

(*  initialize  B  =  C  =  R  =  0  the  same  size  as  A  *) 
Constlmage  ( R , A . row , A . col , 0 } 

B  :  =  C  :  =  R 

for  x  :=  MinVallmage  (A)  to  MaxVallmage  {A}  do  begin 
B  : =  { A  =  x )  (  +  )  R 

R  : =  R  v  {  {  { B vC-C | >  0 } *  x } 

C  :  =  B  v  C 
end  ( *  for  x  *  ) 

(*  write  the  resulting  image  to  an  external  file  *) 
Putlmage  ( R output . img ’ ) 
end  . 

Figure  5.7:  Image  Algebra  Implementation 
of  a  Local  Mode  Filter 


equal  to  the  present  gray  level  of  the  iteration.  Then,  each 
pixel  having  a  frequency  greater  than  the  present  maximum 
frequency  of  that  pixel,  held  by  C,  is  updated  with  the 
present  gray  level  of  the  iteration.  Thus,  each  pixel  is 
updated  with  the  a  new  gray  value  as  long  as  its  frequency  is 
greater  than  the  present  maximum  frequency  for  that  pixel. 

The  final  step  of  each  iteration  updates  the  maximum  frequency 
at  each  p l xe 1  . 

Once  again,  a  comparison  of  the  two  implementations  shows 
the  programming  advantages  of  the  image  algebra.  As  shown, 
the  image  algebra  version  contains  10  lines  of  executable 


(non-comment)  code  compared  to  20  lines  in  the  PASCAL  version. 
This  translates  to  50%  less  programming  required  by  the  image 


algebra  implementation.  Further,  the  object  code  of  the  IA 
program  needs  2K  bytes  of  storage  compared  to  the  3K  bytes  of 
storage  required  by  the  object  code  of  the  direct  PASCAL 
implementation.  Like  the  other  two  algorithms,  changing  the 
size  and  configuration  of  the  neighborhood  is  easier  with  the 
IA  version,  too. 

Inspection  of  the  executable  code  size  and  the  execution 
time  of  these  two  programs  shows  the  disadvantages  of  this 
version  of  the  image  algebra.  The  executable  code  size  of  the 
IA  program  requires  18. 5K  bytes  of  storage  compared  to  only  3K 
bytes  of  storage  needed  by  the  PASCAL  program.  The  execution 
time  on  a  256x256  image  is  203.82  seconds  for  the  IA  program 
and  132.22  seconds  for  the  PASCAL  program. 

Image  Algebra  vs  High  Level  Language  Routines 

As  can  be  seen  by  the  previous  examples,  the  power  and 
simplicity  of  the  image  algebra  is  tremendous.  The  image 
algebra  provides  a  concise  notation  for  algorithm  development 
that  removes  the  engineer’s  attention  to  details  required  when 
programming  in  PASCAL  or  some  other  high  level  language.  This 
notation  allows  the  image  processing  engineer  to  concentrate 
on  building  image  processing  tools  rather  than  on  writing  and 
debugging  computer  programs. 


mograra 

source 

object 

execute 

16x16 

256x256 

mean 

1  .  5 

2 . 5 

3 . 0 

2.03 

7 .44 

92 . 52 

median 

2 . 5 

3.0 

3 . 0 

2 . 38 

12.14 

167.16 

mode 

2.0 

3 . 0 

3 . 0 

2 .25 

9.95 

132.22 

mean 

1 . 0 

1  .  5 

18.5 

5 .85 

11.63 

103.41 

median 

1  .  5 

2 . 5 

19 . 0 

30. 12 

41.52 

231.36 

mode 

1 .0 

2.0 

18.5 

29.41 

39.39 

203.82 

Table  5.1:  Code  Size  and  Execution  Time  for  PASCAL  and  Image 

Algebra  Example  Algorithms 


Futhermore,  accommodating  arbitrary  neighborhood  sizes 
and  configurations  is  much  easier  to  program  into  the  image 
algebra:  merely  alter  the  configuration  of  the  templates. 

Even  though  the  IA  provides  increased  programming 
productivity  for  image  processing  routines,  the  size  of  the 
executable  code  and  slower  execution  times  are  significant 
detriments  to  this  version  of  the  AFITIA.  Table  5.1  shows  a 
summary  of  the  storage  requirements  and  execution  times  for 
the  two  versions  of  each  image  processing  algorithm.  The  top 
set  of  data  corresponds  to  the  PASCAL  implementation  of  each 
algorithm,  and  the  bottom  set  of  data  corresponds  to  the  image 
algebra  implementation.  The  table  shows  that  both  the  source 
code  and  object  code  sizes  are  smaller  in  the  IA  programs,  but 
the  executable  code  size  is  much  larger.  This  is  due  to  the 
linking  of  the  entire  library  of  IA  operations  to  each  IA 
module  whether  or  not  each  operation  is  needed. 
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The  execution  times  on  16x16,  64x64,  and  256x256  binary 
images  show  that  the  IA  implementations  are  much  slower  than 
the  direct  PASCAL  implementations.  The  execution  times  of  the 
PASCAL  implementations  display  the  expected  na  increases  in 
execution  time  where  n  is  the  image  dimension.  Even  though 

Table  5.1  does  not  display  na  increases  in  execution  time  for 

the  IA  programs,  this  does  not  imply  that  the  IA  algorithms 
are  more  efficient  for  larger  images.  The  increased  execution 

time  due  to  the  larger  images  is  hidden  by  the  high  execution 

overhead  from  all  of  the  procedure  and  function  calls. 

The  AFIT  image  algebra  operations  and  preprocessor  should 
be  modified  to  reduce  both  of  these  problems.  One  method  of 
reducing  the  code  size  is  to  program  the  preprocessor  so  that 
only  those  basic  IA  operations  used  by  the  program  are  linked 
to  the  object  code.  Reducing  the  execution  time  of  the  IA 
programs  will  be  more  difficult  due  the  overhead  involved  in 
procedure  and  function  calls. 

As  was  shown  in  the  previous  examples,  the  implementation 
of  an  image  processing  algorithm  in  the  image  algebra  can  be 
very  different  than  its  direct  implementation  in  PASCAL.  The 
algorithm  must  be  formulated  differently  for  the  two  languages 
because  of  the  difference  in  elemental  operations.  Because 
the  IA  forces  an  image  processing  engineer  to  approach  the 
programming  from  a  new  and  very  different  viewpoint,  it  is 
difficult  to  determine  which  method  is  easier  to  implement  at 


VI .  Observations  and  Recommendations 

The  image  processing  environment  designed  and  implemented 
in  this  project  is  a  useful  and  highly  flexible  development 
tool.  Through  the  use  of  a  preprocessing  program,  it  is 
capable  of  translating  image  processing  algorithms,  written 
with  image  algebra  operators  and  PASCAL  control  structures, 
into  executable  programs.  The  entire  process  can  be  automated 
through  the  use  of  a  VAX  command  file. 

The  AFIT  image  algebra  language  and  preprocessor  allow 
the  construction  of  image  processing  algorithms  with  minimal 
development  time  and  programming  effort.  The  actual  program 
code  can  be  reduced  by  45-70%,  and  the  low  level  programming 
details  of  each  algorithm  are  managed  by  the  IA  language  and 
preprocessor.  The  AFITIA  supports  real  and  integer  valued 
images  and  templates,  and  the  syntax  of  the  IA  source  file  is 
relatively  simple  and  flexible.  The  AFITIA  support  of  all 
PASCAL  control  structures  provides  an  IA  language  capable  of 
performing  any  i mage- to- i mage  transformation. 

One  major  problem  with  this  design  appeared  and  remains 
unresolved:  complete  software  independence.  The  dependencies 
of  the  preprocessor  target  language  have  not  been  removed  from 
the  AFITIA  language.  Since  the  ultimate  IA  language  will  be 
translated  to  many  different  computer  languages,  the  image 
algebra  language  must  define  and  support  program  control 


statements  common  to  the  more  popular  high  order  languages. 
Further  research  will  provide  valuable  guidance  about  which 
statements  should  be  supported  and  the  final  form  of  the  image 
algebra  language. 

Some  other  shortcomings  of  this  version  of  the  AFITIA  are 
a  lack  of  support  for  complex  images  and  templates,  the  large 
size  of  the  executable  code  for  eacn  image  algebra  algorithm, 
and  their  slower  execution.  The  correction  of  these  problems 
through  a  more  intelligent  linking  procedure  and  faster  code 
for  the  IA  operations  is  highly  recommended. 

Another  topic  of  interest  stemming  from  this  project  is 
the  inherent  parallelism  of  the  image  algebra  operators. 
Subsequent  research  could  determine  which  operations  are 
amenable  to  parallel  execution  in  either  vector  (1-D  array)  or 
image  ( 2-D  array)  processing  architectures. 

Overall,  the  blending  of  the  image  algebra  developed  at 
the  University  of  Florida  with  the  AFIT  preprocessor  makes 
significant  progress  toward  a  simple,  powerful,  and  universal 
image  processing  tool. 


Appendix  A:  AFIT  Image  Algebra  Operations 


[«a«i  roaataK  '  ri  iop«r.«Bv' )] 

Module  Real  _l  A.operati  oaa; 

( ■■■■immmiiiimtioiimiiummiiimiiaiimmmiiuiuitimin 
(I  Tbi  *  file  coatai  at  all  of  t  be  ban  c  i  mage  and  tempi  ate  operations  fro*  t) 
(t  t  be  I  mage  Algebra  developed  for  the  AFATL#  Egl  i  a  AFI,  FL.  Tbe  routines  s) 
(s  in  this  file  ere  written  for  real  valued  I  nafes  and  templates.  A  t) 

(s  number  of  AFIT  extensions  are  included.  t) 

(xxsxsxxxxxxsssxsttsxxxxsxtxssssssssxsxxstxsssssxsssssssstssssxssssstsissxsx) 


(s  tl  obal  constants* 

types,  and 

van  abl  es  *> 

const 

Maxi  aageRo w 

= 

22 ; 

<  X 

nail  aua  i  aage  row  di  aensi  on 

X) 

Maxi  aageCol  usd 

= 

32; 

(X 

aaxiaua  i  safe  col  uu  diaenston 

X) 

MaxTeapCel  1 

ioo; 

(X 

aaxi  aua  nuaber  of  teapl  ate  pi  xel 

s  X) 

Naae  ten 

= 

so; 

(  X 

aaxiaua  length  of  file  naaes 

X) 

t  ype 

I  aageType 

real  ; 

(  X 

data  type  of  i  aage  fray  values 

X) 

Teapl  ateType 

= 

real  ; 

(X 

data  type  of  teapl  ate  weights 

X) 

I  aaf e  =  record 

(X 

pascal  implementation  of  image  operand  x) 

row 

i nteeer; 

(X 

i  aage  row  di  aensi  on  t) 

col 

i nteger; 

(X 

i  aage  col  uao  di  aensi  on  t) 

gray 

array  ( 

1.  .  Maxi  mageRow,  l.  .  Maxi  mageCol  umm.  )  of 

I  mageType 

end; 

Tempi  ate  =  record 


nun 

cf  ( 


( t  i  mage  record  X) 

(s  PASCAL  i  *pl  ewetatk  os  of  template  operand  x) 
(X  number  of  pixels  wt  t  bi  n  tbe  co  nf  i  (ur  at  i  on  X) 
MaxTeapCel  I  .  )  of  record  (X  conf  i  (ur  at  i  on  X) 

x) 
x) 


X) 


Oper  St  r 
Name  St  r 


I  Rl,  I  R2,  i  R3i  i  R4,  l  PS 
t  Rl  ,  l  R2  ,  t  R3,  l  R4,  t  RS 
dot  vai  •  mi  nval  .  maxval 


i nteger; 
array  ( .  l. 

r  :  integer;  (x  row  offset  of  pixel 

c  :  integer;  (X  col  uma  offset  of  pixel 

w  :  Tempi  atetype;  (X  vwight  of  pixel 
end;  (X  array  record  X) 

end;  (X  template  record  X) 

vary!  ng(  3 ]  of  c bar ; 
varyi  ng(  Name  Leo]  of  ebar; 

Image;  (X  intermediate  image  operands  s) 

Template;  (X  intermediate  template  operands  s) 
ImageType;  (s  intermediate  scalar  operands  x) 


(isstituiitmiitsstisttitsttissiittsiiitititmttstttftttitiimitittititt) 

(B  UNARY  and  BINARY  I  HAGE  OPERATIONS  S) 

( tiiisiiiitttmiisiismtistmtsimiiitsstttttsitisttttmtttttttiiiitiit) 


function  iaageEqual  (  A,  B  :  Image)  :  boolean; 

(a  return  true  if  A  =  B,  point  wise  t) 
var  i  i j  :  i  nt  ege  r ; 

equal  bool  ean; 
begi  n  (s  function  ImageEqual  s) 

if  (  A.  row*  B.  row)  and  (  A.  col  *  B.  col  )  then  equal  true  else  equal  :*  false; 

»  :  *  I; 

while  (equal  =  true)  and  (i  <  =  a.  row)  do  begin 

j  :  =  i; 

while  (equai  =  true)  and  (j  <-  A.  col  )  do  begin 
if  A.  gr  ay( .  i  ,  j .  )  <>  B.  gra  y< .  i  ,  j  .  )  then  equal  :*  false; 
j  :  -  j  *  l  i 

end;  ( t  whi  I  e  j  a ) 
i  :  =  i  ♦  1; 

end;  (a  whi  lei  a ) 

ImageCqual  =  equal; 

end;  (a  function  ImageEqual  a) 

procedure  Constlmage  (  var  A  :  I  sage;  row,  col  integer;  value  :  ImegeType), 
(a  initialize  i  nag*  A  to  a  coustaui  value:  a  j  value  a) 
var  i  • )  :  integer; 
begin  (a  procedure  Constlange  a) 

A.  row:*  row;  A.  col  :  *  col; 
for  i  :  «  1  to  row  do 

for  j  1  to  col  do  A.gray(.i,j.)  value; 

end;  (a  procedure  Coiitluafe  a) 

function  HaiValtmage  (A  Image)  ImageType; 

(a  return  the  ain  ms  pi  zel  value  >o  t  be  »  nafe  a) 
var  i , j  integer; 

nil  val  I  Mf  e  T  yp#  , 
begin  (a  procedure  HaiVallmege  s) 
nil  val  A.  gra  y(.  1,1.). 

for  •  I  to  A.  row  do 

for  j  l  to  A.  rol  do 

if  A.grayf.t'j.)  »  seival  than  nanii  A-gray(  i.j-), 

Ha  i  v  a  I  l  aa(*  naival  . 


end.  (t  pro<»dure  NaiVil  t  aafr  I  ) 


•  /  V 


fuactioa  M  >«■!  ia>(<  (A  :  Iaage)  :  laageType; 

(t  rttan  tht  m  ni  bub  pi  i«l  value  11  iki  i  aage  f) 
var  ■  > i  :  i oteger; 

■  aval  laageType; 
bt(i  a  («  procedure  M  aval  Iaage  *) 

■  aval  :  ■  A.  gray!  .1,1.); 
for  i  :  ■  1  to  A.  row  do 
for  i  :  >  1  to  A.  col  do 

if  A.  g  r  a y(  .  i  ,  j .  )  <  ■  aval  tboa  ■  aval  :*  A.  gray! .  i  ,  J .  ) ; 

M  aval  1  aago  :  =  ■  aval  ; 

and;  (<  procodaro  M  aval  Iaage  t) 

procedure  AbsVal  Iaage  (  A  :  luge;  var  C  :  laage); 

(a  absolute  value  of  t  be  iaage:  c  -  I  A!  a) 

var  i  i  j  iateger; 

begi  a  (a  procedure  Absval  Iaage  a) 

C.  row  :  s  a.  row;  C.  col  *  A.  col  ; 
for  i  1  to  A.  row  do 

for  j  :  «  1  to  A.  col  do  c.  gray! .  i  ,  j  .  )  :  =  aba  (A.  gray!. i  ,  j<  )); 
ead;  (a  procedure  Absval  Iaage  a) 

procedure  laageAdd  <  A,  a  :  Iaage;  var  C  :  laage); 

(a  poiatMse  iaage  additloa:  C  «  A  ♦  •  a) 
var  i  ,  i  :  i stager; 
begi  a  (a  procedure  laageAdd  a) 
i  f  A.  row  >  :  I.  row  t  bea  C.  row  :  »  A.  row  el  ee  c.  row  :  •  b.  row, 

if  A.  col  >  =  1.  col  tbea  C.  col  :  ■  A.  col  else  c.  col  :<  >.  col  ; 

for  i  I  to  C.  row  do 

for  1  :•  I  to  C.  col  do  begi  ■ 

if  (i  <>  A.  row)  aad  (■  <  =  i.  row)  aad  (j  <>  A.  col  )  aad  (j  <»  l.col) 

thea  C.  gray!  .  ,  ,  i .  )  :■  A.  gr  a  y<  .  i  ,  j  .  )  «  l.  gr  ay(  .  i  ,  > .  ) ; 

If  (i  >  A.  row)  or  (J  >  A.  col  )  tbea  C.  grtyl  .  i  ,  ).  )  I.  gray!  .  i  ,  j  .  ) ; 

if  (■  >  a.  row)  or  (j  >  l.col)  tbea  C.  gray!  .  i  ,  ).  )  A.  gray!  .  i  ,  j  .  ) ; 

ead;  (a  for  j  da  a) 
ead;  (a  procedure  laageAdd  a) 

procedure  1  aage tabt r ac t  1  A.  ■  luge;  var  C  :  Iaage), 

(a  poiatwtee  iaage  subtraction  c  A  -  1  a) 
var  i,j  iateger. 

begin  (I  procedure  I  aage tubt r ac t  a) 
if  a.  row  >  -  a.  row  tbea  c.row  A.  row  else  c.row  :  s  s.row, 

if  A.  col  >  ■  ■.  col  tbea  C.rol  A.  col  elae  C.col  :  *  l.col. 

fori  I  to  c.row  do 

f  o  r  j  I  t  o  c.  i  o  i  do  beg  t  a 

if  <i  <  '  A.  row)  and  I  i  (  =  1.  row)  aad  (a  «  »  A. col)  and  Ij  <  -  l.col) 

tbea  C. gray!. i, a.)  A.grayf.i.j.)  -  1.  gray!.  i,).>; 

if  (i  >  A.  row)  or  (j  >  A. col)  tbea  C.greyl.i.j.)  :  =  -l.grarl.i.j.l, 

if  I  i  >  l.row)  or  (j  >  l.col)  tbea  c.grart.i.j.)  A.  grayi  .■.).), 

ead.  la  for  j  do  a) 
ead.  (a  procedure  I  aatetubl  rirt  a) 


procedure  I  mageMul  tt  pi  y  (  A,  •  :  Image;  var  C  :  Image); 

(t  point  Mae  i  aa|e  mu  I  tl  pi  i  cati  on:  C  *  Ax  pa) 
ear  I , J  :  i nt eg er ; 

b*(i  g  (a  procedure  I mageMul  tt  pi  y  a) 
if  A.  row  >  =  >.  row  then  C.  row  :>  A.  row  el  aa  C.  row  a.  row; 
if  A.  col  ><  a.  col  then  c.col  :*  A. col  el  ae  C.col  :»  a. col; 
for  i  :  *  1  to  C.  row  do 
for  j  :  >  l  to  C.col  do  bafi  a 

if  (i  <=  A.  raw)  and  (i  <•  B.  row)  and  (j  <•  A.  col)  and  (j  <*  a. col) 
then  C.  gray<  .  i  ,  i.  )  :>  A.  gray( .  i  ,  j.  )  a  a.  lr  ay( .  i  ,  J .  ) 
cl  ae  C.  a r ay( .  i  ,  j.  )  :  «  o; 
end;  (a  for  i  do  a) 
end;  (a  procedure  I  mageMul  ti  pi  y  a) 

procedure  ImageDivide  (  A.  a  :  Image;  var  C  :  lata(e); 

(a  point  Mae  i  eafe  di  vi  ei  on:  C  =  A  /  B  a) 
var  i >  j  :  i nteger; 
begin  (a  procedure  ImageDivide  a) 
if  A.  row  >  =  B.  row  then  C.  row  :  »  A.  row  elae  C.  row  :  -  B.  row, 

if  A.  col  >:  a.  col  then  c.col  A.  col  else  C.col  :  =  a. col; 

for  i  :  =  1  to  C.  row  do 
for  j  :  »  l  to  c.  col  do 

if  (i  <  =  A.  row)  and  (i  :=  a.  row)  and  (j  <>  A.  col  )  and 
(j  <  =  a. col)  and  (  B.  gray( .  i  ,  j  .  )  <>  0) 

then  C.  gray< .  i  ,  j. )  :  *  A,  g  r  a  y(  .  i  ,  j .  )  /  a.  gray! .  i  ,  j .  ) 

el  so  c.  gray( .  I  ,  j .  )  :  *  o; 

end;  (a  procedure  ImageDivide  a) 

procedure  I  mage  Exponent  <  A.  B  :  Image;  var  C  :  Image); 

(a  point  Mae  exponent!  ati  on  of  an  image  by  aa  image:  C  :  =  Ataa  »  exp(Baln(A))  a) 
var  i  ,  j , k  :  i at  eger ; 

begin  (a  procedure  I  mageEiponent  a) 

C.  row  :  =  A.  row,  c.  col  :  =  A.  col  ; 
for  I  1  to  C.  row  do 
for  j  :  =  1  to  c.  col  do 

if  A. gray!. i, j.)  <>  0  then  C.  grayf .  i  ,  ) .  )  :*  1  n(  abet  A.  gr  a  y(  .  i  ,  )  .  ) ) ) 

elae  C.  gr  a  y(  .  i  ,  ;  .  )  :  -  0; 

I  mageMul  1 1  pi  y  (  C,  a,  C)  ; 
for  i  1  to  c.  row  do 

for  j  :  =  i  to  C.col  do  c.gray(.i,j.)  :  =  exp(C.gray(.i,j.)); 
for  i  1  to  C.  row  do 
for  j  :  •  1  to  C.col  do  begin 
if  (A. gray!. i, j.)  <  0)  and  I  C.  grayl  .  i  ,  j.  )  (>  0) 

then  c.  grayl  .  i  ,  j.  )  :=  1  '  C.  grayl  .  i  ,  j  .  I ; 

if  A. grayl. i, j.)  -  o  then  c.  g  r  a  yl  .  i  ,  a .  )  :  =  o; 
end;  (a  for  )  do  a) 
end;  IS  procedure  I  mage  Eipooe  nt  a) 


procedure  laageMai  (  A.  8  :  lease;  »»r  C  :  !u|e); 


*>(••  A  aid  8:  C  *  aai  (A,  8)  t) 


be (iu  (*  procedure  laageMai  t) 
if  A.  row  n  8.  row  thea  C.  row  :  =  A.  row  el  ae  C.  row  :»  8.  row; 

if  A.  col  >  =  8.  col  thea  C.  col  :>  A.  col  cite  C.  col  :  =  8.  col  ; 

for  i  :  =  l  to  C.  row  do 
for  j  :  -  1  to  C.  col  do  be(i  a 

if  (i  <=  A.  row)  aad  (i  <*  8.  row)  aad  (j  <  -  A. col)  aad  (J  <  =  8.  col  ) 

thea  if  A.  gray( .  i  ,  j.  )  >=  8.  gray( .  i  ■  j.  )  tbea  C.  gray( .  i  ,  j  .  )  :*  A.  gray( .  i  ,  j  .  ) 

el  ae  C.  gray!  .  t  ,  j .  )  :=  8.  gray< .  i  ,  j  .  ) ; 

if  (i  >  A.  row)  or  (j  >  A.  col) 

then  if  8.  gray<  .  i  ,  j  .  )  >  0  thea  C.  (ray( .  i  ,  j  .  )  :*  8.  gray<  .  i  ,  j .  ) 

else  C.  grayf  .  i  ,  j.  )  :  =  o; 
if  (i  >  8.  row)  or  (j  >  8.  col) 

then  if  A.  gray(  .  i  ,  j.  )  >  0  then  C.  gr ay<  .  i  ,  j .  )  :=  A.  (ray(  .  i  ,  j .  ) 

el  se  C.  gray< .  i  ,  j .  )  :  *  o; 


end;  (s  for  j  do  s) 
end;  (s  procedure  laageMai  t) 

procedure  lugeMn  (  A,  8  :  I  wage;  ear  C  :  I  safe); 
(*  return  the  poiatwige  ■  ai  aaa  of  i  safes  A  and  I 


C  =  at  a  (  A,  8)  *) 


begin  (t  procedure  XaageMn  X) 
if  A.  row  >  s  8.  row  thea  C.  row  : 
if  A.  col  >*  8.  col  thea  C.  col  : 
for  i  :  =  i  to  C.  row  do 


A.  row  el  ae  C.  row  :  =  8.  row, 
A.  col  else  c.  col  :>  8.  col  ; 


for  j  :  =  1  to  C.  col  do  begl  a 

if  (i  <=  A.  row)  aad  (i  <=  8.  row)  aad  ()  <=  A.  col  )  aad  (j  <»  8.  col  ) 

then  if  A.  *ray<  .  i  ,  j.  )  <=  8.  grayf .  i  ,  i .  )  thea  C.  grayC .  i  ,  j.  )  :*  A.  *ray( .  i  ,  ) .  ) 

else  c.  *ray(  .  i  ,  j.  )  :*  8.  fray<  .  i  ,  j  .  ) ; 

if  ( i  >  »  row)  or  ( J  >  A.  col  ) 

thea  If  8.  *ray<  .  i  ,  j  .  )  <  o  thea  C.  (ray(  .  i  .  j  .  )  :=  8.  gray! .  i  ,  j  .  ) 


el  se  C.  gray(  . 


,  j. )  :  =  o; 


if  (i  >  8.  row)  or  (j  >  8.  col) 
then  if  A.  gray(  .  i  ,  j  .  )  <  0  thea  C.  g  r  a  y(  .  i  ,  j  .  ) 


thea  C.  gray(  .  i  ,  ;  .  )  :■  A.  gr  ay< .  i  ,  j .  ) 
el  se  c.  gray<  .>,).)  :  o; 


end;  (a  for  >  do  t) 
end;  (S  procedure  ImageMa  s) 


procedure  Cbarlaage  (  A  :  Inle;  operator  :  OperStr;  ■  :  laalt;  var  C  : 
(a  retort  a  bleary  iatle  wt  t  b  piiela>l  for  A  t  <■<*»*>>»>*  ]  ■  true  t) 
(t  Thl  a  procedure  doea  not  uae  i  mage  algebra  priaativea  i  a  order  to 
(a  reduce  t  be  waorr  aad  computational  requi  retaaat  a  of  the  routine 
var  it j  :  i ntefer; 
begin  (a  procedure  Cbarlaage  a) 

C.  row  :  =  A.  row;  C.  col  <  A.  col  ; 

If  operator  ■  *  <  *  tbea 
for  i  :  *  l  to  A.  row  do 
f 01  i  :  *  1  to  A.  col  da 

If  A.  grayf .  i  ,  j.  )  <  ».  grayf . i , j. )  then  C.  grayf .  i  ,  j.  )  :»  1 

el  ae  c.  gray!  .  I  t  j.  )  :  *  0; 


if  operator  =  ’  <»•  tbea 
for  i  :  s  I  to  A.  row  do 
for  j  :  »  l  to  A.  col  do 


A.  gray( .  i 

■  J.  )  <*  ••  grayf  .  1 

i  J  -  )  t  bea 

C.  grayf .  I 

1  j  •  ) 

:  *  1 

el  «• 

C.  grayf .  i 

,  i- ) 

:  =  0; 

f  operator  =  * 

» •  t  be  n 

for  i 

:  2  1  to 

A.  row  do 

for 

te¬ 

ll 

b* 

e* 

o 

A.  col  do 

»  r 

A.  grayf .  j 

•  j.)  •  1.  grayf.  i 

i  i  j.  )  tbea 

C.  grayf .  i 

.  j.  ) 

:  *  1 

•  1  •• 

C.  grayf  .  I 

.  J- ) 

:  2  Oi 

f  opi 

trator  *  * 

> » *  t  be  n 

for  i 

i  :  -  1  to 

A.  row  do 

for  J  :  »  1  to  A.  col  do 

if  A.  gray< .  i  ,  j.  )  >•  >.  grayf .  i ,  ) .  )  t  boa  c.  grayf .  I  ,  i .  )  :*  i 

el  ae  c.  gray!  .  I  ,  J.  )  :  *  o; 

i  f  operat  or  •  *  >*  t  bea 
for  I  ;  »  l  to  A.  row  do 
for  J  :  a  I  to  A.  col  do 

if  A.  gray<  .  I  ,  J  .  )  >  1. grayf. I  •>.)  t  bea  c.  gray!  .  I  ,  j  .  )  :>  I 

el  ae  c.  gray!  .  i  ,  )  .  )  :  *  0; 

end;  (a  procedure  Cbarlaage  a) 

procedure  IinageDot  (  A.  ■  :  laage;  var  Sua  :  ImageType); 

(a  dot  product  of  two  i  aagea:  reault  :«  A  (dot)  b  t) 
var  i  ,  j  ,  RowMai,  Col  Nat  :  integer; 
begin  (a  procedure  laagebot  a) 

if  A. row  >*  (.row  tbea  aowaai  :  A.  row  el  ae  eowetai  :  I.  row, 

if  A. col  >  »  a. col  tbea  Col  Mat  :  =  A. col  elae  Cel  Hat  :  »  ».  col  . 

Sua  :  0; 

for  i  i  to  PowMaa  do 

for  )  I  to  ColMai  do  begin 

if  (i  <  A.  row)  and  (  i  (  ■  a. row)  aad  (  j  <  >  A. toll  and  I  j  *  a.  roll 

tbea  Sua  :  *  Saa  t  A.gray(.i,j.)  a  I.  gran  .  ■  ,  j.  I , 

ead;  ( a  for  )  do  a) 


Beaul  t  :  1  Saa; 

ead;  (f  procedure  laagebot  a) 


(I 
( * 
(* 


UK  All  Y  tad  BINARY  TEMPLATE  OPERATIONS 


*) 

t) 

*) 


procedure  Coaf  i  gToapflCoast  (  vir  A  :  Toaplato;  Col  i  ittgtr;  v»|  ut  :  To  apt  at  oTy  po ) ; 
(i  imtializ  alt  pi  i»l  •  11  t  ha  I  i  Col  horizontal  toaplato  to  val  ao  t) 

(s  t  ho  coator  is  aisaatd  to  bo  t  bo  physical  coator  of  t  ho  configuration  s) 
var  j i cat  :  \  at ogo r ; 

bogi  a  (t  procodur*  Co  af  »  g  ToaplCo  01 1  s) 

(s  iaitializo  t  bo  toapl  ato  coaf  i  g  ur  at  i  o  a  to  zoro  t) 
for  j  :*  l  to  MaxToapcol  I  do  bogi  a 

A.  cf g( . j  .  ) .  r  :  *  o; 

A.  cf  g(  .  j  .  ) .  c  :  =  0  , 

A.  c  f  g(  .  j  .  ) .  w  :  «  o; 

rou,  (i  for  )  do  t) 
cat  :  s  0 ; 

for  j  :  *  -(Col  di  v  2 )  to  (Col  di  r  2 )  do  bo  g  i  a 
cat  :  *  cot  ♦  l ; 

A.cfg(.cat.).r  *  o; 

A.  c f  g(  .  cat .  ) .  c  :  *  )  ; 

A.  c  f  g(  .  c  at .  ) .  w  ■  «  val  ao; 

•  ad;  (t  for  j  do  S) 

A.  aaa  :  *  cat; 

oad;  (t  proroduro  Coaf  i  gToapRCoaat  >) 

procodaro  Coaf  i  gToapvco  ■■  t  (  oar  a  Toapl  ato,  Row  :  latogor;  val  ao  :  Toapl  at  oTypo )  ; 
(s  i  ai  h  al  i  to  ail  pi  tol  •  ia  t  bo  Row  t  1  vortical  toaplato  to  val  no  t) 

(•  i  bo  coator  ia  aoaoaod  to  bo  t  bo  physical  coator  of  t  bo  configuration  t) 

var  i  , c  at  i at ogor ; 

bogi  a  (•  procoduro  Co  af.  g  Toapico  as  t  t) 

<s  initial'*#  t  bo  toaplato  r  o  af  i  g  ur  at  i  o  a  to  ioro  «) 
for  i  i  to  PSaiToapCoi  l  do  b*gt  a 

a.  f  f  g  (.».).  r  o . 

A-  f  f  gl.  i  )  .  r  0  . 

A.  (  fgl  i  .  )  .  v  0. 

oad.  (t  for  i  do  •) 

»  at  o. 

for  i  i  Row  dt  v  2)  id  {  *ov  41  i  2)  do  bogi  a 

<  at  c  at  •  l  . 

A - <  f  g (  (  at.  )  r  »  ; 

a.  c  f  g (  r  at .  ) .  c  o. 

a.  c  f  gt  rat  )  w  val  uo 

•  ad.  i  •  for  )  do  ■) 

A  aua  <  at 

oad.  (•  proiodwro  >  o  a#  i  g  To  ap  v  «  «  a  •  t  •  > 


procoduro  Coaf  i  gToapMooroCoast  (vir  A:  Toapl  ato;  low,  col  :  i  itc(«r;  vil  at  :  Toapl  at oTypo ) 
{*  iiititlit*  all  pi  zol  •  ia  I  ha  Row  by  Col  Moor*  tcapl  at*  to  val  uo  t) 

(«  t  bo  coat  or  ia  assuatd  to  bo  t  bo  physical  coator  of  t  bo  coaf  i  gurati  oa  s) 
oar  i • j i cat  :  i at ogor ; 

bogi  a  (t  procedure  Coaf  i  gTo ■pNooroCoaat  a) 


( 1  1  ■■  » 1 

al  i  zo 

t  ho  toapl  ato 

coaf  1  (oral i  oa 

for  i  = 

l  to 

MaxToapCol  1 

do  bogi  a 

A.cfdf . 

i  .  ) .  r 

:  »  o; 

A.  cf  ((  . 

.  >.c 

:  *  0; 

A.  cf«(  . 

i  .  ).  « 

:  »  o; 

•  id;  (i 

col  :  •  0 

for 

l 

do  t) 

for  >  :»  -(Row  di  v  2)  to  (Row  di  v  2)  do 

for  j  :*  -(Col  di  »  2)  to  (Col  div  2)  do  bogia 
cat  :  =  cat  ♦  l ; 

A.  c f  g(  .  cat .  ) .  r  :  *  i  ; 

A.  cf  g(  .  cat .  ) .  c  :  *  j ; 

A.  cf  g(  .  cat .  ) .  w  :  -  val  uo; 

•  ad;  (  *  for  j  do  •) 

A.  iva  :  *  cat  ; 

•  ad;  (>  procoduro  Coaf  i  gToapHooroCoast  t) 


procodar#  Coaf  i  gToapVNCoaat  (oar  A  :  Toapl  ato;  rad  :  latogor;  val  no  :  Toapl  at  oTypo ) ; 
(s  laitiaiiio  a  voa  Nouasaa  toapl  ato  coaf  i  gar  at  i  oa  of  radius  rad  to  val  ao  a) 

(a  tbo  coator  is  assaaod  to  bo  t  bo  physical  coator  of  t  bo  coaf  i  gar  at  i  oa  t) 
var  i  ,  auabor  :  iatogor; 

bogi  a  («  procoduro  Coaf  i  gToapVNCoaat  t) 

(s  i  at  ti  ah  if  t  ho  toaplato  co  of  i  gar  at  i  oa  to  zoro  t) 
for  i  1  to  HaiToapcol  I  do  bogia 

A.  c  f  g(  .  i  .  ) .  r  :  *  o;  A.cfg(.i.).c  :*  0; 

A.  Cf g(  .  »  .  ).  w  :  '  0, 

oad;  (  l  for  i  do  t ) 
auabor  :  0; 

for  i  -rad  to  -  »  do  bogia 

auabor  :  >  auabor  ♦  I ; 

A.  c  f  g(  .  auabor  .).r  >  i  ;  A.  cf  g(  .  auabor  .  )  .  c  :  o  ; 

a.  c f  g (  .  suabo r  ■  )  •  w  val  uo, 
oad;  (  s  for  i  do  • ) 
for  >  -rad  to  rad  do  bogia 

suabo  r  •  suabo  r  •  i  . 

A.  r  f  g(  .  auabor  .).  r  0,  A.  r  f  g(  .  auabor  .).  c  •  »; 

A.  f  f  g (  .  suabo  r  .  )  .  w  val  uo  , 

ond.  <•  for  i  do  t) 
for>  I  to  rad  do  bo  f  i  a 

suabo  r  •  suabo  r  ♦  i , 

A .  c  f  g  (  .  a  u  abo  r  ,  )  .  r  i,  A.rfg(.au  abo  r  .  )  .  c  0  ; 

A.  r  f  f  (  .  suabo  r  .  )  .  w  val  uo ; 

oad.  (  *  for  i  do  I) 

A.  aua  suabo r  . 


•  ad.  (■  procoduro  r#  af  i  t  To  mm  v  uco  ns  t  «) 


■»% 


ffTWSfrr 


*• 
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procedure  TeapAdd  (  A,  ■  :  Teapl  ate;  var  C  :  Template); 

(a  poiDtMic  addition  of  templates:  C  *  A  ♦  B,  for  A  0  B  «) 

var  cot ;  cot  A»  cot  I  :  integer; 

i iteraect  :  bool eao; 
begin  (>  procedure  TeapAdd  t) 

c  *.  *  a; 

c  ot  :  *  A.  no  a; 

for  cBtl  :  *  I  to  I.  aua  do  begin 
i nteraect  : =  f al ee; 

for  cntA  :  *  1  to  A.  bob  do  (s  add  intersecting  teaplate  pixels  t) 

if  (  A.  cf  g( .  ent  A.  ) .  r  =  B.  cf  g(  .  ent  B.  ) .  r)  and  (  A.  cf  g( .  ent  A.  ) .  c  *  B.  cf  g< .  ent  B.  ) .  c )  and 

(intersect  *  false)  tbeo  begin 
C.  cf  g(  .  ent  A.  ) .  w  :  5  A.  cf  f(  .  ent  A.  ) .  w  B.  cf  g(  .  cot  B.  ) .  w, 
i  nt e r sect  : =  true; 
end;  (t  if  intersect  true  t) 

if  intersect  *  false  tben  begin  (s  append  noni  nt  ersect  *  ng  pixels  s) 


cat  :  -  cat 

C.  ct  t(  .  cat 

♦  i; 

).e  : 

=  *.cf<( 

c  nt  B.  ) .  c ; 

C.  cf  t(  .  cat 

).  r  : 

=  1.  cf*( 

c  nt  B.  ) .  r ; 

C.  cf  t(  .  cat 

)-  w  ; 

--  >.  cf<( 

ent  B.  ) . 

•  ad;  ( t  if 

intersect  false  tben  8) 

•  od;  (»  for 

ent  B 

do  3) 

C.  bub  :  *  ent ; 

end;  (s  procedure  TeapAdd  t) 

procedure  TeapSubtract  (  A.  B  :  Teaplate;  var  C  :  Teaplate); 

(s  point  wise  subtraction  of  teaplates:  C  *  A  -  1,  for  A  U  B  S) 

var  ent ,  ent A, ent B  :  integer; 

i nteraect  :  bool eao; 
begin  (S  procedure  TeapSubtract  s) 

C  :  =  a; 

ent  :  -  A.  bub. 

for  cntl  :  1  to  B.oua  do  be  gin 

i  nt  e  r  sec  t  false; 

for  cat  A  :  I  to  A.  nua  do  (t  subtract  intersecting  teaplate  pixels  S) 

if  (  A.  cfg(  .  ent  A.  ).  r  B.rfgf.cntB.  ).r)  and  (  A.  c  f  g  (  .  c  nt  A.  ) .  c  =  B.  c  f  g(  .  c  nt  B.  ) .  c  )  and 

(intersect  false)  then  begin 
C.cfg(.cntA.  ).*:•  A.  (f  g(  .  cut  A.  )  .  w  -  B .  cfg(  .  cntB.  ).  w, 
intersect  true; 

end;  (S  if  intersect  true  then  i) 

if  intersect  false  then  begin  (S  append  no  nt  nt  e  r  sec  1 1  ng  pixels  8) 
ent  ent  ♦  l ; 

C.  cf  f(  .  rst.  I  .  c  B.tfgl.rniB.  ).c; 

c.  cff(  .  cm.  )  .  r  I.  *  f((  .  rntB.  ).  r; 

C.  cfg(  .CBt.  ).  w  :  -  I .  c  f  g  (  .  cntl.  ).  w, 

end;  (S  if  intersect  false  the  a  t) 
end;  (8  for  cat!  do  t) 

C.  sua  :  ent , 

end.  (8  procedure  TeapSubtract  <) 


.*  / 
.**/ 
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procedure  TeapMul  ti  ply  (  A*  ■  :  Teapl  ate;  var  c  :  Teaplate); 

(a  poiatvaae  aal  1 1  pi  i  cat  i  oa  of  teapl  atea  :  C  =  A  i  I,  for  A  ft  B  t) 

var  cat ,  cat  At  cat  a  :  integer; 
begi  a  (*  procedure  TeapMul  ti  ply  t) 
cat  :  =  0; 

for  cotA  :<  1  to  A.  aua  do 

for  call  :  -  l  to  1.  aua  do  (a  aultiply  >  atertecti  of  teaplate  pi  xel  a  a) 

if  (  A.  cf  K  .  cat  A.  ) .  r-a.  cf  (( .  cat  B.  ) .  r)  aad  (  A.  cf  g( .  cat  A.  ) .  c  =  B.  cf  g<  .  cat  a.  ) .  c)  then  bed  a 

cat  :  =  cat  ♦  l; 

C.  cf  ((  .  cat .  ) .  r  :«  A.  cf  (( .  cat  A.  ) .  r; 

C.  cf  1( .  cat .  ) .  c  :»  A.  cf  g(  .  cat  A.  )  •  c; 

C.  cf  C(  •  cat .  ) .  w  ;s  A.  cf  (( .  cat  A.  ) .  w  a  a.  cf  (t .  cat  a.  )  .  w, 

cad;  (a  if  a) 

C.  aua  :  s  cat; 

end;  (a  procedure  TeapMul  ti  pi  y  a) 

procedure  TeapDi  vi  de  (  A,  8  :  Teaplate;  var  C  :  Teapiate); 

(a  poiatvoce  diviaioa  of  teapl  atea:  C  =  A  /  I,  for  A  n  I  a| 

var  cat  ■  cat  A.  cat  B  :  i  atefer; 
twain  (a  procedure  TeapDi  vi  de  a) 

cat  :  =  o; 

for  cat  A  :  =  1  to  A.  aua  do 

for  enta  :  =  1  to  a.  aua  do  (a  divide  iateraectiaf  teaplate  pi  iel  a  a) 

if  (  A.  cf  (( .  cat  A.  ) .  r  =  a.  cf  g( .  cat  a.  ) .  r )  aad  (  A.  cf  f(  .  cat  A.  ) .  c  *  a.  cf  (( .  cat  a.  ) .  c)  then  beat  a 

cat  :  =  cat  ♦  l; 

C.  cf  ((  .  cat .  ) .  r  A.  cf  (( .  cat  A.  ; .  r; 
c.  cf  (( .  cat.  ) .  c  :•  A.  cf  g(  •  cat  A. ) .  c; 
if  a.  cl  t( .  cat  a.  ) .  w  <>  o 

then  C.  cf  g<  •  cat. ) .  w  :  =  A.  cf  f<  .  cat  A.  ) .  *#  /  B.  cf  (( .  cat  a.  ) .  w 
el  ae  C.  cf  g(  .  c  at .  ) .  w  :  =  o; 
ead;  (a  if  a) 

C.  aua  5  ■*  cat; 

end;  (  a  procedure  TeapDi  vi  de  a) 

procedure  TeapMai  (  A.  a  :  Teupljte;  var  C  :  Teaplate); 

(a  poiilMae  aa  it  aua  of  t  eapi  atea  :  C  =  aal  (A,  a),  for  A  u  a  a) 

var  c  at  i  c at  A,  c ot a  :  integer; 

i  alerted  boolean; 
begin  (a  protedure  TeapMai  a) 

C  A, 

rot  A.  aua; 

fa  for  each  piiel  in  teaplate  a  si 
for  i  it  I  1  to  I  oua  do  begin 

interned  false. 

for  dll  1  to  A.  aua  do  (a  ana/*  interceding  teaplate  piieln  a) 

if  fA.rfgf.rntA.i.c  I.  ifgl  .rail.  I  .i  I  aad  (A.rfgt.cntA.).r  B.rfgf.cotB.  ).r)  aad 
fialersect  falae)  then  begin 

if  A.rfgf  r  at  A.  I .  w  <  l.cfgl  tail.  ).w  t  be  n  c.  r  f  gf  .  t  at  A.  I  .  w  a.  r  f  gf  .  .at  B.  1 .  w, 
i  aleraert  true. 

end.  (a  if  interned  true  then  I) 


A  -  1  0 


if  )  attricct  >  fall*  tbea  bed  a  (1  append  bo  bi  at  •  r  i*c  1 1  ■(  pi  i*l  s  1) 

cat  :  =  cat  ♦  l ; 

C.  cf  (( .  eat .  ) .  r  :  =  8.  cf  (( .  cat  1.  ) .  r; 

C.  cf  (( .  cat . ) .  c  :*  I.  cf  (( .  cat  8.  ) .  c; 

if  1.  cf  ((  .  cat  1.  ) .  w  >  0  tb«a  C.  cf  K  .  cat  •  )  .  w  :*  A.  cf  f(  -  cat  8.  ) .  w 
•  la*  C.  cf  *(  .  cat .  ) .  w  :=  o; 

«ad;  (x  if  i  ateraect  false  tb«a  i) 

•  ad;  (1  for  cata  do  1) 

C.  ana  :  =  cat; 

cad;  (t  procedure  Temp  Hat  t) 

procedure  Tempi*  a  (  A.  1  :  Template;  var  C  :  Template); 

(t  poiatwtie  a  ti  me  of  templates  :  c  =  mb  (A, 8),  for  A  U  8  a) 

ear  cot ,  cat  A,  cat  1  :  i  atefer; 

intersect  :  booleaa; 
tiefi  s  (s  procedure  Tempi*  a  s) 

C  :  =  a; 

cat  :  -  A.  aum; 

(1  for  each  pi  sel  10  template  *  a) 
for  cat!  :  a  1  to  8.  aum  da  b*(l  a 
i  atersect  :  =  f  al  se; 

for  catA  :  3  i  to  A.  aum  do  (t  maxi  mi  le  i  ateraecti  a(  template  pixels  S) 

if  C  A.  cf  ((  .  cat  A.  ) .  c=  8.  cf  (( .  cat  l.  ) .  c)  aad  (  A.  cf  (( .  cat  A. ) .  r  «  >.  cf  f( .  cat  a.  ) .  r )  aad 

(intersect  >  false)  tbea  befl  a 

if  A.  cf*( .  cat  A.  ) .  w  >  l.  cf  (( .  cat  >.  ) .  w  then  C.  cf  ((  .  cat  a.  ) .  w  :  =  8.  cf  f( .  cat  1.  ) .  w; 

I atersect  :  =  true; 

cad;  (S  if  i  atersect  true  tbea  t) 

if  intersect  *  false  tbea  bed  a  (s  append  nooi  at  •  r  sec  1 1  nf  pixels  t) 

cat  ;  :  cat  ♦  J; 

C.  cf((  •  cat .  ) .  r  :  »  8.  cf  ((  .  cat  8.  ) .  r  ; 

C.  cf  ((  .  cat .  ) .  c  :  3  8.  cf  (( .  cat  8.  ) .  c; 

if  I.  cf«( .  cat  8.  ) .  U  <  0  tbea  C.  cf  |( .  cut .  ) .  w  :=  8.  cf*( .  cat  a.  ) .  w 
else  c.  cf  ((  .  c  at .  ) .  w  :*  o; 

end;  (I  if  intersect  false  tbea  t) 
end;  (  t  for  cat  a  do  t) 

C.  aum  t  cat; 

end;  (X  procedure  Tempi*  a  t) 

procedure  TempCi  rrl  eil  ui  (  A.  8  Template;  var  c  Template); 
var  i  ,  c  at ,  c  at  A,  c  at  B  integer, 
dupl  irate  bool  v  a; 

l*|i  B  ft  procedure  *e  mpf'i  r  r  I  e  PI  us  tl 
i  at  o . 

(I  for  each  pi  tel  in  template  8  tl 
for  cats  1  to  8.  sum  do 

ft  compute  8<bt,byl  t  Afat.ay)  for  template  A  centered  tl 
ft  on  Sfbt.by)  aad  turn  the  product  i  at  o  Cf  at  •  bt ,  a  y  •  by  )  tl 
for  catA  ;  to  A.  sum  do  be  da 
duplicate  f  al  te  . 


(s  f  i  ad  t  be  pixel  id  the  template  coo/i  lurati  od  if  it  exists  s) 


for  i  :  =  1  to  CDt  do 

if  (  C.  cf  g( .  i  .  ) .  r  3  a.  cf  g(  .  cat  A.  ) .  r  ♦  B.cfg(.cntB.).r)  and 
(  C.  cf  g(  .  i  .  ) .  c  =  A.  cf  g< .  cot  a.  ) .  c  ♦  B.  cf  g(  .  cot  B.  ) .  c) 
then  begi  a  (s  pixel  already  is  template  C  configuration  s) 

C.  Cf  <(.».)•  w  :s  C.cf((.t.)>w  ♦  A.  Cf  ((  .  COt  A.  )  .  w  t  B.  cf  g  (  .  c  it  B.  )  .  w, 

dupl  i  c at e  :  *  true; 

end;  (t  if  pixel  found  t) 

if  duplicate  =  false  then  begin  (s  append  aoai  at er tec 1 1  ag  pixels  s) 
cat  :  *  cat  ♦  1; 

C.  cf  g( .  c  at .  )  •  r  :=  A.  cf  g( .  cot  A. ) .  r  ♦  B.  cf  g( .  cat  B.  ) .  r; 

C.  cf  g( .  cot .  )  •  c  :=  A.  cf  g( .  cnt  A.  )  •  c  ♦  B.  cf  g( .  cnt  B.  ) .  c; 

C.  cfg(.cnt.).w:=  A.cfg(.  cntA.  ).  w  t  B.cfg(.catB.  ) .  w; 
end;  (t  if  duplicate  false  s) 
end;  (s  for  cnt  A  do  t) 

C.  nua  :  -  cnt; 

end;  (i  procedure  TeapCi  rcl  e PI  us  t) 


procedure  TeapCi  rcl  e  Max  (  A,  B  :  Template;  vir  C  :  Template); 
var  i  »  cat  i  cnt  A*  cat  I  :  integer; 

duplicate  :  boolean; 
begin  (t  procedure  TeapCi  rcl  eMax  s) 

cat  : =  0; 

(t  for  each  pixel  in  template  B  s) 
for  catB  :=  l  to  B.  nua  do 

(s  compute  B<bx,by)  s  A(ax,ay)  for  template  A  centered  s) 

(s  OB  B(bx.by)  and  maxi  m  xe  t  be  product  in  C(  a  t  ♦  bi ,  a  y  ♦  b  y )  s) 
for  citA  :  s  l  to  A.  sum  do  begi  a 
dupl  i  c  at  e  :  «  false 

(s  find  the  pixel  in  the  template  c  o  of  i  g  ur  at  i  o  a  if  it  exists  S) 
'or  t  l  to  cat  do 

If  (  C.  cf  g(  .  i  .  ) .  r  *  a.  cf  g(  .  cat  a.  ) .  r  ♦  B.  c  f  g<  .  c  at  B.  ) .  r  )  and 

(C.rfg(.i.).e  A.  cfg<  .  c  at  A.  ) .  c  ♦  ».  c  f  g(  .  c  at  B.  )  .  c  ) 

then  begin  (t  pixel  already  ia  template  C  coif  i  g  ur  at  i  on  i) 
if  A.  cf  g(  .  cat  A.  > .  w  t  B.  c  f  g(  .  e  at  I.  ) .  w  >  C.cff(  .  i  .  ).  w 

then  C.  cf  f  (  .  i  .  ) .  w  A.  rffl  .  cat  A.  ).w  i  B.cfgf.tatB. 
dupl  i  r  at  e  t  r  ue  , 

•  ad.  (t  if  pixel  found  •) 


if  duplicate  fa! 
cat  cat  •  i . 

f'.cfgf  rat.  ).  r 

r.ilfl  rot  1  < 

•  f  ( i  -at 

Mid.  I  if  dupl  >  i 
end.  <1  for  nt  4 

nua  -  nt  , 


■  e  then  be  g  ■  a 

A.  t  f ((  .  c  at  A. 
*.  '  f|(  i  nl  A 
4.  •  f  C  t  -  n  t  4 
ate  fain#  tv 
do  t  I 


ft  appeed  nom  nirrin  t*  i|  pi  i  I  • 

r  *  t.  •  f (i  .  •  st  I.  I  .  r 

»  •  I  i  ((<  si  I  i  i  . 

to  I  R.  .  f  (I  <ntll  I  to 


a  « > 


end.  (I  procedure  Te  mp <  i  rrl  #  H* 


A 


proc* dur*  TfBBCi  rcl  *IB  a  (  A.  ■  :  T*^>;  at*;  »ar  c  :  racial*); 

»ar  i  ,  cat ,  cat  A,  cat  ■  :  latagar; 

dap)  i  cat*  :  boot  via; 
b*t>  a  (a  procvdur*  Ttapci  rcl  •  **  a  () 

c at  :  >  o; 

(l  for  *acb  pi  i*l  i*  t*aplat*  •  t) 
for  cat  *  :  *  I  to  I.  bub  do 

(a  coapgit  Kbi.br)  »  A(aa.ay)  for  t*aplat*  A  c*at*r*d  a) 
(a  oa  B(  ba.br)  aad  ■  n  a  11  t  b*  prodact  i  a  C(  at  ♦  bi ,  a  y*  br )  *) 
for  catA  :  •  l  to  A.  aoa  da  b»*i  a 
dapl  i  c  at  »  :  «  f  al  a*; 


(  t  ft  ad 

t  ba 

P» 
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■  a 

t  ha  t 

•  opl  at  a  c  oaf 
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.  C 
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w  <  C. 
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du  pi 
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t  r 

ua  . 

•  nil. 
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cat 
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<  1  g< 

at  B.  i 

C  . 
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'  * 

a  I  • 
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for 

n 
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do 

* 
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>  r  •  *«  n  « 
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if  duplicate  =  false  then  begi  a  (s  appead  no  ai  at  e  r  aec  1 i  og  pixels  s) 


cat  :  = 

cat 

♦  i ; 

C.  cf  g(  . 

ent 

) .  r 

--  a.  cf  *( 

c  at  A.  )  -  r  ♦ 

*.  cf  if 

c  nt  B.  ) .  r ; 

C. cf g( . 

cat 

) .  c 

-  A.  Cf  if 

c  nt  A.  ) .  c  ♦ 

».  cf  if 

cat  B.  ) .  c; 

C.  cf  g<  . 

c  nt 

) .  w 

:  =  A.  cf  if 

cot  A.  }.  W  ♦ 

».  cf  tf 

c  at  B.  ) .  w, 

end;  (  x 

i  f 

dupl 

icate  false  s) 

od;  (  x 

for 

ent 

A  do  >) 

C •  bum  :  *  cat; 

end.  (t  procedure  Te apsquare Mai  >} 

procedure  TeswSquarettn  (  A,  *  :  Teaplate;  v*r  C  :  Teaplate); 
var  i  ,  c at ,  c at  A. c at  1  :  integer; 

duplicate  boolean; 
begi  a  (  t  procedure  TeapSquaret*  n  *> 
cat  o , 

ft  for  *ir  h  o»  *  e  I  in  t  e  nol  at*  B  t ) 
for  cotl  :  1  to  B.  nua  do 

(t  coapotr  B(bSiby)  ♦  A(  ai.  ay)  for  template  A  centered  S) 

(t  on  B(bx,by)  and  a  di  ■  zt  the  sua  in  C(  ax  ♦  bx  *  a  y*  by)  S) 
for  cat  A  :  l  t  o  A.  oua  do  begi  n 

dupl  »  c  at  e  :  =  false’ 

(s  find  the  pi  iet  in  the  teaplate  co  of »  g  ur  at  i  o  t  if  it  exists  t) 

for  i  =  l  to  cot  do 

if  (  c.  c  f  g(  .  i  .  ) .  r  =  A.  cf  g(  .  cat  A.  ) ,  r  ♦  i.  c  f  g(  .  c  nt  B.  ) .  r )  and 

(  C.  cf  g(  .  i  .  ) .  c  -  A.  cf  g(  .  ent  A.  )  •  c  ♦  B.  cf  g(  .  ent  B.  ) .  c) 

then  begin  (t  pixel  already  in  teaplate  C  configuration  X) 
if  A.  cf  g(  .  cat  A.  ) .  w  ♦  B.  cf  g(  .  cot  B.  )  •  *  <  C.cfg(.i.).w 
then  C.cfg(.i.).w  :=  A.  c  f  g(  .  c  nt  V  ) .  v  ♦  B.  cf  g( .  ent  B.  ) .  w; 
dupl  \  cate  :  -  true; 
end;  (X  if  pixel  fc-nd  x) 
if  duplicate  -•  false  then  begin 
ent  :  -  ent  *  1 ; 

C.  cf  (( .  ent .  ) .  r  a.  cf  g(  .  c  nt  a.  ) .  r  ♦  B.  cf  g(  •  ent  B.  ) .  r  ; 

C.  cf  gC  .  c  nt .  ) .  c  A.  c  f  g(  .  c  nt  A.  ) .  c  ♦  B.  c  f  g(  .  c  nt  B.  ) .  c ; 

C.  cf  g(  .  ent .  ) .  w  a.  cf  g(  .  ent  A.  ) .  w  ♦  B.  cf  g(  .  ent  B.  ) .  w, 

end;  (X  if  duplicate  false  X) 
end;  (X  for  cli  A  do  X) 

C.nua  :  r  ent; 

end;  (X  procedure  TeapSquareMi  n  x) 

procedure  TeapScalarMultiply  (A  :  Teaplate;  value  Te  apl  a  t  e  T  »  pe  ;  var 
(>  aulti  plication  of  an  «aa(e  by  a  scalar  c  A  x  value  t ) 
var  ent  integer; 

begin  (*  procedure  TeapStalarMultiply  t  ) 

c  .  a; 


for  ent  1  to  A.nuado  C.  r  f  f(  .  cot  .  )  .  v 

end,  (X  procedure  Teapot  a  I  arHul  t  «  pi  >  t) 
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(*  BINARY  IMAGE  -  TEMPLATE  OPERATIONS  t) 

( iiiimitiimtiititiiiiiismiittiitmtumtiMumtuiiumtitmiiit) 

procedure  I  aageTeapCi  rcl  ePl  us  (  A  :  lu(e;  ■  :  Tenet  ate;  var  C  :  laage); 

(*  convolve  A  wi  t  h  template  ■  t) 
var  Row,  Col  ,  i  ,  j,  ent  :  i  nte(er; 

Sub  :  IaageType; 

begin  (*  procedure  I  aafeTeapCi  rcl  ePl  us  *) 

C.  row  ;s  A.  row;  c.  col  :=  A.  col; 
for  Row  :  =  l  to  A.  row  do 
for  Col  :  =  l  to  A.  col  do  begi  n 

(*  sua  the  products  of  i  aa(e  pixels  and  overl  ayed  template  pixels  s) 

Sum  :  =  o; 

for  ent  :  =  I  to  B.  odb  do  begi  n 
i  :  =  B.  cf  (( ■  ent . )  -  r;  j  B.  cf  t(  •  ent .  ) .  c; 

if  ( (  Rowoi  >  =  1 )  and  (  Rowoi  <  =  A.  row)  and  (Col*j>  =  l)  and  (  Col  *i  <  =  A.  col  )  ) 
then  Sua  :  =  Sub  ♦  A.  gray! .  Rowoi  ,  Col  *j .  )  x  B.  cf  E( .  ent . ) .  w; 
end;  (S  for  ent  do  l) 

(*  set  the  new  pixel  value  equal  to  the  sub  x) 

C.  (ray( .  Row,  Col  .  )  :=  sub; 
end;  (t  for  Col  do  t) 
end;  (s  procedure  I  aafeTeapCi  rcl  ePl  us  X) 

procedure  I  aageTeapCi  rcl  eMax  (  A  :  Iaafe;  B  :  Template;  var  C  :  I  mate); 
var  Row,  Col  ,  i  ,  j,  cat  :  integer; 

Maxval  :  IaageType; 

begin  (t  procedure  I  aageTeapCi  rcl  eMax  *) 
c.  row  :  =  A.  row;  c.  col  :  *  A.  col  ; 

Const  I  mage  (  C,  A.  row,  A.  col  ,  0) ; 

(*  for  eacb  pixel  position  s) 
for  Row  :  =  l  to  A.  row  do 
for  Col  :  =  l  to  A.  col  do  begin 

(t  set  aaxval  to  the  product  of  present  i  aage  pixel  and  teapl  ate  pixel  center  s) 

ent  :  =  l; 

while  (  B.  cfg(.  ent.  ).  r  <>  o)  and  C  B.  cfg( .  ent . ) .  c  <>  0)  do  ent  :=  ent  ♦  l; 
if  (  B.  cf  g( .  ent .  ) .  r  -  o)  and  (  B.  cf  g( .  ent .  ) .  c  *  0) 
then  Maxval  :=  A.  gray< .  Row,  Col  .  )  x  B.  cf  g( .  ent .  ) .  w  (*  template  center  defined  *) 
else  Maxval  =  o;  (s  teapl  ate  center  undefined  *) 

(s  search  for  the  local  taaxi  aua  product  s) 
for  ent  :  =  i  to  B.  ana  do  begi  n 
i  :=  B.  cf  g(  .  ent.  )  .  r;  j  :=  B.  cf  g( .  ent.  )  .  c; 

if  ( (  Rowti  >  :  i)  and  (RowticA,  row)  and  (Coltj>  =  l)  and  (  Col  ♦  j  <  =  A.  col  ) )  then 
if  (  A.  gray( .  Rowti  ,  col  ♦  j  .  )  x  B.  cf  g( .  ent,  ) .  w)  >  Maxval 
then  Maxval  :=  A.  gray( .  Rowoi  ,  Col  .  )  t  B.  cf  g(  .  ent.  )  .  w, 
end;  (t  for  ent  do  a) 

(t  set  the  new  pixel  value  equal  to  the  aaxi  bub  s) 

C.  gray( .  Row,  col  .  )  :  =  Maxval; 
end;  (x  for  Col  do  x) 
end;  (X  procedure  I  aageTeapCi  rcl  eMax  s) 
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proceduie  I  aageTeapCi  rcl  eM  a  (  A  :  {mage;  (  :  Tempi  at e;  var  C  :  laage); 
var  Row,  col  >  i  ■  i,  cat  :  i  nteger; 

M  aval  :  I  mage  Type; 

begi  n  (*  procedure  I  aageTeapCi  rcl  till  *) 

C.  row  :  =  A.  row;  c.  col  :*  A.  col; 

Coast  Inage  (  C,  A.  row,  A.  col  ,  0) ; 

(t  for  each  pixel  poaitloa  <) 
for  Row  ;  *  1  to  A.  row  do 
for  Col  :>  l  to  A.  col  do  btfla 

(*  set  ad  aval  to  the  product  of  preseat  iatft  pixel  aad  template  center  pixel  x) 

cat  :  *  1; 

while  (  a.  cf  g(  .  cat .  ) .  r  <>  0)  aad  (  B.  cf  g(  .  cat .  ) .  c  <>  o)  do  cat  :  =  cat  *  l; 

if  (  B.  cf  g(  .  cat .  ).  r  =  0)  aad  (  B.  cf  g( .  cat .  )  .  c  =  0) 

then  Mi  aval  :  =  A.  gray( .  Row,  col  . )  X  B.  cf  g( .  cat.  ) .  w  (X  template  center  defined  X) 
else  Mi  aval  ;  =  o;  (x  teaplate  center  undefined  x) 

(X  search  for  the  local  ai  ni  aua  product  x) 
for  cat  l  to  I.  ana  do  begi  a 
i  :  -  B.  cf  g( .  cut .  )  .  r; 
j  :=  B.  cf  g( .  cat.  ).  c; 

if  ( (  Rowei  >  =  l)  and  (  Rowei  <  =  A.  row)  and  (Col*j>  =  l)  aad  (  Col  ♦  i  <  *  A.  col  )  )  then 

if  (  A.  gray( .  Rowei  ,  Col  ej. )  X  B.  cf  g(  .  cat.  )  .  w)  <  M  aval 

then  M  aval  :*  A.  gray( .  Rowei  ,  Col  ej.  )  X  B.  cf  g(  .  ent .  ) .  w; 
end;  (x  for  cat  do  X) 

(X  set  the  aew  pixel  value  equal  to  the  m  ni  aua  x) 

C.  gray( .  Row,  Col  .  )  Mi  aval  ; 
end;  (X  for  Col  do  X) 
end;  (x  procedure  I aageTeapCi  rcl  eda  x) 

procedure  I  aageTeapSquareMax  (  A  :  laafe;  B  :  Teaplate;  var  C  ;  laage); 
var  Row,  col  ■  i  ,  j,  cat  :  integer; 

Maxval  :  Inage  Type; 

begin  (X  procedure  KaageTeapSquareMax  X) 

C.  row  :  =  A.  row;  c.  col  :  »  A.  col  ; 

Coastlaage  (  c,  A.  row,  A.  col  ,  0) ; 

(X  for  each  pixel  position  x) 
for  Row  ;s  1  to  A.  row  do 
for  Col  :  =  1  to  A.  col  do  begin 

(x  set  aaxval  to  the  sua  of  present  image  pixel  and  teaplate  pixel  center  x) 

ent  :  =  t; 

while  (  B.  cf  g( .  cat .  ) .  r  <>  o)  and  (  B.  cf  g( .  cat .  )  .  c  <>  o)  do  cat  :=  cat  t  l; 

if  (  B.  cf  g( .  cat.  ) .  r  =  o)  and  (  B.  cf  g( .  cat.  ).  e  =  0) 

then  MaxVal  :=  A.  gray( .  Row,  Col  .  )  ♦  B.  c  f  g( .  c  at .  ) .  w  (X  teaplate  center  defined  x) 

else  Maxval  :=  A.  gray( .  Row,  col  .  ) ;  (X  teaplate  center  undefia'd  x) 

(X  search  for  t  be  local  naxi  aua  sua  X) 
for  cot  :  =  1  to  B.  aua  do  begi  a 
i  B.  cf  r( .  cat. ).  r;  )  ;=  b.  cf  g( .  cat . ) .  c; 

if  ((Row*i>  =  i)  aad  (  Rowei  (  =  t,  row)  aad  (Col»j>  =  2)  and  (  Col  ♦  j  <  »  A.  col  ) )  then 
if  (  A.  gray( .  Rowei  ,  Col  . )  ♦  B.  cf  g(  .  cat.  )  .  w)  >  MaxVal 
then  Maxval  :=  A.  gray( .  Rowei  ,  Col  ♦  i .  )  *  B.  cf  g( .  cut.  ) .  w; 
end;  (X  for  cat  do  X) 
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(x  set  tbe  mw  pixel  value  equal  to  t  be  aaxi  mum  t) 

C.  gray(  .  Row,  col . )  :»  Haxval  ; 
end;  (s  for  Col  do  *) 
eod;  (t  procedure  I  mage Temp Square  Max  *) 

procedure  I  aageTempSquareM  a  (  A  :  I  stage;  I  :  Template;  war  C  :  Image); 
war  Row,  col  ,  i  >  j,  cat  :  integer; 

F*  aval  :  I  mage  Type; 

begin  (x  procedure  I  mageTempSquareM  n  x) 

C.  row  ;s  A.  row,  C.  col  :=  A.  col  ; 

Const  I  mage  (  C,  A.  row,  A.  col  ,  0) ; 

(x  for  eacb  pixel  position  X) 
for  Row::  l  to  A.  row  do 
for  Col  :  =  1  to  A.  col  do  begin 

(x  set  as  oval  to  tbe  sub  of  present  i  stage  pixel  and  template  center  pixel  x) 
ent  :  =  l; 

while  (  B.  cf  g(  .  cat .  )  •  r  <>  0)  and  (  B.  cf  g( .  ent .  )  .  c  <>  0)  do  cat  :*  ent  *  l; 

if  (  B.  cf  g( .  ent. ) .  r  =  0)  and  {  B.  cf  g( .  ent .  )  .  c  =  0) 

then  rt  nVal  :=  A.  gr  ay(  .  Row,  col  .  )  ♦  B.  cf  g(  .  ent .  ) .  w  (x  template  center  defined  x) 

else  M  aval  :=  A.  gray( .  Row,  col  .  )  ;  (x  template  center  undefined  a) 

(x  search  for  the  local  minimum  sum  x) 
for  ent  :  =  l  to  B.  num  do  begin 
i  :  -  B.  cf  g( .  ent. ) .  r; 
j  :  =  B.  cf  g( .  ent . ) .  c; 

if  ((Rowti>  =  l)  and  (  Rowsi  CA.  row)  and  (Col+j>  =  l)  and  (  Col  ♦  j  <  =  A.  col  )  )  then 
if  (  A.  gray(  .  Rows!  ,  col  ♦  j .  }  ♦  B.  cf  g( .  ent .  )  .  w)  <  M  aval 
then  M  nval  :  =  A.  gray( .  Rowsi  ,  Col  ♦  j.  )  ♦  B.  cf  g( .  cat .  )  .  w, 
end;  (X  for  ent  do  X) 

(X  set  the  new  pixel  value  equal  to  the  minimum  x) 

C.  gray( .  Row,  col  .  )  :  =  Ml  aval  ; 
end;  (x  for  col  do  x) 
end;  (X  procedure  I  mageTempSquareM  n  X) 
end.  (X  Module  Real  _l  A.Operatl  ons  X) 
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Appendix  B:  AFIT  Image  Algebra  Preprocessor 


program  I  mageAl  gebraPreprocessor  ( i  aput ,  out  put ) ; 

II  Tbit  program  reads  a  Mia  containing  t  ha  Image  AJ  gebra  description  of  a) 
((  aa  I  mage  processing  algorithm,  translates  it  to  VAX  PASCAL  source  coda  s) 
(a  aad  procedure  cal  I  s.  and  outputs  the  results  to  aa  external  file  for  a) 
(a  subsequent  compilation.  This  program  is  written  in  VAX  PASCAL.  a) 

(a  Capt  Christopher  J  Titos,  Air  Force  Institute  of  Technology,  GE-86D  a) 

(a  global  constants,  types,  and  variables  a) 


const  MaxNode  = 

50; 

c* 

Maxi  mum 

number  of 

nodes  i 

D 

t  he  bi  nary  tree 

s) 

Buf  fer  Leo  = 

136; 

(* 

•axi  MUM 

character 

1  ength 

of 

I  A  command  1  i  ne 

t) 

NameLen 

so; 

(* 

Maxi  MUM 

character 

1 engt  h 

Of 

vari  abl  e  names 

*) 

MaxNames 

100; 

(* 

Maxi  MUM 

number  of 

vari  abl  e 

names  of  one  type 

*) 

Default  Image  =  ’real';  (a  default  data  type  of  i  mage  gray  levels  a) 

Defaul  t  Tempi  ate  =  ’real';  (a  default  data  type  of  template  weights  a) 

type  BufferStr  =  varyi  ng[  Buffer  Lea]  of  char; 

NameStr  =  var  yi  ng(  Name  Lea]  of  char; 

Intermedoper  =  record  exist  :  boolean;  end; 

Tree  =  array  ( .  l.  .  Max  Node.  )  of  record  (a  binary  tree  implemented  with 
value  ;  NameStr;  (a  a  double  linked  list  structure 

parent  :  integer;  (a  using  an  array  for  data  storage 

I  chi  I  d  :  i  nteger; 
rchi  1  d  :  I  nteger; 
end;  (a  tree  record  a) 

Variable  List  «  array  ( .  l. .  MaxNames. )  of  varyi  ng[  NameLen]  of  char; 

var  i  HI,  i  R2,  i  US,  i  M,  i  RS,  111,  tR2,  tR3,  tM,  tBS  :  Intermedoper; 

Name,  I  ndent ,  Temp  Name  :  NameStr; 

buf  fer,  oi  dbuf fer  :  BufferStr; 

I  mage  Li  st ,  Tempi  ateLi  at ,  seal  ar  Li  st  :  Vari  abl  e Li  st ; 

Numl  mage  Names,  NumTempI  ateNames,  NumScal  arNames  :  integer; 

I  nFi  I  e,  TmpFi  I  e,  Pf  i  I  e  :  text; 

Expression  :  Tree; 

I  ,  j  ,  nodept  r,  nodeent ,  code  :  integer; 

I  mageType,  Tempi  atetype  :  varyi  ng[7]  of  char; 

heel  aredType  :  varyi  ng(«]  of  char; 


procedure  Aaai  f  nl  at er me dResul  t  (operltype  :  char;  operator  :  Namestr; 

var  Beful  t  :  Namestr); 

(a  assign  as  intermediate  reaul  t  for  t  be  operatioa  a) 
befi  a  (a  procedure  Aaai  fal  ntermedResul  t  a) 

Re  a  til  t  :  =  ’  • ; 

if  (operator  <>  ■:=’)  aad  (operator  <>  tbea  befin 

(a  aaai f n  aa  I  stormed!  ate  result  for  the  operatioa  a) 
if  operltype  =  *  l*  tbea  befi  a 

(a  operatioa  is  i  mife-imafCi  i  mafe-templ  ate,  i  made-seal  ar,  01  oaary  a) 
(a  i  mafe;  the  result  is  aa  i  mate;  aaai  (a  intermediate  i  mate  operands  a) 
i f  i  Rl.  eai at  =  f al  se 

tbea  befi  a  Result  :  =  MRi*;  i  Rl.  exist  :  =  true;  end 
else  if  iR2. exist  =  false 

then  befi  n  Result  :  =  'it!';  i  R2.  exist  true;  end 
else  if  i  R3.  exist  -  false 

tbea  bef  i  n  Result  :  =  *  i  R  3  *  ;  i  R3.  exist  :  =  true;  ead 
else  if  i  R4.  exist  =  false 

tbea  befi  a  Result  :  =  1  i  R4’ ;  i  Rt.  exist  :=  true;  ead 
else  if  IRS. exist  =  false 

tben  beRk  a  Result  :  =  'IRS';  i  RS.  exist  :  =  true;  ead 
el  se  befi  a 
Resul  t  :  -  •  ?'  ; 
wri  tel  a; 

wri  tel  a  ('ERROR  ao  more  intermediate  i  mates  available  ERROR'); 
wri  tel  a  ('error  expression  caa  not  be  evaluated  ERROR'); 

ead; 

ead;  (X  if  i  mate  operatioa  X) 
if  operltype  =  *t'  tbea  befi  a 

(a  operation  is  template-template  or  template-scalar;  t  be  result  s) 

(a  is  a  template;  rel  eaae  intermediate  template  operands  a) 

I  f  t  Rl.  exi  at  =  f  al  se 

then  befi  n  Result  :  =  'tRl';  tRl.  exist  :*  true;  end 
else  if  tR2.  exist  =  false 

thou  befi  u  Result  :  =  *tR2’;  tR2. exist  :  =  true;  end 
else  if  tR3.  exist  =  false 

tben  befi  n  Result  :  s  •  tR3';  tR3.  exist  :  -  true;  end 
else  if  tR4.exist  *  false 

tben  beti  a  Result  'tR4';  IR4.  exist  :=  true;  end 
else  if  IRS.  exist  *  true 

tben  befi  n  Result  'tRS';  tRS. exist  :=  true;  end 
el  se  befi  u 
Resul  t  :  =  '  ; 

wri  tel  n; 

wri  tel  n  ('ERROR  so  more  intermediate  templates  available  ERROR'); 
wri  tel  u  (’ERROR  expression  cau  not  be  evaluated  ERROR'); 

end; 

cud;  (a  if  template  operation  a) 
ead;  (a  if  operator  tbea  a) 
end;  (a  procedure  Assi  fal  atermedResul  t  a) 


procedure  *iii(gleiolti  ( operand! .  operands,  operator  :  NaaeStr; 
var  Pea ul  t  :  NaaeSt  r; 
var  expreaa:  Tree;  aodeptr:  i  ate(er) ; 

(a  update  *  be  ezpreiti  ob  tree  to  reflect  the  reaul  t  a) 

(*  of  t  be  operation;  rel  eaae  any  i  ateraedi  ate  reaul  ta  t) 
begi  a  (<  procedure  Aaai  gnResul  ta  t) 
if  operator  »  •  tbea  Reaul  t  :  >  operand2; 
if  operator  -  tbea  Reaul  t  :»  'dotval'; 

((  update  t  be  expreaaioa  tree  to  reflect  t  be  operation  reaul  t  a) 
expreaa( .  no  dept  r. ) .  val  ue  :=  Reault; 


(* 

rel ease  unused 

oter 

nodi  ate  operators 

a) 

i  r 

((operandl^ 

• »  m 

> 

or 

(  oper  and2  - 1 

i  Rl* )) 

and 

(  Resul  t  <  >*  i  Rl*  ) 

t  ben 

i  Rl .  exi  at 

i  r 

((operandl: 

’  i  R2 

> 

or 

( oper and2  =  ' 

i  R2* )  ) 

and 

( Retul t  <  >  * i  R2* ) 

t  ben 

i  R2 .  exi  at 

•  r 

( ( operandl: 

*  i  R3 

> 

o  r 

( oper  aod2  =  ’ 

i  R3*  ) ) 

and 

(  Reaul  t  <  >  *  i  R3*  ) 

t  hen 

i  R3.  exi  at 

i  r 

( ( operandl: 

•  i  R4 

> 

or 

(operands:1 

ill'll 

and 

( Resul  t  <>* i  R4»  ) 

t  hen 

i  R4 .  exi  at 

i  r 

( ( operandl : 

•  »  R5 

> 

or 

( operand2= * 

i  RS* )  ) 

and 

( Resul t  <>  *  1 RS* ) 

t  ben 

i  RS.  exi St 

i  f 

( ( operandl: 

*  t  Rl 

) 

or 

( oper and2  =  1 

t  R 1 '  )  ) 

and 

( Resul t  <>*  t  Rl* ) 

t  ben 

t  R 1 .  exi  st 

i  r 

( ( operandl : 

•  t  R2 

> 

or 

( oper and2 = » 

tR2*  )) 

and 

(  Resul  to*  tR2*  ) 

t  hen 

t  R2 . exi at 

i  t 

( ( operandl : 

•  t  R3 

> 

or 

(operands:* 

t  R31 )  ) 

and 

(  Resul  t  <>*  t  R3*  ) 

t  hen 

t  R3.  exi  at 

i  f 

( ( operandl: 

'  t  R4 

) 

or 

( operand?:* 

tR4* )) 

and 

(  Resul  t  o*  t  R4*  ) 

t  hen 

t  R4.  exi St 

»  r 

((operandl: 

»  t  R3 

> 

or 

( operand2=* 

IRS'  )) 

and 

( Resul t  <>*  t RS* ) 

t  ben 

t  RS.  exi  st 

i  f 

(  operandl : 1 

?*  )  or 

(  operand2:'  T* ) 

t  ben 

expreaa! .  BO  dept  r.  ).  val  ue  : 

=  >  ?•  ; 

express! .  aodeptr.  ) .  I  chi  I  d  :=  0; 
expreaa( .  aodeptr.  ).  rchi  I  d  :  =  o; 
end;  (a  procedure  Aati  gnResul  t  a  a) 

procedure  Bi  aar  yOper  at  i  on  (ear  expreaa  ;  Tree;  aodeptr  ;  integer); 

(a  write  the  binary  operation  procedure  calla  to  t  be  output  file  a) 
ear  operator,  operaadt.  operands,  Reaul  t,  teapoper  :  NaaeStr; 

operltype,  oper2type,  teaptype  :  ebar; 
begin  (a  procedure  Bi  aar  yOper  at  I  on  a) 

(a  retrieee  t  be  operauda  froa  t  be  binary  tree  a) 
operandl  :=  expreaa( .  expreaa( .  aodeptr.  ).  I  cbi  1  d.  ).  val  ue; 
operator  ;  =  expreae(  .  node pt  r .  ) .  eal  ue; 

operand!  exprese( .  expresa( .  node  pt  r .  ) .  rc  bi  I  d.  ) .  val  ue; 

(a  drum  nr  t  be  operand  types  a) 
operltype  ’  >;  operXtype  :  =  •  * ; 
for  i  u  Mo  Nual  aageNaaea  do  begin 
if  operandl  :  1  mage  Li  at  ( .  i  .  )  then  operltype  :  =  •  i  • ; 
if  operand2  =  I  aage  li  at  ( .  i  .  )  then  oper2type  :  = 
end;  (a  for  i  do  x) 

for  i  1  to  NuaTeapI  ateNaaes  do  begi  n 
if  operandl  -  Teapl  at  e  Li  at  (  .  i  .  )  tbea  operltype  »  t  •  ; 

if  operand2  =  Teapl  ate  Li  at  (.  i  .  )  tbea  oper2type  :=  » t  •  ; 

end;  (a  for  i  do  a) 
for  i  1  to  NuaScal  arxaara  do  begi  n 
if  operandl  -  Seal  ar  Li  at  ( .  i  .  )  tbea  operltype  :  =  ’s'; 

if  operands  -  Seal  ar  Li  at  ( .  i  .  )  tbea  oper2typ  :=  'a*; 

end;  (a  for  i  do  a) 
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if  operltype  >  ’  1  the*  begin 

op«r  it  ype  :  *  ’  : 

for  i  :  *  1  to  I  engt  h(  oper  andl )  do 
if  (operltype  *  ’s')  and  (operandl(i]  in  [,-,i'0,..'9t  J) 
t  bon  oporlt  ypo  z  -  'o' 
ol  >«  oporltypo  :  *  •  '  ; 
oad;  (t  if  oporltypo  tbca  a) 
if  oporltypo  >  •  *  theo  bo(i  a 
oporlt  ypo  :  >  ’  o'  i 
for  i  :  *  l  to  I  eagt  h(  o  per  a  add)  do 
if  (opentype  *  ’*’)  aad  (operaad2[i  ]  ia  C  • *  0* . .  ’  9’ ] ) 
t  boa  oper2t  ypo  :  =  *  s* 
el  «#  oper2t  ypo  : =  ’  • ; 
oad;  (2  if  oper2type  then  a) 
i  f  operlt  ype  =  1  *  t  boa 

wri  tel  a  ( 1  ERROR  ’, operandl,’  ia  not  aa  image,  template  or  (cal  ar  oporaad  ERROR*  ) 
if  oper2t  ypo  =  *  ’  t  boa 

wri  tol  a  (’ERROR  ’,operand2,’  is  aot  aa  i  aafti  template  or  scalar  operand  ERROR’) 

if  (operandl  <>  ’2’)  aad  (operand2  <>  ’  ?' )  thea  begin 
write  (  Pf  i  I  e,  I  adeat) ; 

(S  arrange  t  be  operands  ia  the  ;  roper  order:  i-a,i-t,t-s  *) 

if  ((operlt ype  =  ’s')  aad  (opentype  =  *i’))  or 

((oporltypo  =  •  t  ’  )  aad  (oper2type  *  *i*))  or 

((opentype  =  ’s’)  and  (oper2type  =  ’t’))  then  begi  a 

teapoper  :  =  operaadl; 
teaptype  :=  operltype; 
operaadl  :  >  operaad2) 
operltype  :  =  opentype; 
oper and2  :  =  teapoper; 
opentype  :=  teaptype; 
ead;  (S  arraafe  operands  a) 

As  si  cal  ateraedResul  t  (  oper  It  ype,  operator,  Re  sul  t ) ; 


if  (operltype  =  ’  t  ’  )  and  (oper2type  =  ' 
(s  i  aafe-i  aafe  binary  operations  a) 
if  operator(l)  ia  ’a’, 


)  then  begi  n 


’  >’  ] 


then  case  operator(U  of 

i  ,’  wr  iteln  (Pfile,’l  aafe  Add  ( ' ,  o  per  a  n  d  1 ,  ’  ,  '  »  o  pe  r  a  n  d2 ,  *  ,  '  ,  Res  ul  t ,’);’) , 

’-’  wri  tel  a  (  Pfi  I  e,  ’  I  aagesubtract  (’,  oper  andl oper  and2, ’,’>  Resul  t  ,');•) ; 


’a’  :  if  ( 1  engt  h(  operator)  >  I)  and  (operator[2)  =  ’a') 

then  writeln  (  Pfi  I  e,  *  I  aafeExponent  ( ’  ,  oper  andl ,  ’  ,  ’  ,  oper  a  nd2 ,  ’  ,  '  ,  Re  sul  t , 


else  writeln  (  Pf  i  1  e ,  ’  I  inage mil  t  i  pi  y  (’.operandl, 


o  pe  ran  d2,’,’, Res  ult, ’);’); 

' /*  writeln  (  Pf  i  1  e,  '  I  mage  Di  vi  de  (’,  operandl  ,*,*,  oper  and2,  ’,’,  Resul  t  ,*);’) ; 

•  »’  writeln  (  Pf  i  I  e,  ’  l  aageMai  (’,  operandl, operand2, Resul  t ,’);') ; 

:  writeln  (  Pf  i  I  e,  ■  1  aifeM  a  (’,  operandl,  ’,’,  operand2.  ’,*,  Resul  t ,’);’) ; 

:  writeln  (  Pf  i  1  e,  '  I  aagecot  (’,  oper  andl oper  and2,  ’,  do  t  ral  );•) ; 


writeln  (  Pf  i  I  e,  o  pe  r  and2 ,  ’  :=  *,  o  per  andl ,’;’) ; 
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otherwise  writelo  (  PFi  I  e>  1  Chari  safe  (’, operaodl, 

',’’’, operator opera od2,  ’,’,  Resul  t ,’  ) 

end  (<  case  operator  s) 
el  se  begi  n 
wri  tel  d; 

writelo  ('ERROR  operator  '■  operator!  '  does  not  exist  ERROR’); 

writelo  (’ERROR  for  i  u(e-i  aa(e  operations  ERROR’); 

operaodl  :=  ’?>;  operand2  :=  *?’; 
end; 

end;  (s  1  aa(e-i  aa(e  binary  operations  S) 

f  (operltype  =  ’  i  ’ )  and  (oper2type  =  *t’)  then  bet,  n 
(  si  aa(e- 1 eapl  ate  binary  operations  s) 

if  (operator!  1)  in  [’(*i’(*])  and  (operator(2)  in  (•*’,**’,’•*]) 
t  ben  be(i  n 

if  operator(l)  =  '(*  then 
case  operator[2]  of 

'o'  wri  tel  n  (  Pf  i  I  e,  *  I  aafeTeapei  rcl  ePI  ns  operaodl, 

’ i ’ i operand2i Resul t ,’);’) ; 

'V  :  wri  tel  n  (  Pf  i  I  e,  ’  I  aafeTeapCi  rcl  eMax  (', operaodl, 

'  ,  operand2, '  ■ Resol  t ; 

’*’  :  writelo  (  Pf  i  I  e,  '  1  aafeTeapCi  rcl  eM  n  (’, operaodl, 

operaod2,  Resul  t ,’);’) ; 

end;  (s  case  operator  t) 
if  operator!  I)  =*[’  then 
case  operator(2]  of 

’V  ;  writelo  (  Pf  i  I  e,  >  I  safeTeapSquareMax  (’, operaodl, 

’  >  ’  ,  operand2,  ’ ,  * ,  Resul  *«■);*); 

:  writelo  (  Pf  i  1  e,  ’  I  aateTeapSquareM  n  (’, operaodl, 
ope rand2, Resul t, ’);’) ; 

end;  (S  case  operator  s) 
end  (s  if  operator  then  s) 
el  se  befi  n 
wri tel  n; 

writelo  (’ERROR  operator  ’.operator,’  does  not  exist  ERROR’) 

wri  tel  n  (’ERROR  for  i  aa(e-teapl  ate  operations  error’); 
operaodl  :=  ’?';  operand2 
end; 

end;  (s  i  ■age-teapl  ate  binary  operations  a) 
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if  (operltype  =  ’t')  and  (oper2type  =  ’t*)  then  be(i  n 
(X  tempi  ate-teapl  ate  binary  operation!  !) 
if  (operator!  l]  in  [>♦*,*-*  ,’*',•/•  ,’ v ] )  or 

((operator(  l]  in  (’(’,’(’])  and  (operator(2)  in  (>♦>,’  *’  ,’*’)) ) 

t  hen  bef  i  n 

if  operatorfl]  =  ’(’  then 
cane  operator!  2}  of 

»♦»  :  wri  tel  n  (  Pf  i  I  e.  ’  TempCi  rcl  ePI  us  (  ’  ,  ope  r  a  ndl ,  * ,  *  ,  o  pe  r  and2 ,  ’  ,  ’  ,  Res  ul  t , 

•  v>  :  wri  tel  n  (  Pf  i  I  e,  •  TeapCi  rcl  eMai  ( * ,  operandl,  * , • ,  operand2>  * ,  ’ ,  Resul  t, ' 

**’  :  wri  tel  n  (  Pf  i  I  e,  ’  TeapCi  rcl  e  M  n  ( '  ,  operandl  >  *  operand2,  Penal  t ,  1 

end;  (x  case  operator  x) 
if  operatorfl]  =  *(•  then 
case  operator(2)  of 

’V’  :  wri  tel  n  (  Pf  i  I  e,  ’  TempSquareMax  (’,  o  per  andl optr  and2,  Pesul  t ,  1 

■  wri  tel  n  (  Pf  i  I  e,  *  TempSquareMi  n  ( ' ,  o  per  a  ndl,  *  ,  ’  ,  oper  and2 ,  *  ,  ’  ,  Resul  t .  1 

end;  (X  case  operator  X) 


if  (operator(l]  in  (’(’,’ f  *])  =  f  al  se  then 
case  operatorfl]  of 

>  ♦>  :  wri  tel  n  (  Pf  i  I  e,  1  TempAdd  (’,  operandl)  >,'i  operand2.  ’ ,  Resul  t ,’)  ; 

•->  ;  wri  tel  n  (  Pf  i  I  e,  •  TeapSubt  r  ac  t  (*,  operandl,  •  ,  operand2.  *,’»  Resul 

•X'  :  wri  tel  n  (  Pf  i  I  e, 1  TeapMul  ti  pi  y  (’,  operandl,  ’  .  ■,  operand2,  Resul 

’ /’  :  writeln  (  Pf  i  I  e,  ’  TewpDi  »i  de  (’,  operandl, oper  and2, Resol  t , 

■  v»  :  writeln  (  Pf  i  I  e,  '  TeapMax  (’,  o  per  andl  oper  and2  Resul  t ,’) ; 

>*>  :  writeln  (  Pf  i  I  e,  •  TeapM  n  (’,  oper  andl o  pe  r  and2 Resul  t ,’)  ; 

:  writeln  (  Pf  i  I  e,  operand2,  '  :=  oper  andl 1 ) ; 
end;  (x  case  operator  X) 


t ,»>;•) ; 
t ,•);■) ; 


el  se  bedi  n 
wri  tel  n; 

writeln  (’ERROR  operator  operator,’  does  not  exist 

writeln  (’ERROR  for  t  eapl  at  e- 1  e«pl  a  t  e  operations  ERROR’); 
operandl  :=  ’  ?’ ;  operaodZ  :=  ’  ?’ ; 
end; 

en  ;  (X  template-template  binary  operations  X) 


ERROR’ ) ; 
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if  <  oper  1  type  =  *  i  * )  and  (o  per  2  type  =  ’a’)  then  begin 
(t  i  aage-acal  ar  binary  operations  s) 
if  operator(l]  in  (**•,’/•  =’,*>’  ] 

then  case  operator(l)  of 

's'  if  ( 1  east  h<  operator)  >  l)  and  (operator(2]  =  's') 
t  be  n  be  f  t  n 

writeln  (  PFi  I  e,  ’  const  I  aage  (',  Resul  t operandl, ’.  row,  ’  , 

o pe randl , ’ . col ,  ’  , ope rand2, •);•); 

write  (  Pf  i  I  e,  I  ndent) ; 

wri  tel  n  (  Pf  i  I  e, '  I aa(e  Exponent  (’, o  pe  randl,’, ’,  Re  s  ul  t, ’,’,  Re  suit, ’);’); 
end 

else  begi  n 

wri  tel  o  (  PTi  I  e,  •  Cons 1 1  mage  <  •  ,  Res u!  t ,  •  ,  ’  ,  oper andl ,  >  .  r ow,  ’  , 

ope  randl,' .col  ,’  .operand!,  '  );  1  ) ; 

write  (  Pf  i  I  e,  I  ndent ) ; 

wri  id  o  (pfile.’lmageHultiply  (', o  pe  randl, Re  s  ul  t, Re  suit, 
end; 

'  /'  begi  n 

wri  tel  n  (  PFi  I  e,  '  const  I  mate  (',  Re  sul  t oper  andl row,  ’  , 

operandl, '.col  ,',operand2, 

write  (  Pf  i  I  e,  I  ndent ) ; 

writeln  (  Pf  i  I  e,  ’  I  mageDi  vi  de  (’, o  pe  randl, ',',  Re  suit,’,',  Re  suit, 
end; 

ot  her  wi  se  begi  n 

writeln  (  PFi  1  e,  '  Const  I  mage  (',  Resul  t operandl,  row,  ’  , 

ope  randl, ’.col,’, operands, 

wri  te  (  pf  i  I  e,  I  ndent) ; 

writeln  (  PFi  I  e,  ■  Chart  mage  <’,  operandl, 

’,''', o  pe  rator, ''',’>  Re  suit,’,’,  Re  suit, 
end;  ( t  case  el  se  S) 
end  (S  case  operator  S) 


else  befi  n 
wri  tel  n; 

writeln  ('ERROR  operator  ', operator,’  does  not  eii  st  ERROR’); 
writeln  ('ERROR  for  iaa(e-scalar  operations  ERROR’); 
operandl  :  =  •?';  ope  ran  d2  :  =  '?'; 
e  nd; 

end;  (S  i  made-seal  ar  binary  operations  t) 
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if  (operlt  ype  =  ’t’)  aid  (oper2type  =  ’s')  tben  bed  a 
(s  template-scalar  binary  operations  t) 
if  operator[l]  in  [’s'] 
tben  case  operator(l)  o t 

:  writela  (  Pf  i  1 1,  ’  TeapScal  ar  Mul  t  i  pi  y  ( 1 .  operandl ,  *  ,  1  ,  operand2,  *  , 
end  (s  case  operator  a) 
el  se  be  a  i  n 
wri  tel  n; 

writeln  (’ERROR  operator  ', operator,’  dons  not  exist  ERROR’); 
writeln  ('ERROR  for  template-scalar  operations  ERROR’); 
operandl  :  =  ’  ?’ ;  operand!  :  =  >7’; 
end; 

end;  (s  template-scalar  binary  operations  S) 

if  (operltype  =  ’s’)  and  (oper2type  =  ’s’)  tben  bed  n 
(s  scalar-scalar  binary  operations  X) 
if  operator  -  ’  :  =  • 

tben  writeln  (  Pf  i  I  e ,  oper  and2 ,  1  :=  ’  ,  operandl,  ’  ;  ’  ) 
else  bed  n 
wri  tel  n; 

writeln  ('ERROR  operator  *, operator,’  does  not  exist  ERROR’); 
writeln  ('ERROR  for  scalar-scalar  operations  ERROR’); 
operandl  :  =  ’  ?’  ;  operand2  :  =  '7*; 
end; 

end;  (X  scalar-scalar  binary  operations  X) 
end;  (t  if  operands  <>  7  tben  X) 

AssignResults  (o  pe  randl, epe  ran  d2, operator.  Re  salt, ex  press,  no  deptr); 
end;  (X  procedure  Bi  naryOperati  on  x) 

procedure  Unar  yOper  at  i  o  n  (  var  eipress  :  Tree;  nodeptr  :  inte(er); 

(X  write  t  be  unary  operation  procedure  calls  to  t  be  output  file  X) 
var  operator,  operandl,  operand2,  Resul  t  ;  NameStr; 

oper  1 1  ype,  operxt  ype  :  ebar; 
bed  n  (x  procedure  Unar  yOperat  i  on  X) 

(X  retrieve  the  operand  and  operator  f  ron  the  binary  tree  x) 
operator  :  =  express) .  nodept  r.  ) .  val  ue; 
if  express),  nodeptr.  )■  I  c  hi  I  d  <>  0 
then  operandl  e  x  pr  ess( .  expr  e  a  s  ( .  node  pi  r .  ) .  I  c  hi  1  d.  ) .  val  ue 
else  operandl  :  ^  ex  press!  .  e  x  pre  s  s(  .  no  de  pt  r .  ) .  r  c  hi  I  d.  ) .  val  ue ; 
operand!  : =  *  ' ; 


I 

I 

’,Result,’);’)|  i 
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(x  deter*  oe  the  operand  type  x) 


operltype  :  =  '  o per 2 type  :=  ’  '  ; 

for  i  :  =  1  to  NumI  HfeNaMi  do  if  operand!  =  I  mage  It  s  t  (  .  i  .  )  then  operltype  ’  i  '  ; 
for  i  :  =  1  to  NumTempI  ateNames  do 
if  operandl  =  Tempi  at  e  Li  st  ( .  i  .  )  then  operltype  :=  *  t  •  ; 
for  i  :  =  i  to  NumScal  arNames  do  if  operandl  =  Seal  a  r  Li  s  t  (  .  i  .  )  then  operltype  :  =  ’s'; 
if  operltype  =  *  '  then  begin 
oper l t ype  : =  '  s* ; 

for  i  :=  1  to  I  eng  t  h(  operandl  )  do 
if  (operltype  =  ’s')  and  (operandl(i]  in  [  * ,  *  O’  .  .  1  9*  ] ) 
t  hen  oper  1 1  ype  :  =  *  s' 
el  se  oper It  ype  : =  *  * ; 
end;  (X  if  operltype  then  x) 
if  operltype  =  *  '  then 

wn  tel  n  (  *  ERROR  operandl.1  is  not  an  image,  template,  or  scalar  operand  ERROR') 

if  operandl  <>  » ?'  then  begin 
write  (  Pf  i  1  e ,  I  ndent ) ; 

as  si  gnl  ntermedResul  t  (operltype,  operator,  Re  sal  t ) ; 
case  operltype  of 

*  i  *  begin  (x  unary  image  operations  X) 
if  operator  =  »!'  then 

writeln  (  Pf  i  I  e,  *  Abs  val  I  mage  (',  operandl,  Resul  t  ,');*) ; 
if  operator  =  »-»  then  begin 
writeln  (  PFi  1  e,  '  Const  I  mage  (', Result, 

'  *  ',  operandl,  * .  row,  operandl,  » .  col  ,  -1) ;  »  ) ; 
wri  te  (  Pf  i  I  e,  l  ndent ) ; 

wri  tel  n  (  Pf  i  1  e,  '  I  mageMul  t  i  pl  y  (', oper a ndi Resul t Resul t ; 
end; 

end;  (x  case  unary  image  operations  x) 

't»  begin  (*  unary  template  operations  X) 
if  operator  =  *  i  '  t  hen 

writeln  (  Pf  »  1  e,  '  Abs  Val  Temp  (*,  operandl, Resul  t ; 
if  operator  =  1  - '  then 

wri  tel  n  (  Pf  i  1  e,  '  Temp  sc  a  I  arHuI  tipi  y  (*,o  pe  randi, *, -l, Re  suit, ');*); 
end;  (x  case  unary  template  operations  s) 
end;  (x  case  operand  type  x) 
end;  (X  if  operandl  <>  ?  then  x) 

As  si  gnResul  ts  ( operandl, operand2,operator,  Resul  t,  express,  nodeptr); 
end;  (x  procedure  Unar yOper at  i  oo  X) 
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procedure  IAParser  (buffer  :  BufferStr;  var  express  :  Tree; 

ear  nodeptr  :  integer); 

(X  pane  the  I A  expressions  and  build  a  bi  aary  tree  of  operaada  and  operator!  t) 
var  i,j,next  :  i  nteger; 
letter  :  cbar; 

assi  gnoper,  absval  f  I  ag,  EndOfOperand  :  boolean; 
bed  i  (*  procedure  IAParser  s) 
absval  f  1  ag  :  -  f  al  se; 

(t  if  an  assignment  exists  in  the  expression!  evaluate  the  right  side  s) 
assi  gnoper  false; 
if  i  ndex(  buf  f  er,  ■ :  =’ )  <>  o 
t  ben  begi  n 
assi  gnoper  :  =  true; 

I  APar  ser(  subst  r(  buf  f  er,  i  ndex<  buffer,  * :  -•  )  *2, 1  engt  h(  buffer ) -i  adex(  buf  fer,':=*)-l). 
express i nodeptr); 

end  (s  then  x) 

else  begin  nodecnt  :=  1;  nodeptr  :  =  i;  end; 

i  :  =  1 ; 

if  assignoper  =  false  then  begin  (X  translate  the  1A  expression  x) 
while  i  <  I  engt  b(  buf  f  er)  do  begin 
1 etter  : =  buf f er( i ] ; 
next  : =  i  ♦  l; 

if  (letter  in  ('  a* u*  ,*  w* ..*  x*  ,*  0* ..*»']  )  or 

((letter  =  •->)  and  (  buf  fer[  next  ]  in  l'0». .'»»))) 
then  begin  (X  image,  template,  scalar,  or  variable  operand  x) 
nodecnt  ;=  nodecnt  »  l; 
if  express(  .  nodeptr. ).  val  ue  =  *’ 
then  express(  .  nodeptr. ) .  I  chi  I  d  :=  nodecnt 
else  express! .  nodeptr.  ).  rc  hi  I  d  :  =  nodecnt; 
express! .  nodecnt .).  parent  :  =  nodeptr; 
j  :  =  next ; 

EndOfOperand  :=  false; 

while  (j  <  I  engt  h(  buf  f  er) )  and  (EndOfOperand  =  false)  do 
if  (  buffer!  j  1  in  [•  a’ ..»  u> ,’ w*  0*  ..’»•]  )  or 

( (  buf  f  er[  j  ]=’.’)  and  (bufferU»l]  in  ('o*.  .'*>))) 
t  hen  j  :  -  j  *  l 
else  EndOfOperand  :=  true; 
express! .  nodecnt.  ).  val  ue  :  =  s  ubs  t  r(  buf  f  er ,  i  ,  j -i  ) ; 
i  :  =  j ! 

express! .  nodecnt.  ).  I  chi  I  d  :=  o; 
express! .  nodecnt .).  rc  bi  I  d  o; 

end;  (X  if  image,  template,  scalar,  or  variable  operand  X) 
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A  -aJU  k  aJjds  ^  -*v- 


if  letter  «  •(’  thee  be(i  a  (*  start  pareatbeeii  S) 
nodecot  :  =  nodeent  ♦  1; 
if  express! .  nodept  r.  ) .  »al  ne  =  •  ’ 
then  eipreae( .  eodeptr.  ).  I  chi  1  d  :  =  nodeent 
el  ee  exprete( .  no  dept  r.  ) .  rchi  I  d  :  <  nodeent; 
express! .  nodeent .).  parent  :  >  nodeptr; 
nodeptr  :•  nodeent; 
i  ;  =  i  »  1 ; 

end;  (*  if  start  parenthesis  s) 

if  letter  *  *)*  then  be(i  a  (t  finish  parenthesis  t) 
if  express! .  nodept  r.  ).  parent  <>  o 
then  nodeptr  :  =  express!  .  nodeptr.  ).  parent 
el  se  beti n 

nodeent  :  =  nodeent  ♦  1; 

express! .  nodept  r.  ).  parent  :=  nodeent; 

express! .  nodeent.  ).  I  chi  I  d  :  =  nodeptr; 

nodeptr  :  =  nodeent; 

end; 

■  :  =  i  *  l ; 

end;  is  if  finish  parenthesis  s) 

if  fletter  in  [■  <  *  .  •  I  * .  *  ♦  *  .***.•/•,  •  .*  .**•,’ *  •  ,•:*,•='  ,•<•,»>•  J )  or 
filetter  =  •-•)  and  !  buff  er(  next  ]  in  [’ a*  u‘ w*  j‘ ,’(’)) ) 
then  begin  is  algebraic  operator  encountered  S) 
if  express! .  nodeptr.  ).  *al  ue  <>  *  *  then  beti  n 
is  chained  algebraic  operator  encountered  -  set  a  new  node  s) 
nodeent  :»  nodeent  ♦  t; 
if  express! .  nodeptr. ).  rchi  I  d  <>  0 
then  begi  n  is  chained  binary  operator  encountered  t) 
if  express! .  nodept  r.  ).  parent  <>  o 
then  beti  n  is  chained  operator  in  parentheses,  alter  pointers  s) 
Is  place  the  new  node  between  the  present  and  parent  node  S) 

Is  link  the  new  node  to  the  present  node's  parent  S) 
if  express!  .  express! .  nodept  r.  ).  parent .).  I  chi  1  d  =  nodeptr 
then  express! .  express!  ■  nodept  r.  ).  parent  .)•  I  chi  I  d  :  =  nodeent 
else  express! .  express! .  nodept  r.  >•  parent .).  rchi  I  d  :*  nodeent; 
express! .  nodeent .).  parent  :  =  express!  .  nodept  r.  ).  parent ; 
end;  It  then  alter  pointers  S) 

IS  link  the  new  node  to  the  present  node  s) 
ex  pres  s!  .  no  de  pt  r .).  parent  :  -  nodeent; 
e i pr ess!  .  nodec nt .  ) .  I  c hi  I  d  :  =  nodeptr; 
end  is  chained  binary  operator  s) 
else  befi  a  (t  chained  unary  operator  encountered  S) 

■  f  express!  .  nodept  r.  ).  »al  oe  =  •' 
then  express!  .  node  pt  r .  ) .  I  c  hi  I  d  :  -  nodeent 
else  e  i  pr  ess!  .  nodept  r .  ) .  r  c  hi  I  d  :*  nodeent; 
express!  .  nodeent .),  parent  :  :  nodeptr; 
end;  is  chained  unary  operator  S) 


(a  n*>ve  *0  the  new  node  *) 
nodeptr  :•  nodecnt; 

end;  (t  chained  algebraic  operator  eaconatered  a) 


(a  place  the  operator  Into  the  eapty  node  a) 
if  express! .  nodept  r.  ).  val  ae  =  ••  then  begin 
If  (letter  *  ’[’)  or  (letter  =  '(') 
t  hen  begi  a 

eipreul .  nodeptr.  ).  «nl  tie  :»  anbntr(  buf  f  er,  i  ,  3) ; 
i  =  I  *  3;  (a  a  hi  p  the  operator  charatera  a) 
end  (a  if  letter  then  a) 
el  ae  begi  n 

if  ((letter  =  'a')  and  ( buf  fer(  next  ]  *  *a>))  or 

((letter  in  {’>’.*<•])  and  {  buf  f  erf  neat )  in  p*», •>*})> 

t  hen  begi  n 

expreaa(  .  nodeptr.  ).  val  ue  :«  aubat  r(  buffer,  i  ,  2) ; 
i  :  =  i  ♦  *; 

end 

el  ae  begi  n 

if  absvalflag  =  falae 
then  exprean( .  nodept  r.  ) .  val  ue  :  -  letter; 
i  :=  i  ♦  1}  (a  ship  the  operator  character  a) 
if  (letter  =  *!’)  and  (absvalflag  =  false) 
then  absvalflag  :=  true 
else  abaval  fl  ag  :*  falae; 
end;  (a  if  letter  else  a) 
and;  (a  if  letter  else  a) 
end;  (a  place  operator  into  the  eapty  node  a) 
end;  (a  if  algebraic  operator  a) 

if  letter  =  '  '  then  i  :  =  i  ♦  l; 

end;  ( a  whi  lei  do  a) 

(a  return  to  the  top  node  a) 
while  expresa( .  nodept  r. ).  parent  <>  0  do 
nodeptr  :  =  express!  .  nodept  r.  ).  parent ; 
end;  (a  if  aasi  gnoper  false  a) 

(a  if  assign  operator  exists,  place  it  at  the  top  of  the  tree  a) 
if  ansi  gnoper  *  true  then  begin 
nodecnt  :  =  nodecnt  *  l; 
express! .  nodeptr. ).  parent  :■  nodecnt; 
express( .  nodecnt .  ) .  I  chi  I  d  nodeptr; 
express!  .  nodecnt .).  val  ue  :  = 
nodeptr  :  =  nodecnt; 
nodecnt  :  «  nodecnt  t  i; 
express!  .  nodept  r.  ).  rchi  I  d  nodecnt; 
express! .  nodecnt .).  parent  :  =  nodeptr; 

express!  .  nodecnt  .).value  :  *  su  bstr(  buffer,  l.index(  buffer,  • :  =•)-!); 
end;  (a  if  asaignoper  then  a) 
end;  (a  procedure  lAParser  a) 


procedure  IATraaolator  ( var  expreaa  :  Tree;  aodeptr  :  integer); 

(a  traaal  ale  t  be  IA  expressions  by  evaluating  the  bi  aary  tree  a) 
begin  (a  lATraaalator  a) 
if  eapresa(  .  aodeptr.  )  .  I  chi  I  d  <>  0 
the  a  I  ATraaal  atorl  eipresa,  eipressl .  nodept  r .  ) .  I  chi  I  d)  ; 

if  express!  .  aodeptr.  ).  rchi  I  d  <>  o 
t  be  a  I  ATraaal  ator(expresatea  preset,  aodeptr.  ).  rchi  I  d); 
if  express! .  aodept  r.  ) .  val  ue  <>  *■ 
tbea  begin  (a  execute  t  be  operator  a) 
if  (( express!  .  aodept  r.  ).  1  chi  I  d<>0)  aad  (express! .  aodeptr.  ).  rcbi  1  d=0) ) 
((  express! .  nodept  r.  ).  I  chi  I  d=0)  aad  (express!,  aodeptr.).  rchi  ldoo)) 
tbea  UnaryOperati  on  ( express,  aodept  r) ; 
if  (  express!  .  aodept  r.  ).  1  chi  I  d  <>  0)  aad  (  ex  press! .  aodept  r .).  re  bi  I  d  <> 
then  Bi  naryOperat  i  oa  ( express,  aodept  r>; 
end  (a  operator  execution  a) 
else  if  express! .  aodeptr. ).  parent  <>  0  then  begin 
(a  reaove  aa  unnecessary  pair  of  parentheses  a) 
if  express!  .  express! .  aodept  r. ).  parent .).  I  chi  I  d  =  aodeptr 
tbea  (a  link  tbe  parent  node's  left  child  to  this  node's...  a) 
if  express! .  aodeptr.  ).  I  chi  I  d  <>  0 
tbea  begin  (a  ...left  child  a) 

express!. express!,  no  deptr.).  pa  rent. ).lc  hi  Id  :  =  express!,  no  deptr.). 
express! .  express! .  aodept  r.  ) .  1  cbi  1  d.  ) .  parent  :  =  express! .  aodept  r.  ) . 
end 

else  begin  (a  ...right  child  a) 

express!  .  express! .  aodeptr. ) .  parent.  ) .  I  cbi  I  d  :  *  express!  .  aodeptr.  ). 

express! .  express!  .  aodept  r.  ) .  rchi  I  d.  ).  parent  :  =  express! .  aodeptr.  )  . 

end; 

if  express! .  express! .  aodept  r. ).  parent .}.  rcbi  I  d  =  aodeptr 
then  (a  link  the  parent  node's  right  child  to  this  node's...  a) 
if  express!  .  nodept  r.  ) .  1  cbi  I  d  <>  0 
then  begin  (a  ...left  child  a) 

ex  pressf. ex  press!,  nodept  r.  ).  parent.  ).  rchi  I  d  :  =  express!,  aodept  r.  ). 

express!  .  express!  .  aodept  r.  ) .  1  cbi  I  d.  ) .  parent  :  =  express!  .aodeptr.  ) . 


I  chi  I  d; 


rchi  I  d; 


el  se  begi  a  (  a 


,  ri  ght  chi  1  d  a) 


express!. ex  pr  east,  no  deptr.).  pa  rent.). rcbi  Id  :  *  express!,  aodeptr.). 
express! .  express! .  aodeptr.  ) .  rchi  1  d.  } .  parent  :  =  express! .  aodeptr.  ) , 
end; 

end;  (a  reaove  parentheses  a) 
end;  (a  lATranslator  a) 

begin  (t  aai  a  prograa  I  aage  _  Al  gebr  a_  Preprocessor  a) 
reset  ( i  nput ) ; 
rewri  t  e  ( out  put ) ; 


I  cbi  I  d; 


rcbi  I  d; 


(s  intialize  variable!  and  set  intermediate  operands  to  false  t> 
i  HI.  exist  :  =  false;  IR2. exist  ;  >  false;  i  M3,  exi  st  :=  false; 

I  R4.  exist  false;  i*S. exist  :  =  false; 

tRi. exist  :  »  false;  tR2. exist  :»  false;  tR3.  exist  :  =  false; 

tR4. exist  :*  false;  tR3.  exist  :»  false; 

baf fer  : =  ■  * ; 

(s  initialise  the  lists  of  i  aage,  template,  and  scalar  variables  t) 

InafeU  st{  .  x.  )  ;=  *  i  Rl  *  ;  I  Bate  Li  st(  .  2.  )  :*  'ins’; 

I  aifeli  st  ( .  3.  )  :=  ’iR3>;  I  aageLi  st( .  4. )  :  =  >i*4>; 

I  male  Li  st( .  S.  )  :  =  *  i  R3» ; 

Kual  naftNanes  :  =  S; 

for  i  :=  Muni  aageNaaese  1  to  Ha  i  N  a  Be  s  do  I  Bate  Li  st  ( .  i  . )  :=  •*; 

Teapl  ate  Li  st( .  l.  )  :  =  ’tRl’;  Teapl  ateli  st( .  2.  )  :  =  'tR2’; 

Teapl  ateLi  st( .  3.  )  :=  1 1«3' ;  Teapl  at  e  Li  st  ( .  4 .  )  :  =  »iR4*; 

Teapl  ateLi  st(  .  5.  )  :  =  MRS'; 

NuaTeapI  ateNaaes  :  =  3; 

for  i  :  =  NuaTeapI  ate  Names  el  to  MaxNaaes  do  Teapl  at  eLi  st  (.  i  .  )  ;  =  >•; 

Seal  ar  Li  st  (  .  l.  )  :  =  ‘dotvai’;  Seal  ar  Li  st  ( .  2.  )  :  =  ’•oval’; 

Seal  ar  Li  st  (  .  3.  )  :=  *  aaxval  ’ ; 

NumScal  arNaaes  ;  =  3; 

for  i  :=  NuaScal  arNaaes*  l  to  MaxNaaes  do  Seal  arLi  st(  .  i  .  )  :  = 

(S  retrieve  an  existing  I  made  AlRebra  file  to  be  translated  S) 
open(  1  nFi  I  e, ’ t  ransl  at .  i  a*  ,  unknown) ; 
reset  ( I  nFi  I  e)  ; 

(s  chanfe  all  letters  in  the  file  to  lower  case;  s) 

(<  remove  any  comments  and  blank  lines  from  the  file  *) 
open(  TapFi  I  e, ’ t  ransl  at.  tap1 ,  unknown)  ; 
rewri  te  (  TapFi  I  e) ; 

wri  tel  n  ('  ...  t ransl  ati  nf  letters  to  lower  case;  reaovi  at  coannents  >); 

readl  n  (In  FI  I  ei  buffer) ; 
while  (eofflnFi  I  *)  =f  al  se)  do  befi  n 
(*  remove  comments  from  the  line  s) 
while  i  ndex(  buf  f  er,  1  (  s’ )  <>  a  do  be(i  n 
buffer  :=  '  *  ♦  buffer  ♦  1  ’;  (s  required  for  proper  operation  of  substr  s) 
if  i  ndex(  buf  f  er,  •  s)  •  )  <>  a 

then  buffer  :*  subs  t  r  (  buf  f  er ,  l,  i  ndex(  buf  f  e  r(  ’  (  s* )  - 1 )  ♦ 
substr(bufrer,index(buffer,,t),)t2,lendth(buffer)-index(buffer,,x),)-i) 
else  buffer  :=  subs  t  r(  buf  f  er  >  1 ,  i  ndex(  buf  f  er,  *  (  s'  )  - 1 ) ; 
end;  (s  while  comemnt  do  s) 

(S  remove  trailing  blanks  s) 
i  =  I  en{t  h(  buffer) ; 

while  (i  >  1)  and  (buffer(i)  =  '  1 )  do  i  :  =  i  -  1; 
buffer  substrl  buf  f  er,  1,  i  ) ; 

(s  translate  all  letters  to  lower  case  s) 
for  I  :  =  1  to  I  enft  b{  buf  f  er)  do 

if  ord(  buf  f  er(  i  1 )  in  (63. .90)  then  buffer(i)  :=  chr(  or  d(  buf  f  er  [  i  I )  ♦  32) ; 


(*  print  Don-bl  ank  lines  to  TapFi  I  e  S) 
if  i  >  1  then  writeln  (  TapFi  I  e,  buf  f  e  r ) ; 
readl  n  ( I  nFi  I  e>  buf  f  er) ; 
end;  (s  while  not  tot  t ) 

(s  write  the  last  line  s) 
writeln  (  TapFi  1 1,  buf  f  er) ; 
close  (  TapFi  1  e); 
cl  ose  (  I  nFi  I  e) ; 

(s  build  the  PASCAL  file  from  '  t  r  ansi  at .  t  ap'  t) 
open(  I  nFi  I  e>  1  transl  at.  t  ap’  ,  unknown) ; 
reset  ( I  nFi  I  e) ; 

open!  pf  i  I  e,  >  t  ransl  at .  pas'  ,  unknown) ; 
rewri  te  (  Pf  i  I  e) ; 

(S  find  the  profraa  name  if  it  exists  s) 
readl  n  ( I  nFi  1 1,  buf  f  er) ; 

if  ( i  ndei(  buf  fer.  *  const  ’)=<>)  and  ( i  ndex(  buf  f  er, ’ t  ype  ■  )-0)  and 
( i  ndex(  buf  f  er  ■  ’  var  *  )  =  0 )  and  ( i  nder<  buf  fer,  1  bed  n  *)=0) 
then  Nant  :  =  buffer 
el  se  Name  :  =  '  unnasad*  ; 

(s  determine  the  data  type  of  the  i  aafe  and  teapl  ate  operands  t) 

IaafeType  =  Def  aul  tl  aafe; 

Teapl  ateType  :=  Def  aul  tTeapI  ate; 

(t  search  the  declaration  section  s) 
while  i  ndex(  buf  f  er.  *  bed  n1  )  -  0  do  bed  » 
readln  ( I  nFi  I  e.  buf  f  er) ; 
if  1  ndex(  buf  fer,  ’  type  » )  <>  0  then 
(X  find  the  data  type  declarations  if  they  exist  t) 

while  (( i  ndex(  buf  fert  •  const  ’)  =  0)  and  {  i  ndex(  buf f er> ’ var  ’  )  =  0)  and 

( i  ndex(  buf  f  er>  '  befi  o' )  =0))  do  beds 
(X  deterai  ne  the  data  type  of  the  i  sites  and  templates  s) 
if  i  ndex(  buf  f  er ,  '  i  t  ype* )  <>  0  then  bed  n 
if  i  ndex(  buf  f  er,  ’  i  ntefer'  )  <>  0  then  IaafeType  :=  'intefer'i 

if  i  ndex(  buf  fer.  ’  real  ’  )  <>  o  then  IaafeType  :<  'real*; 

end;  (s  if  itype  then  x) 
if  i  ndex(  buf  f  er ,  ’  1 1  ype*  )  <>  0  then  befi  n 
if  i  ndex(  buf  f  er,  •  i  ntefer*  )  <>  0  then  Teapl  ateType  'intefer’; 

if  i  ndex( buf fer( ’ real  1  )  <>  0  then  Teapl  ateType  :=  'real*; 

end;  (X  if  ttype  then  X) 
readln  (  I  nFi  I  e.  buf  f  er) ; 

end;  (x  while  not  'const  * , '  var  •  >  or  '  befi  n'  X) 
end;  (x  while  not  •  bef  i  n'  x) 

(x  include  the  proper  i  aaf  e  al  febra  and  input/output  routines  X) 
if  (IaafeType  =  'intefer')  and  ( Teapl  at  eType  3  'intefer') 


t  be  d 

wri tel 

n  (  pf  i  1  e,  ’  ( i  nhen  t  (  ' 

•  i  i  ao per .  eov'  *  *  i  i  o.  env* 

*)) 

■ ) 

el  se 

wr i  tel 

n  (  Pf  i  1  e , '  ( i  nhe r i  l  ( ’ 

*  ri  aoper.  env'  *  »  ri  o.  env* 

•  >  i 

’ ) ; 

wr  i  t  el 

o  (  Pf  i 

1  e,  '  prof  raa  ’  ,  Name,  ' 

( i  nput  t out  put ) ; 

'  ) ; 

B-  1  5 


wri  i«li  ( 1 


.  procnsi  o(  user-defined  types,  coutitlii  and  variables'); 


(t  find  and  append  aay  user  defined  constants  here  t) 
reset  ( I  nFi  1  e) ; 

(s  search  the  entire  declaration  section  for  constants  t) 
read!  n  ( I  nFi  I  e,  buf  f  er ) ; 

while  i  ndex(  buf  fer,  '  befi  n*  )  =  o  do  befi  n 
read)  n  ( I  art  I  e,  buf  for) ; 
if  i  ndei(  buf  fer,  •  const  ’)  <>  o  then 
(t  process  each  additional  constant  declaration  line  s) 
while  ( ( i  ndei(  buf  f  er. ’ t  ype  ')  *  0)  and  <  i  ndez(  buf  f  er ,  '  var  '  )  *  0)  and 

( i  udei(  buff  er» '  befi  a*  )  *  0))  do  befi  a 
(s  strip  the  reserved  word  'type'  f ros  the  buffer  s) 
if  i  adex(  buf  f  er,  1  const  ')  <>  0 
then  buffer  :  =  s  ubst  r(  buf  fe  r .  i  nde*(  buf  f  er  ■  '  const  ')*5, 

I  enf  t  h(  buf  fer)  -i  ndei(  buffer,  'const  ’)-«); 

(s  move  all  blanks  froa  the  line  s) 
ol  dbuffer  : =  buffer; 
buffer  : =  ’ ' ; 

for  i  l  to  I  enft  b(  ol  dbuffer)  do 
if  oldbuffer(i)  <>  '  '  then  buffer  :  =  buffer  ♦  ol  dbuf  f  er[  i  ) ; 

(s  write  the  const  declaration  as  is  and  ensure  each  line  ends  with  a 
if  i  ndexf  buf  f  er,  ' ;  '  )  =  0  then  wri  tel  n  (Pfile,’  const  '.buffer,';’) 

else  writeln  (Pfile,'  const  '.buffer); 

readl  n  ( I  nFi  I  e,  buf  f  er ) ; 

end;  (S  while  not  ’  t  ype'  ,  '  var* ,  or  ’  befi  n* ,  do  t) 
end;  (S  while  not  ’  befi  a*  do  t) 

(t  find  and  append  aay  user  defined  types  here  S) 
reset  ( I  nFi  I  e) ; 

(S  search  the  entire  declaration  section  for  types  S) 
readl  n  ( I  nFi  I  e,  buf  f  er ) ; 

while  i  ndex(  buf  f  er,  '  befi  o'  )  =  0  do  befi  n 
readl  a  ( I  nFi  I  e,  buf  f  er ) ; 
if  i  ndei(  buf  f  er, ’ t  ype  ')  <>  0  then 
(t  process  each  additional  type  declaration  line  s) 

while  (( i  ndei(  buf  fer,  '  const  ')  =  0)  and  ( i  ndez(  buf  f  er,  ’  var  '  )  =  0)  and 

( i  ndez(  buf  fer ,' befi  n' )  =  0))  do  befi  n 
(s  write  the  type  declaration  as  is  if  it  is  not  i  t  ype  or  ttype  t) 
if  ( i  ndei(  buf  f  er,  ’  i  t  ype'  )  =  0)  and  ( i  nde  K  buf  f  e  r ,  ’  1 1  ype’  )  =  0)  then  bef  i  n 
(s  strip  the  reserved  word  'type'  froa  the  buffer  a) 
if  i  ndes(  buf  f  er, ' t  ype  •)  <>  0  then 
buffer  :=  s  ubs  t  r  (  buf  f  e  r ,  i  nde  i(  buf  f  er ,  ■  t  ype  ’)tt, 

I  enft  h(  buffer)-i  ndeii  buffer, ’ t  ype  •  )-J)  ; 

(»  reaove  all  blanks  froa  the  line  S) 
ol  dbuf  fer  : -  buffer; 


for  i  :  =  I  to  I  enf  t  blot  dbuf  fer)  do 
if  oldbuffer(i)  <>  •  •  then  buffer  buffer  *  ol  dbof  f  e  r  I  i  1  ; 
if  i  ndei(  buf  f  er ,  ’  ;  '  )  *  0  then  writeln  (Pfile,'  type  '.buffer,';') 

else  writeln  (pfile,'  type  ', buffer). 


& 


end;  (s  if  type  declaration  X) 
readl  a  ( 1  nFi  I  e ,  buf  f  er ) ; 

end;  (x  while  not  ’coml  ■ ,  ’ear  '  ,  or  '  be(i  n*  do  t) 
end;  (t  while  not  ’  begi  n*  t) 

(t  find  and  append  any  oxer  defined  variable!  here  X) 
rexet  ( I  nFi  I  e) ; 

(X  inarch  the  entire  declaration  aection  for  vari  abl  ex  x) 
read)  n  ( I  nFi  I  e,  buf  f  er ) ; 
while  i  ndex(  buf  f  er,  ’  befi  n*  )  =0  do  befi  n 
readl  n  ( I  nFi  I  e.  buf  f  er ) ; 
if  i  ndex(  buf  f  er  i  ’  var  ’)  <>  0  then 
(X  process  each  additional  constant  declaration  line  X) 

while  (( i  ndex( buf fer, *  const  ■)  =  o)  and  ( i  ndex(  buf  f  er ,  *  t  ype  •)  -  o)  and 
( i  ndei(  buf  f er ,  *  befi  n*  )  =  0))  do  befi  n 
(X  strip  the  reserved  word  *vai'  frost  the  buffer  s) 
if  i  ndex(  buf f er, ’ var  ')  <>  0 
then  buffer  :  =  subst  r(  buf  f  er,  i  ndex(  buf  f  er,  1  var  *)v3, 

len(th(  buf  fer )  -  i  ndex(  buffer,  ’var  ’  )  -2 ) ; 

(X  renove  all  blanks  froa  the  line  X) 
oldbuffer  :=  buffer; 
buffer  : =  ’  ’  ; 

for  i  :  =  l  to  I  enft  h(  ol  dbuf  f  er )  do 
if  oldbuffer(i)  <>  1  1  then  buffer  :  :  buffer  t  ol  dbuf  f  er  [  i  ) ; 

(x  write  the  variable  declaration  as  is  and  ensure  each  line  ends  with  a  s) 

if  i  ndex(  buf  f  er,  1  ;  ■  )  =  o  then  wri  tel  n  (Mile,’  var  ’  ,  buf  fer,  * ;  * ) 
else  writeln  (Mile,*  var  >, buffer); 

(X  if  the  line  contains  any  i  aafe,  teaplate,  or  scalar  declarations,  s) 

(s  append  these  naans  to  the  appropriate  naan  list  X) 


if  ( i  ndex(  buffer,  ' 

:  1  aa(e  ’ 

) 

<> 

0) 

or 

( i ndex( buffer, ' 

i  aafe' 

) 

<  > 

0) 

or 

( i  ndex( buffer, ' 

;  teapl  ite* ) 

<> 

0) 

or 

( i ndex( buffer, ’ 

teapl  ate' 

1  ) 

<  > 

0) 

or 

( i  ndex(  buffer,  ’ 

:  I  nte(er’ 

) 

<> 

0) 

or 

( i  ndexf  buffer,  ’ 

i ntefer ’ 

) 

<  > 

0) 

or 

( i  ndex( buffer, ' 

:  real  * 

) 

<  > 

0) 

or 

(1  ndexl  buffer,  ’ 

real  • 

) 

<  > 

0) 

t  boa 

(X  deterai  ne  the  type  of  variable  declarations  X) 

Decl  aredType  :  =  subs t  r(  buf  f  er ,  i  nde i(  buf  f  er,  ’ :  ’  )  ♦  1 ,  «) ; 

(X  extract  the  variable  naaes  froa  the  buffer  X) 

i  ;  >  l; 

while  i  <  i  ndex(  buf  f  er,  ’ :  1  )  do  befi  n 

J  ■  i  ; 

while  (buffer(j]  <>  and  (bufferfj)  <>  ’,’)  do  j  :  -  J  *  1 ; 

if  ( j-  i  )  >  0  then  be  (in 
if  bed  aredType  =  ’  i  Bag’  then  be(i  n 
Nuad  aafeNaaes  :=  Nual  aafeNaaes  »  i; 

I  aa(e  Li  st  ( .  Nual  aafexants.  )  subst  r(  buf  f  er,  i  ,  i -i  ) ; 
end;  (X  if  ■  aa(e  declarations  then  X) 
if  Decl  aredType  -  'leap'  then  be(i  n 
NuaTeapI  ateNaans  :=  NuafTeapI  atevaaes  ♦  i; 

Teapl  ate  Li  st( .  NoaTeapI  ateNaars.  )  :  =  subst  r(  buf  fer ,  i  ,  j -i  ) ; 
end;  (X  if  teaplate  declarations  then  t) 
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if  ( Decl aredType  >  ’into’)  or  ( D«el  aredType  =  'real')  tbeo  b«(i  a 
Nuascal  arNaaes  :•  NuaScal  irHaati  *  l; 

Seal  ar  Li  st( .  NnaScal  arNaaes.  )  :*  subat  r(  bof  f  tr,  i  ,  j-i  ) ; 
aad;  (s  if  scalar  declarations  thaa  t) 
cad;  (t  if  j-i  t  bca  x) 
i  :  *  i  *  1; 

cad;  (x  while  buffer(i)  do  x) 
eod;  (X  if  I  Bate,  template,  or  xcal  ar  declaration  X) 
r cadi  a  ( I  an  I  »,  baf  f  cr) ; 

cad;  (x  while  aot  *  coast  •  ■  1 1  ype‘ ,  or  '  btfi  ■' ,  do  X) 
cad;  (x  while  aot  *  begi  a’  do  X) 


(x  traaslate  the  ia  file  to  PASCAL  procedure  calls  x) 

(x  find  the  begi  aai  ng  of  the  executable  code  sectioa  x) 
while  i  adei(  buf  f  er,  '  begi  a‘  )  =0  do  readl  a  ( I  a  Fi  I  e,  buf  f  e  r ) ; 


wr i  tel  n 

(  PFi  1  e,  • 

begin  (X  prograa  1  t  Naa,  1 

X) 

•  ) 

wri  tel  n 

(  Pf i 1  e,  * 

reset  ( i nput ) ; 

’) 

wr i  tel n 

(  Pf i 1  e,  » 

rewri  te  (  oat  put ) ; 

•  ) 

(X  fi  ad  aad  repositioa  any  i  avari  ant  teaplate  definitions  here  X) 
while  eof  ( I  nFi  I  e)  -t  al  se  do  begin 

while  ( ( I  adex(  buf  f  er  > 1  i  avari  ant 1  )  =0)  or  ( i  ndex(  baf  f  cr ,  *  teapl  at  e‘ )  =0) )  and 
(eof  ( I  nFi  I  e)  =f  al  se)  do  readl  a  ( I  nFi  I  e,  buf  f  er) ; 
if  ( i  adex(  buf  f  er  i  '  i  nvari  aat 1  )  <>  0)  and  ( i  ndex(  buf  f  er ,  *  t  empl  at  e* )  <>  0)  then  bef  i  n 
(X  extract  the  indentation  froa  the  declaration  line  X) 

1  ndent  :  «  "; 
i  :  =  l; 

while  buffer(i)  =  ’  ■  do  begin  Indent  :  =  Indent  *  '  * ;  i  :  =  i  el;  end; 

(x  extract  the  teaplate  aaan  froa  the  declaration  line  x) 
i  :  =  i  ndex(  buf  fer,  ’  teapl  ate*  )  ♦  s; 
while  bufferfi]  =  *  '  do  i  :  =  i  *  l; 

j  :  =  i  ; 

while  (j  <  I  eagtb(  buf  fer))  and  (  buffer(j]  <>  •  ’ )  do  j  :  =  j  ♦  i; 

TeapNaae  :*  substr(  bof  fer,  i  ,  J-i  ♦  1) ; 

(X  read  past  the  begin  stateaeat  X) 
readl  a  ( I  nFi  I  e,  buf  f  er) ; 

(X  read  the  first  line  of  the  teaplate  definition  X) 
readl  a  ( I  aFi  I  e,  buf  f  e r ) ; 

i  ■  »  i; 

(X  process  each  additional  line  of  the  teaplate  definition  X) 
while  i  ndex(  buf  f  er.  •  end’  )  =  0  do  begin 
(X  write  the  necessary  declarations  for  each  template  cell  X) 
write  (Pfile,lndenti  Temp  Name, '.cfg(.',j,’.).r:=’); 
if  ( i  ndex(  buf  f  er  i  '  ,  C  ) -i  odex(  buf  f  er  > '  (  r '  )  -  1 )  >  2 
then  write  (  pf  I  I  e ,  subst  r  (  buf  f  e  r ,  i  odei  (  buf  f  er ,  '  (  r  ’  )  *  2, 

i  odex<  buf  f  er ,  ’ ,  c’ )  -i  ndei(  buf  f  er ,  •  (  r *  )  -a ) ,  ’  ;  ’  ) 
else  write  (Pfile,*  0;'); 
write  (Pfile,’  >,  TeapNaae,  •.  efg<  j  ,■•).  c:  *•) ; 
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if  ( i  ndex(  buf  f  er,  ’  )  =•  ) -i  ndex(  buf  f  er .  '  ■  c»  )  )  >  2 
tt«i  writ*  (  pf  i  I  e,  aubatrf  buf  f  er,  I  odex(  buf  f  er, '  ,  c*  )  »2, 

i  udex(  buffer,  ' )  =  ’  )  -i  ndex(  buf  f  er, ' ,  c' )  -2)  ,  '  ;  *  ) 
cl  *•  writ*  (Pflle,’  o ;  •  > ; 
writ*  (Pfilc, '  ’ ,  TeapNaae,  ’ .  cf  C( .  ’  ■  a, 1  .  )  .  «R  *• ) ; 

if  i  ndeK  buf  f  «r,  ’ ;  ’ )  «  0 

th*o  writcla  (  Pf  1 1  e,  aobst  r(  bof  f  er,  i  ndex(  buf  f  er ,  '  =' )  ♦  I , 

1  eaft  h<  buf  f«r)-i  udcx(  buf  f*r,  •=•)),*;•  ) 
else  writeln  (  Pf  I  I  *,  cub*  t  r  (  buf  f  er ,  i  ndei(  buf  f  er,  *  * 1  )  ♦  1 , 

i  *dex(  buf  f  er,  ’ ;  • )  -i  ndex(  buffer,  ’  = '  )  - 1 ) ) ; 

i  :  *  i  *  x; 

recdl  n  ( I  of  i  I  e,  buf  f  er) ; 
end;  (x  while  'end'  do  X) 
writeln  (  pf  i  I  e,  I  adent ,  TeapNaan,  • .  do*  :  = 
end;  (X  then  invariant  template  definition  X) 
end;  (x  while  not  eof  X) 

(x  find  the  be<i  nni  n(  of  the  executable  portion  of  the  routine  x) 

reset  ( I  nFi  I  e) ; 

read)  n  ( I  nFi  I  e,  buf  fer) ; 

while  i  nde  x  (  buffer,  1  bed  i  n* )  =  0  do  read)  n  ( I  nf  i  I  e,  buf  f  er)  ; 

(x  read  epch  line  of  tbe  file  and  translate  accord!  nil  y  x) 
writeln  ('  .  .  .  processi  nf  the  executable  st  ateaents*  )  ; 

readl  n  (  I  nFi  I  e,  buf  fer) ; 

while  i  ndex(  buffer,  >  end.  ’  )  =  0  do  bed  a 
(x  if  the  line  doesn't  contain  any  PASCAL  contructs  or  explicit  x) 

(x  I A  procedures,  then  translate  the  XA  expression  to  PASCAL  x) 


ia  procedures, 

then  translate  the  XA  expression  to  PASCAL 

( ( i  ndex( buffer 

*  and  ' 

)  = 

0) 

and 

( i  ndex(  buffer 

1  array* 

)  = 

0) 

and 

( i  ndex(  buff  er 

*  begi  n* 

)  = 

0) 

and 

( i ndex( buffer 

’caie  * 

)  * 

0) 

and 

( i odex( buf  fer 

’const  * 

>  * 

0) 

and 

( i ndex( buffer 

•  do  * 

)  = 

0) 

and 

( i  ndex(  buf  fer 

1  down  to  1 

)  = 

0) 

and 

( i  nd*x(  buffer 

•else  * 

)  * 

0) 

and 

( i ndex( buf  fer 

»  end* 

>  * 

0) 

and 

( i ndex( buf  fer 

»  for  * 

>  * 

0) 

and 

( i ndex( buffer 

*  f  uacti on  1 

)  = 

0) 

and 

( i ndex( buffer 

•if  * 

>  * 

0) 

and 

( i  ndex(  buf  fer 

* i nput 1 

)  = 

0) 

and 

( i  nde  x ( buf  fer 

*  mod  * 

)  = 

0) 

and 

( i ndex( buf  fer 

1  not  * 

>  * 

0) 

and 

( i ndex( buffer 

'  of  » 

>  * 

0) 

and 

( i adex( buffer 

•or  * 

)  = 

0) 

and 

( i  adex(  buf  fer 

' oat  put  * 

)  * 

0) 

and 

( i ndexl buffer 

* procedure  * 

>  = 

0) 

and 

( i  ndex(  buf  f  er 

*  program  » 

>  « 

0) 

and 

( i  ndex(  buf  fer 

*  record  * 

)  = 

0) 

and 
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( i adexf buffer, 

1  atrt  a*’ 

) 

= 

0) 

aad 

( i od ext buffcr( 

■  t  he  a  ' 

) 

= 

0) 

aad 

( i adexf bufferi 

•type  • 

) 

* 

0) 

aad 

( i adex( bufferi 

•  to  • 

) 

3 

0) 

aad 

( i ndex( bufferi 

•until  • 

) 

= 

0) 

aad 

( i  BdexC bufferi 

•  var  ’ 

) 

0) 

aad 

( i  odex(  bufferi 

•  wtai  1  e  > 

) 

0) 

aad 

( i  ndex(  buffer! 

•  Mi  th  • 

) 

0) 

and 

( i  adex{ buffer! 

•  bool  eaa* 

) 

= 

0) 

aad 

( i adex( bufferi 

’  byte* 

) 

0) 

aad 

( I adex( buffert 

’  char* 

) 

= 

0) 

aad 

( i adex( bufferi 

’  cl  ose* 

) 

0) 

aad 

( ■  ndex( buf  f  er, 

*  i  ate*er* 

) 

* 

0) 

and 

( i  adexf  buffer, 

'open  * 

) 

0) 

aad 

( i  adexf  buf  f er, 

*  read* 

) 

0) 

aad 

( i  adexf  buffer. 

•  real  ' 

) 

= 

0) 

aad 

f i adexf  buffer, 

’  reset* 

) 

= 

0) 

and 

f  i  adexf  buf  f  er, 

•  rewri  te* 

) 

0) 

aad 

f  i  adexf  buffer, 

’text  * 

) 

= 

0) 

aad 

f i adexf  buffer. 

’  true’ 

) 

= 

0) 

aad 

f  i  adexf  buffer, 

'  wri  te’ 

) 

= 

0) 

aad 

( ■  adexf  buffer, 

*  I  mat*’ 

) 

0) 

aad 

( i adexf  buffer, 

*  coaf  i  (temp* 

) 

0) 

and 

f  i  adexf  buffer. 

'  tempi  ate* 

) 

= 

0)) 

then  bogi  a  (f  translate  the  line  froa  IA  to  PASCAL  t) 

(s  initialize  the  indentation  for  the  line  S) 
i  ndent  :  =  *  • ; 

i  :  =  i; 

while  buffer(i)  =  *  *  do  begin 
i  ndent  :  =  i  ndent  ♦  •  *  ; 
i  :  =  i  ♦  1; 
end; 

(*  renove  all  blanks  froa  the  buffer  t) 
oldbuffer  : =  buffer; 
buf  f  er  : =  *  * ; 

for  i  l  to  I  engt  b(  ol  dbuf  f  er )  do 
if  oldbuffer(i)  <>  '  1  then  buffer  :=  buffer  ♦  ol  dbuf  f  er  [  i  ] ; 

(t  delete  the  seat -col  on  froa  the  end  of  the  line  s) 

if  i ndei( buf f er , 1 ; • )  <>  0  then  buffer  subst r( buf f er , l , i ndei( buf f er , * ; ' ) 

(s  translate  parentheses  to  ooo-coaant  deliaiters  for  writing  to  Pfile  s) 
oldbuffer  :=  buffer; 

for  i  :=  1  to  I  e  ngt  h(  ol  dbuf  f  e  r )  do  begin 
if  oldbuffer[i]  =  *(’  then  oldbuffer(i)  :=  '(*; 
if  oldbuffer[i]  =  »)•  then  oldbufferfi]  :=  *)’; 
end;  (s  for  i  do  t) 

(s  print  t  be  IA  expression  to  PFi  I  e  *) 

writeln  (  Pf  i  I  e,  I  ndent ,  1  (  txxxx  ’  ,  ol  dbuf  f  er,  *  ttsst)'); 

(t  pad  the  string  with  one  blank  on  the  end  t) 
buf  f  er : =  buf  f  er  ♦  *  * ; 
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(X  initialize  the  expression  binary  tree  to  null  X) 
for  i  :  =  1  to  MaxNode  do  befi  n 
Expressi  on( .  i  .  ) .  val  ue  := 

Expressi  on( .  i  .).  parent  :  =  o; 

Expressi  out .  i  .).  I  cbl  I  d  :=  o; 

Expressi  oo(  .  i  .).  rcbi  I  d  :  =  o; 
end; 

(X  evaluate  t  be  expression  and  translate  the  line  froa  1A  to  PASCAL  X) 
lAParser  (  buffer ,  Expressi  on,  no  dept  r ) ; 

I  ATransI  ator  (Expressi  on,  no  deptr); 
end  (x  translation  froa  IA  to  PASCAL  X) 
else  befi  n  (X  PASCAL/1A  construct  or  leapt  ate  definition  x) 

(x  translate  parentheses  to  non-coaaent  del  inters  X) 
for  i  :=  1  to  I  en£t  b(  buf  f  er)  do  be(i  n 
if  buf  f  er(  i  }  =  •(’  then  bufferfi]  :  =  •<•; 
if  bufferfi]  =  •)’  then  buffer(i)  :  =  •)'; 
end;  (X  for  i  do  X) 
if  i  ndex(  buf  f  er  ■  ’  var  i  ant '  )  =  0 
then  befi  n  (X  PASCAL/IA  construct:  write  the  line  as  is  x) 
write  (  Pf  i  I  e,  buf  f  er ) ; 

(X  ensure  the  lines  end  with  a  unless  the  buffer  ends  X) 

(X  with  the  reserved  words  1  be(i  n’ ■  '  do*  ■  ’then1,  or  ’else1  X) 
if  ( i  ndex(  buf  f  er  ■  * ;  • )  =0 )  and 

( i  ndex(  buffer,  ’  befi  n’  )  <  >  I  enft  b(  buffer)  -  4)  and 
( i  ndex(  buffer,  •  do’ )  <>l  engt  b(  buf  f  er)  -  I )  and 
( i  ndex(  buffer,  '  then’  )  <  >1  enft  h(  buf  f  er )  -3)  and 
( i  ndex(  buffer,  •  el  se'  )  <  » I  en*t  h(  buffer)  -  3) 
then  wrl  tel  n  (Pfile,';') 
el  se  wri  tel  n  (  Pf  i  I  e) ; 
end  (X  PASCAL/IA  construct  X) 

else  be<i  n  (X  teapl  ate  definition  encountered  -  skip  it  X) 
readl  n  ( I  nFi  I  e,  buf  f  er) ; 

while  i  ndex(  buf  f  er,  *  end'  )  =  0  do  readl  n  ( I  nFi  I  e,  buf  f  er ) ; 
end;  (X  else  teswlate  definition  x) 
end;  (X  else  PASCAL/IA  construct  or  teaplate  definition  x) 
readl  n  ( I  nFi  I  e,  buf  f  e r ) ; 
end;  (i  while  not  EOF  x) 

wri  tel  n  (Pfile,'  end.  (I  profraa  ’ ,  Name,  '  X)’); 
close  (  Pf  i  I  e )  ; 

end.  (X  procedure  I  Bate  a!  te  br  a_  Pr  e  pr  oce  s  sor  x) 


b’ 
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Appendix  C:  AFIT  Input  /  Output  Operations 


[ inher  i  t (  ’ riaoper ,env'  )  ,  environment!  ’ rio . env ’  )  ] 

Module  Rea l_IA_IO_Ope rat  ions  ( i nput , output )  ; 

procedure  Getlmage  (var  A  :  Image;  FileName  :  NameStr )  ; 

(*  prompt  and  retrieve  the  requested  image  from  disk  *) 
var  i , j  :  integer ; 
row , co I  :  real ; 

InFile  :  file  of  real; 
begin  (*  procedure  Getlmage  *) 
if  FileName  =  ’’  then  begin 
write  (’read  image  from  file:  ’); 
readln  (FileName); 
end;  (*  if  FileName  *) 
i f  index ( F i leName , ’ . ’ )  =  0 
then  FileName  :=  FileName  +  ’ . img ’ ; 

open (InFile, FileName , unknown ) ; 
reset  ( InF i le )  ; 

(*  read  in  the  input  array  *) 
read  ( InFile , row, col ) ; 

A . row  :=  round) row); 

A. col  :=  round(col); 
for  i  : =  1  to  A. row  do 

for  j  : =  1  to  A. col  do  read  ( InFile , A . gray (. i , j .)) ; 
close  (InFile); 

end;  (*  procedure  Getlmage  *) 

procedure  Putlmage  (A  :  Image;  FileName  :  NameStr) ; 

(*  write  the  image  to  an  external  file  *) 
var  i , j  :  integer ; 

OutFile  :  file  of  real; 
begin  (*  procedure  Putlmage  *) 
if  i ndex ( F i leName )  =  0 
then  FileName  :  =  FileName  +  ’.img’; 

open ( OutFile .FileName , unknown ) ; 
rewrite  (OutFile); 
write  (OutFile , A. row, A.col ) ; 
for  i  : =  1  to  A. row  do 

for  j  : =  1  to  A.col  do  write  ( OutF i le , A . gray ( . i , j . ) ) 
close  ( OutF i Le ) ; 
end;  (*  procedure  Putlmage  *) 
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procedure  GetTemplate  (var  A  :  Template;  FileName  :  NameStr) 
(*  prompt  and  retrieve  the  requested  template  from  disk  *) 
var  cnt  :  integer; 

rowof f , colo f f , we ight  :  real; 

InFile  :  file  of  real; 
begin  (*  procedure  GetTemplate  *) 
if  FileName  =  ’’  then  begin 
write  (’read  template  from  file:  ’); 
readln  (FileName); 
end;  (*  if  FileName  *) 
if  index ( F ' leName )  =  0 
then  FileName  :=  FileName  +  ’.tmp’; 
open ( InF i le , F i leName , unknown ) ; 
reset  ( InFi le )  ; 

( *  read  in  the  input  array  * ) 
cnt  : =  0 ; 

read  ( InFile, rowof f,coloff,weight)  ; 
while  ( EOF ( InFi le ) =false )  do  begin 
cnt  : =  cnt  +  1  ; 

A . cf g (  . cnt .  ) . r  :  =  round ( rowof f )  ; 

A . cf g (  . cnt . ) . c  :  =  round ( colof f  )  ; 

A . cfg ( . cnt . ) . w  :=  weight; 

read  ( InFile, rowof f , colof f , weight )  ; 

end;  (*  while  not  EOF  *) 

A . num  : =  cnt ; 
close  ( InFile ) ; 

end;  (*  procedure  GetTemplate  *) 

procedure  PutTemplate  (A  :  Template;  FileName  :  NameStr); 

(*  write  the  template  to  an  external  file  *) 
var  i  :  integer; 

OutFile  :  file  of  real; 
begin  (*  procedure  PutTemplate  *) 
if  i ndex ( Fi leName ,'.’ )  =  0 
then  FileName  :=  FileName  +  ’.tmp’; 

open ( OutF i le , Fi leName , unknown ) ; 
rewrite  (OutFile); 
for  i  : =  l  to  A . num  do 

write  (OutFile, A. cfg(  .i.  ) . r , A . c f g (  .i.  ) . c , A . c f g (  .i  .  )  . w )  ; 
close  (OutFile); 

end;  (*  procedure  PutTemplate  *) 
end.  (*  Module  Real  IA  10  Operations  *) 


Appendix  D:  Automated  Image  Algebra  Translator 


Translate  Image  Algebra 

function  -  This  command  procedure  translates  an  image 
processing  function  from  its  description  in 
the  image  algebra  to  an  executable  file. 

format  -  @TIA  [filename] 

on  error  then  exit 

if  pi  .nes.  ""  then  goto  preprocess 

inquire  pi  "image  algebra  routine  to  be  translated 
preprocess : 

copy  ’pl'.ia  translat.ia 

write  sys$output  "...preprocessing" 

r  preproc 

write  sys$output  "...compiling" 
pas  translat.pas 

write  sys$output  "...linking" 
link  translat , iaoper , io 

copy  translat.exe  ’pl’.exe 

write  sys$output  "...executable  file  built" 

del  translat.*;* 
purge 
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